Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,188,411 C→T 5.8% E2K (GAA→AAA)  yehA ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,188,4110CT 5.8% 220.0 / 2.4 69E2K (GAA→AAA) yehAputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base C (30/35);  new base T (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

TTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAAC  >  NC_000913/2188263‑2188556
                                                                                                                                                    |                                                                                                                                                 
ttATGCTCATATGTAGCCATTTGTGTTTGAGTAATTGCATAATCCCCCTGAACAACCTAGTTGTGCTGAACAGATAACTGACTATAAGATTCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTc                                                                                                                                                   <  2:49155/149‑1 (MQ=255)
   tGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTTCTTAACAGATAACTAAATATAACATGCAGCATAGCTCACAGGCATAACGATCATCATTAATATAAATAGCATTATGCGAATTTCCa                                                                                                                                                 <  1:217773/148‑1 (MQ=255)
       cATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGaa                                                                                                                                            >  1:67242/1‑149 (MQ=255)
         tatGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATATCCATCATCATTAATATAAATAGAATTATGCTAATTTCCATAGacc                                                                                                                                           >  6:167179/1‑146 (MQ=255)
           tGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAaccac                                                                                                                                        <  7:99292/149‑1 (MQ=255)
                  aTTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTa                                                                                                                                 >  5:47834/1‑149 (MQ=255)
                       tgtTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACt                                                                                                                            <  7:2603/149‑1 (MQ=255)
                       tgtTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACt                                                                                                                            <  2:75349/149‑1 (MQ=255)
                           tGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGc                                                                                                                         <  5:130065/148‑1 (MQ=255)
                             aGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCaaa                                                                                                                      >  7:229744/1‑149 (MQ=255)
                                   tGCAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAAt                                                                                                                >  7:146847/1‑149 (MQ=255)
                                   tGCAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAAt                                                                                                                >  5:263887/1‑149 (MQ=255)
                                     cAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAAttt                                                                                                              <  1:295856/149‑1 (MQ=255)
                                     cAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAAtt                                                                                                               <  3:224435/148‑1 (MQ=255)
                                        aaTCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTAAAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAacca                                                                                                                                         <  8:128533/119‑1 (MQ=255)
                                        aaTCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTAAAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAacca                                                                                                                                         >  7:128533/1‑119 (MQ=255)
                                          tCCACCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATACCCATAATCATTAATATAAATAGCATTATGCGAATTTCAATAGAACCACCTGACTATTACTGGCAAATATAATTTATACt                                                                                                         >  5:165412/1‑149 (MQ=255)
                                           cccccTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACtt                                                                                                        >  4:310192/1‑149 (MQ=255)
                                             cccTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCACCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACttt                                                                                                       >  8:273603/1‑148 (MQ=255)
                                              ccTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCTGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTtca                                                                                                     <  3:49051/149‑1 (MQ=255)
                                                tGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCAGCTGACTATTACTGGCAAATATAATTTCTACTTtcat                                                                                                    >  6:74075/1‑148 (MQ=255)
                                                tGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTtcat                                                                                                    >  4:51031/1‑148 (MQ=255)
                                                 gAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTtcatc                                                                                                   >  8:140092/1‑148 (MQ=255)
                                                        cAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCAAAGCTCACAGGCATAACCATCATCATAAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGGAAATATAATTTCTACTTTCATCAATCTCtt                                                                                           <  3:61744/149‑1 (MQ=255)
                                                         aaGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCttt                                                                                          <  7:89403/149‑1 (MQ=255)
                                                         aaGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCttt                                                                                          <  6:47834/149‑1 (MQ=255)
                                                           gttgttCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTaa                                                                                        <  5:233574/149‑2 (MQ=255)
                                                               tgCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGaa                                                                                    <  8:184549/149‑1 (MQ=255)
                                                                 cTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAgg                                                                                  <  8:248924/149‑1 (MQ=255)
                                                                  tGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGt                                                                                 >  2:262449/1‑149 (MQ=255)
                                                                  tGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGt                                                                                 >  1:294980/1‑149 (MQ=255)
                                                                       aGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATtg                                                                            >  5:15543/1‑149 (MQ=255)
                                                                           aaCTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCa                                                                        >  5:158950/1‑149 (MQ=255)
                                                                            aCTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCa                                                                        <  6:158950/148‑1 (MQ=255)
                                                                            aCTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAAATTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCa                                                                        <  1:28692/148‑1 (MQ=255)
                                                                                aCTATAACATGCCGCCTAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATCCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGaaa                                                                   <  4:142073/149‑1 (MQ=255)
                                                                                aCTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGaaa                                                                   <  5:167179/149‑1 (MQ=255)
                                                                                aCTATAACATGCCGCATAGCTCACAGACATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATCGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGaaa                                                                   <  1:197576/149‑1 (MQ=255)
                                                                                  tataACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGt                                                                 >  1:129275/1‑149 (MQ=255)
                                                                                    taACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTcc                                                               <  6:263887/149‑1 (MQ=255)
                                                                                    taACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTcc                                                               <  2:294980/149‑1 (MQ=255)
                                                                                    taACATGCCGCATAGCTCACAGGAATAACCATAATAATTAATATAAATATCATTATGCTAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTcc                                                               <  6:165412/149‑1 (MQ=255)
                                                                                              cATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTa                                                     <  1:304020/149‑1 (MQ=255)
                                                                                              cATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTa                                                     <  3:51031/149‑1 (MQ=255)
                                                                                                tAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTAt                                                    >  1:80883/1‑148 (MQ=255)
                                                                                                   cTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTTCATAGAACCCCCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTACTTGTTGCTTATCAAt                                                >  1:331274/1‑149 (MQ=255)
                                                                                                   cTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTTCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAAt                                                >  1:127550/1‑149 (MQ=255)
                                                                                                      acaGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGc                                             >  8:291092/1‑149 (MQ=255)
                                                                                                      acaGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGc                                             <  2:67242/149‑1 (MQ=255)
                                                                                                      acaGGCATAACCATCACCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAAtt                                               <  7:273603/147‑1 (MQ=255)
                                                                                                       caGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGcc                                            <  6:15543/149‑1 (MQ=255)
                                                                                                        aGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAgg                                                                                  <  1:157837/110‑1 (MQ=255)
                                                                                                        aGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAgg                                                                                  >  2:157837/1‑110 (MQ=255)
                                                                                                             tAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACAtt                                       >  2:207309/1‑148 (MQ=255)
                                                                                                                     atcatTAATATAAATAGCATTATGCGAATTTTCATATAAACACCTGACTATTACTGGCATATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCATTTGCCACATTAACCGATgg                              <  2:331274/149‑1 (MQ=255)
                                                                                                                     atcatTAATATAAATAGCATTATGCGAATTTTCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATgg                              <  2:127550/149‑1 (MQ=255)
                                                                                                                     atcatTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATg                               <  8:146847/148‑1 (MQ=255)
                                                                                                                      tcatTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGc                             >  7:61338/1‑149 (MQ=255)
                                                                                                                       catTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATgg                              >  8:264302/1‑147 (MQ=255)
                                                                                                                       catTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATgg                              <  7:264302/147‑1 (MQ=255)
                                                                                                                        atTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCgg                           >  3:270665/1‑149 (MQ=255)
                                                                                                                        atTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCgg                           <  7:140092/149‑1 (MQ=255)
                                                                                                                              ataAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACtt                      >  7:109970/1‑148 (MQ=255)
                                                                                                                                   tAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAAc                                    <  7:261132/129‑1 (MQ=255)
                                                                                                                                   tAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAAc                                    >  8:261132/1‑129 (MQ=255)
                                                                                                                                      cATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACTTCATTGGTGc             <  5:200239/149‑1 (MQ=255)
                                                                                                                                       aTTATGCGAATTTCAATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACTTCATTGGTGc             >  6:198213/1‑148 (MQ=255)
                                                                                                                                        ttATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACTTCATTGGTGCg            >  3:333282/1‑148 (MQ=255)
                                                                                                                                                 tttCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAAc  >  7:325045/1‑149 (MQ=255)
                                                                                                                                                 tttCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAAc  <  8:229744/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTATGCTCATATGTCGCCATTTGTGTTTGAGTAAGTGCAAAATCCCCCTGAACAACCAAGTTGTGCTGAACAGATAACTCACTATAACATGCCGCATAGCTCACAGGCATAACCATCATCATTAATATAAATAGCATTATGCGAATTTCCATAGAACCACCTGACTATTACTGGCAAATATAATTTCTACTTTCATCAATCTCTTTACCGAAGGTGATTGTGCATGAAAGTCCTTGTTGCTTATCAATTGCCACATTAACCGATGGCGGTACTTCATTGGTGCGAATAAATAAC  >  NC_000913/2188263‑2188556

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