Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,105,716 G→A 6.6% P815L (CCG→CTG)  ydiJ ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,105,7160GA 6.6% 427.0 / 13.5 147P815L (CCG→CTG) ydiJputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (62/74);  new base A (6/4);  total (68/79)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 147 total aligned reads displayed.

TACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGCCGGCA  >  NC_002947/5105569‑5105863
                                                                                                                                                   |                                                                                                                                                     
 acacCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGgcg                                                                                                                                                   <  1:1868273/150‑1 (MQ=255)
   acCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGcc                                                                                                                                                 >  2:1817568/1‑150 (MQ=255)
     cAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAg                                                                                                                                               <  2:951853/150‑1 (MQ=255)
      aGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAACCCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACGCCGCCAgg                                                                                                                                              >  1:1015068/1‑150 (MQ=255)
          tCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATa                                                                                                                                          >  2:1744623/1‑150 (MQ=255)
          tCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATa                                                                                                                                          >  1:1531603/1‑150 (MQ=255)
          tCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATa                                                                                                                                          >  2:1790409/1‑150 (MQ=255)
          tCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATa                                                                                                                                          >  2:1047548/1‑150 (MQ=255)
             tCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGa                                                                                                                                             <  2:1798701/144‑1 (MQ=255)
             tCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGa                                                                                                                                             >  1:1798701/1‑144 (MQ=255)
             tCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaa                                                                                                                                               >  1:953724/1‑141 (MQ=255)
             tCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaa                                                                                                                                               <  2:953724/142‑2 (MQ=255)
               gCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCgg                                                                                                                                     <  1:1614843/150‑1 (MQ=255)
               gCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCgg                                                                                                                                     <  2:1363558/150‑1 (MQ=255)
                 cGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTg                                                                                                                                   <  2:44886/150‑1 (MQ=255)
                     tCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGccc                                                                                                                               >  2:963687/1‑150 (MQ=255)
                     tCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGccc                                                                                                                               >  1:2167651/1‑150 (MQ=255)
                           cccACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCgg                                                                                                                         <  2:86884/150‑1 (MQ=255)
                            ccaccaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGt                                                                                                                        <  2:1691032/150‑1 (MQ=255)
                            ccaccaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaacgg                                                                                                                                            <  1:1710492/130‑5 (MQ=255)
                             caccaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGt                                                                                                                        >  2:1329183/1‑149 (MQ=255)
                             caccaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGt                                                                                                                        <  1:1329183/149‑1 (MQ=255)
                             caccaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACCGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTg                                                                                                                       <  1:1978780/150‑1 (MQ=255)
                                caGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGc                                                                                                                    <  2:1633759/150‑1 (MQ=255)
                                caGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCCGTGGCCCAGGCGGTGGGc                                                                                                                    <  2:1006908/150‑1 (MQ=255)
                                       acacCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAg                                                                                                                 <  2:1489352/146‑1 (MQ=255)
                                        cacCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGa                                                                                                            <  1:2357125/150‑1 (MQ=255)
                                         acCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAg                                                                                                                                               <  1:1867400/114‑1 (MQ=255)
                                         acCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAg                                                                                                                                               >  2:1867400/1‑114 (MQ=255)
                                         acCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGAt                                                                                                           <  1:1946859/150‑1 (MQ=255)
                                             cAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAg                                                                                                       >  2:645259/1‑150 (MQ=255)
                                              aaTCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGGTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACAGCTTGCCGTTGGCGCTGTCCGGCGCCAGGAATACCCGGACTGCCATGCTGGATGCAAGGTCGATGAAgg                                                                                                      >  2:872969/1‑150 (MQ=255)
                                                tCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGaa                                                                                                        >  2:552247/1‑146 (MQ=255)
                                                  ggCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGa                                                                                                  <  1:1790409/150‑1 (MQ=255)
                                                   gCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaac                                                                                                                                              >  2:693916/1‑103 (MQ=255)
                                                   gCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGAAAAGGAAACCCTGTACATAAAGTAATGTGCCGTTGGCGCTGTCGCCAGTCAGGAGTATCTGTTGGCTCAGGCGGTGGGAAAGGTCGATGAACGCCGAc                                                                                                 >  1:541934/1‑150 (MQ=255)
                                                    ccAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGAca                                                                                                <  2:1083397/150‑1 (MQ=255)
                                                    ccAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGAca                                                                                                >  2:1244939/1‑150 (MQ=255)
                                                     cAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGAcag                                                                                               <  2:2167651/150‑1 (MQ=255)
                                                     cAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaac                                                                                                                                              >  1:1372676/1‑101 (MQ=255)
                                                      aGGGCCTTGAGCTGCGTGGCGTTGCGGATAGCTGCCTTGGCAACAGCGCCACGGCAAACATGTAAGTGCACCGGCTCACCGTTGGCGGTCTACGGCGAGACAGATACCACGTGCGCCCGGCGATGGGAAAGGTAGCAGAAGGCCGAcagc                                                                                              >  2:402727/1‑150 (MQ=255)
                                                           cTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAgcgg                                                                                         <  1:1744623/150‑1 (MQ=255)
                                                               aGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGtt                                                                                     >  1:2397232/1‑150 (MQ=255)
                                                                 cTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTc                                                                                   <  1:157952/150‑1 (MQ=255)
                                                                        gCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAg                                                                            <  1:1034895/150‑1 (MQ=255)
                                                                         cGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGc                                                                           <  2:1531603/150‑1 (MQ=255)
                                                                          gTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCg                                                                          >  2:2281194/1‑150 (MQ=255)
                                                                          gTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCg                                                                          >  1:2152861/1‑150 (MQ=255)
                                                                           ttGCGGATCGCTGCCTTGGCAAACGCGCCAAGTAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCgg                                                                         <  2:306805/150‑1 (MQ=255)
                                                                               ggATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTgg                                                                                                                      >  1:475072/1‑101 (MQ=255)
                                                                               ggATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTgg                                                                                                                      <  2:475072/101‑1 (MQ=255)
                                                                               ggATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGAt                                                                                                           >  2:2234292/1‑112 (MQ=255)
                                                                               ggATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGAt                                                                                                           <  1:2234292/112‑1 (MQ=255)
                                                                                gATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGc                                                                           >  2:2365528/1‑143 (MQ=255)
                                                                                gATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGaa                                                                   <  1:525713/150‑1 (MQ=255)
                                                                                 ttcggtGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAg                                                                  >  2:1954707/6‑150 (MQ=255)
                                                                                      tGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGgcgc                                                                                                                                                  <  2:1342356/66‑1 (MQ=255)
                                                                                      tGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTc                                                             >  1:1333435/1‑150 (MQ=255)
                                                                                       gCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTcc                                                            >  2:995598/1‑150 (MQ=255)
                                                                                        ccTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCgg                                                                         >  2:1652515/1‑137 (MQ=255)
                                                                                        ccTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCgg                                                                         <  1:1652515/137‑1 (MQ=255)
                                                                                          ttGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTg                                                                                                                                   <  1:995600/77‑1 (MQ=255)
                                                                                          ttGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAg                                                                                                       <  1:2290790/105‑1 (MQ=255)
                                                                                          ttGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGAGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGGAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGc                                                         >  1:739946/1‑150 (MQ=255)
                                                                                           tGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCa                                                        >  1:1221421/1‑150 (MQ=255)
                                                                                               aaaCGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTGGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAg                                                    >  1:320435/1‑150 (MQ=255)
                                                                                               aaaCGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGAca                                                                                                <  1:1126020/107‑1 (MQ=255)
                                                                                               aaaCGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGAca                                                                                                >  2:1126020/1‑107 (MQ=255)
                                                                                                aaCGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGc                                                   <  2:1221421/150‑1 (MQ=255)
                                                                                                  cgcgCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCAc                                                 <  1:2281194/150‑1 (MQ=255)
                                                                                                    cgcCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGa                                               <  1:2329735/150‑1 (MQ=255)
                                                                                                        aaGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAg                                                                                                                             >  2:2084815/1‑69 (MQ=255)
                                                                                                        aaGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAg                                                                                                                             >  2:2085211/1‑69 (MQ=255)
                                                                                                        aaGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAg                                                                                                                             <  1:2085211/69‑1 (MQ=255)
                                                                                                        aaGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAg                                                                                                                             <  1:2084815/69‑1 (MQ=255)
                                                                                                         aGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCgg                                                                         >  1:1521449/1‑120 (MQ=255)
                                                                                                            aaaCCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCg                                                                          >  1:664220/1‑116 (MQ=255)
                                                                                                            aaaCCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCg                                                                          <  2:664220/116‑1 (MQ=255)
                                                                                                              aCCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCAATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTc                                     >  1:855509/1‑150 (MQ=255)
                                                                                                                ccTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTcgc                                   >  2:1987452/1‑150 (MQ=255)
                                                                                                                       gTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGAt                                                                                                           <  2:1742241/72‑1 (MQ=255)
                                                                                                                        tGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGaa                                                                                >  1:1362167/1‑98 (MQ=255)
                                                                                                                        tGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGaa                                                                                <  2:1362167/98‑1 (MQ=255)
                                                                                                                         ccAACGGCAAGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCg                          >  1:1190304/2‑150 (MQ=255)
                                                                                                                          cAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACGGGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCgg                         >  1:1508693/1‑150 (MQ=255)
                                                                                                                          cAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCgg                         >  2:2127896/1‑150 (MQ=255)
                                                                                                                           aaCGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGNTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGATGGGCCggg                        <  1:402727/150‑1 (MQ=255)
                                                                                                                           aaCGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCggg                        <  1:645259/150‑1 (MQ=255)
                                                                                                                           aaCGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCggg                        <  2:1257383/150‑1 (MQ=255)
                                                                                                                            aCGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCggg                        >  1:1509691/1‑149 (MQ=255)
                                                                                                                            aCGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCggg                        <  2:1509691/149‑1 (MQ=255)
                                                                                                                             cGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAAAAGCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGt                      >  1:1770652/1‑150 (MQ=255)
                                                                                                                              ggCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTc                     <  1:1018353/150‑1 (MQ=255)
                                                                                                                                   gCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTc                <  1:1244939/150‑1 (MQ=255)
                                                                                                                                     cGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGTACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTcgc              >  1:1707101/1‑150 (MQ=255)
                                                                                                                                      gTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGTAAAGTGCTGGCATCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTcgcg             >  2:1270076/1‑150 (MQ=255)
                                                                                                                                        tGGCGCCGTACAGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCa           >  2:407610/1‑150 (MQ=255)
                                                                                                                                             cTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGcc      <  1:962615/150‑1 (MQ=255)
                                                                                                                                              tGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAg                                                                  <  2:1572744/89‑1 (MQ=255)
                                                                                                                                                 aCGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTACACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGCCGGCa  >  2:2356548/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGCCGGCA  >  NC_002947/5105569‑5105863

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