Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,545,568 T→C 7.6% L158P (CTG→CCG)  fimD → fimbrial usher outer membrane porin protein; FimCD chaperone‑usher

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,545,5680TC 7.6% 238.0 / 8.1 79L158P (CTG→CCG) fimDfimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Reads supporting (aligned to +/- strand):  ref base T (37/36);  new base C (3/3);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

CTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGC  >  NC_000913/4545420‑4545715
                                                                                                                                                    |                                                                                                                                                   
ctgACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCt                                                                                                                                                     >  5:364532/1‑149 (MQ=255)
ctgACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCt                                                                                                                                                     <  1:91290/149‑1 (MQ=255)
  gACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGa                                                                                                                                                   >  2:368035/1‑149 (MQ=255)
   acacGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGa                                                                                                                                                   <  1:370235/148‑1 (MQ=255)
          cgcAACTCGCCAGTATGGGGCTGAATACGGCTTCTGGCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCt                                                                                                                                           <  8:404540/149‑1 (MQ=255)
            cAACTCGCCAGTATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCa                                                                                                                                         <  1:1739/149‑1 (MQ=255)
             aaCTCGCCAGTATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAg                                                                                                                                        >  5:227035/1‑149 (MQ=255)
               cTCGCCAGTATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGc                                                                                                                                      <  3:390687/149‑1 (MQ=255)
                    cAGTATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTACCCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTa                                                                                                                                 >  2:263910/1‑149 (MQ=255)
                       tATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTTGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGa                                                                                                                              <  7:271666/149‑1 (MQ=255)
                       tATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGa                                                                                                                              <  3:253621/149‑1 (MQ=255)
                       tATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGa                                                                                                                              <  7:68914/149‑1 (MQ=255)
                         tGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGt                                                                                                                            <  1:368035/149‑1 (MQ=255)
                               tGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATcgc                                                                                                                      >  7:422245/1‑149 (MQ=255)
                                 aaTACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATcgcgc                                                                                                                    >  8:6267/1‑149 (MQ=255)
                                 aaTACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATcgc                                                                                                                      <  8:283502/147‑1 (MQ=255)
                                 aaTACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGACATTAACCACAATGTTCCAGGACTCTACTTCGAATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATcgcgc                                                                                                                    <  4:129150/149‑1 (MQ=255)
                                   tACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATcgcgcgc                                                                                                                  >  7:341392/1‑149 (MQ=255)
                                      ggCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTg                                                                                                               >  7:151545/1‑149 (MQ=255)
                                       gctgctGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTgg                                                                                                              <  3:341874/149‑1 (MQ=255)
                                           cTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTta                                                                                                           >  2:187588/1‑148 (MQ=255)
                                           cTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGAGGTta                                                                                                           >  2:181687/1‑148 (MQ=255)
                                            tGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTtata                                                                                                         >  1:389246/1‑149 (MQ=255)
                                                 ccGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATtcct                                                                                                    >  5:372068/1‑149 (MQ=255)
                                                  cGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATtcctc                                                                                                   >  3:328873/1‑149 (MQ=255)
                                                      aTGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTcctcct                                                                                                 >  8:5706/1‑147 (MQ=255)
                                                       tGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAg                                                                                              >  5:217535/1‑149 (MQ=255)
                                                        gAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCACTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGt                                                                                             >  1:407134/1‑149 (MQ=255)
                                                               ctgGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGTTGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTAt                                                                                          <  6:364532/145‑1 (MQ=255)
                                                               ctgGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATggg                                                                                       <  2:389246/148‑1 (MQ=255)
                                                                tgGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGAt                                                                                     >  7:40769/1‑149 (MQ=255)
                                                                    ggATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCg                                                                                 <  2:180123/149‑1 (MQ=255)
                                                                     gatgatGCCTGTGTGCCATTATCCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTcctcct                                                                                                 <  3:228624/132‑1 (MQ=255)
                                                                     gatgatGACTGTGTGCCATTATCCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTcctcct                                                                                                 >  4:228624/1‑132 (MQ=255)
                                                                        gatgCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGa                                                                                                                              <  4:286164/100‑1 (MQ=255)
                                                                        gatgCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGa                                                                                                                              >  3:286164/1‑100 (MQ=255)
                                                                        gatgCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTa                                                                              <  6:217535/148‑1 (MQ=255)
                                                                        gatgCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTAt                                                                             >  4:299306/1‑149 (MQ=255)
                                                                            ccTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTaa                                                                          <  6:227035/148‑1 (MQ=255)
                                                                                    ccATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCCGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGa                                                                                               <  7:279943/119‑1 (MQ=255)
                                                                                    ccATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCCGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGa                                                                                               <  7:273269/119‑1 (MQ=255)
                                                                                    ccATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCCGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGa                                                                                               >  2:215160/1‑119 (MQ=255)
                                                                                    ccATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCCGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGa                                                                                               >  8:279943/1‑119 (MQ=255)
                                                                                    ccATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCCGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGa                                                                                               <  1:215160/119‑1 (MQ=255)
                                                                                    ccATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCCGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGa                                                                                               >  8:273269/1‑119 (MQ=255)
                                                                                     cATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATc                                                                                                                                              >  1:290323/1‑71 (MQ=255)
                                                                                     cATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATc                                                                                                                                              <  2:290323/71‑1 (MQ=255)
                                                                                       ttAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGGGTGGTTAAATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCt                                                              <  1:181687/149‑1 (MQ=255)
                                                                                       ttAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTAAATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCt                                                              <  1:187588/149‑1 (MQ=255)
                                                                                          aCCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCaa                                                           >  5:403207/1‑149 (MQ=255)
                                                                                            cacaATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTtata                                                                                                         >  4:33929/1‑101 (MQ=255)
                                                                                            cacaATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTtata                                                                                                         <  3:33929/101‑1 (MQ=255)
                                                                                            cacaATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGAGCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGAGCCCGGTATTAATGCCGGATTGCTCAAt                                                          >  6:124415/1‑148 (MQ=255)
                                                                                             acaATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCAttt                                                                                                                                  <  3:190127/75‑1 (MQ=255)
                                                                                             acaATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCAttt                                                                                                                                  >  4:190127/1‑75 (MQ=255)
                                                                                               aaTGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATTATTAATCGCGCGCGTGGTTATATTCCGACTGAGGTATGGGATCCCGGTATTAATGCCGGATTGCTCATTtat                                                       >  4:303727/1‑148 (MQ=255)
                                                                                               aaTGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATtata                                                      <  5:412095/149‑1 (MQ=255)
                                                                                                         gACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCgg                                            <  3:299306/149‑1 (MQ=255)
                                                                                                            gCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCTTTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGCGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGaaa                                         >  3:106216/1‑149 (MQ=255)
                                                                                                                cTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATa                                       <  8:40769/147‑1 (MQ=255)
                                                                                                                cTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAgt                                     <  6:403207/149‑1 (MQ=255)
                                                                                                                 tGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAgtg                                    >  8:93574/1‑149 (MQ=255)
                                                                                                                  gcgcATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATa                                       >  7:312003/1‑145 (MQ=255)
                                                                                                                  gcgcATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATa                                       <  8:312003/145‑1 (MQ=255)
                                                                                                                   cgcATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTa                                  >  1:78508/1‑149 (MQ=255)
                                                                                                                    gcATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTa                                  >  5:127577/1‑148 (MQ=255)
                                                                                                                     cATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACa                                <  3:66641/149‑1 (MQ=255)
                                                                                                                        cTGGATGTTGGTCAGCAGCGTCACAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCTTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGa                              <  8:417834/148‑1 (MQ=255)
                                                                                                                          ggATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATg                                                                        <  4:184826/104‑1 (MQ=255)
                                                                                                                          ggATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATc                           >  3:33600/1‑149 (MQ=255)
                                                                                                                          ggATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATc                           >  8:45606/1‑149 (MQ=255)
                                                                                                                           gATGTTGGTCAGCAGCGACTGAACCTGAGGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAATTAGTGTACAGAATCg                          >  5:287819/1‑149 (MQ=255)
                                                                                                                             tGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGTAAATAGTGTACAGAATCGGa                        <  1:263910/149‑1 (MQ=255)
                                                                                                                               ttGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAAt                                        <  3:70837/131‑1 (MQ=255)
                                                                                                                               ttGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAAt                                        >  4:70837/1‑131 (MQ=255)
                                                                                                                                  gTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTggg                   <  7:6267/149‑1 (MQ=255)
                                                                                                                                   tCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTgggg                  <  6:127577/149‑1 (MQ=255)
                                                                                                                                                 aCCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTAt    >  8:275936/1‑149 (MQ=255)
                                                                                                                                                  ccTGACGAACCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGAATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATg   >  1:373045/1‑149 (MQ=255)
                                                                                                                                                   cTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATg   >  2:31440/1‑148 (MQ=255)
                                                                                                                                                   cTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGc  <  7:93574/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTGACACGCGCGCAACTCGCCAGTATGGGGCTGAATACGGCTTCTGTCGCCGGTATGAATCTGCTGGCGGATGATGCCTGTGTGCCATTAACCACAATGGTCCAGGACGCTACTGCGCATCTGGATGTTGGTCAGCAGCGACTGAACCTGACGATCCCTCAGGCATTTATGAGTAATCGCGCGCGTGGTTATATTCCTCCTGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGC  >  NC_000913/4545420‑4545715

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑