Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,698,852 Δ1 bp 5.9% coding (558/1014 nt) cheBC → chemotaxis response regulator protein‑glutamate methylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,698,8520G. 5.9% 427.7 / 20.5 135coding (558/1014 nt)cheBCchemotaxis response regulator protein‑glutamate methylesterase
Reads supporting (aligned to +/- strand):  ref base G (63/64);  new base . (4/4);  total (67/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

Only 100 of 139 total aligned reads displayed.

GCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAGAA  >  NC_002947/1698703‑1699001
                                                                                                                                                     |                                                                                                                                                      
gCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAg                                                                                                                                                        <  2:1085023/150‑1 (MQ=255)
   gCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCAt                                                                                                                                                     >  1:1019358/1‑150 (MQ=255)
   gCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCAt                                                                                                                                                     <  2:2129520/150‑1 (MQ=255)
      aGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGACCAGCacgat                                                                                                                                                  <  1:1527003/150‑5 (MQ=255)
       ggTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCAtgtg                                                                                                                                                 >  1:1578070/1‑149 (MQ=255)
       ggTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCAtgtg                                                                                                                                                 <  2:1578070/149‑1 (MQ=255)
         tCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGAc                                                                                                                                              <  1:60438/150‑1 (MQ=255)
         tCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGAc                                                                                                                                              <  2:133237/150‑1 (MQ=255)
         tCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGAc                                                                                                                                              <  2:1385560/150‑1 (MQ=255)
           ggTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCa                                                                                                                                            <  2:2163981/150‑1 (MQ=255)
            gTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGGATGTGGACCAg                                                                                                                                           >  2:988752/1‑150 (MQ=255)
            gTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAg                                                                                                                                           <  2:1607716/150‑1 (MQ=255)
            gTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAg                                                                                                                                           <  2:1045788/150‑1 (MQ=255)
             tGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCAgc                                                                                                                                          <  1:821800/150‑2 (MQ=255)
              ggCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGgt                                                                                                                                         <  1:1650502/150‑1 (MQ=255)
              ggCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATAGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGgt                                                                                                                                         <  1:286091/150‑1 (MQ=255)
                  gCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGttt                                                                                                                                     <  2:1506777/150‑1 (MQ=255)
                  gCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGttt                                                                                                                                     <  1:741654/150‑1 (MQ=255)
                  gCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGttt                                                                                                                                     <  1:2259490/150‑1 (MQ=255)
                    aCCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGGCCAGCATGTGGACCAGGTGTTTGc                                                                                                                                   >  1:1949780/1‑150 (MQ=255)
                       gcggcgcGAGGCTTCCCAGCTCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGc                                                                                                                                <  1:734652/147‑1 (MQ=255)
                       gCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGc                                                                                                                                >  1:2317062/1‑150 (MQ=255)
                            gcgAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGGGGCCCGGCAGCACTCGAAGTACTGCTCAAGGGTTTGCCTGCAGCGCTTCCGGCAGCCATCGTGCTCGTCCAGCATGTGGACCAGGCGTTTGCGGCAgtag                                                                                                                           >  2:1048951/1‑147 (MQ=255)
                              gAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATg                                                                                                                         >  2:915098/1‑150 (MQ=255)
                              gAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATg                                                                                                                         >  1:333563/1‑150 (MQ=255)
                                     cctcgcgcAGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGCGTTTGCCGCAGGCATGGCCGAAt                                                                                                                  <  1:58535/146‑1 (MQ=255)
                                     cccAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAAt                                                                                                                  <  1:206722/150‑1 (MQ=255)
                                      ccAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATg                                                                                                                 >  2:1521292/1‑150 (MQ=255)
                                       cAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCAgca                                                                                                                                         >  2:1735959/1‑123 (MQ=255)
                                       cAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCAgca                                                                                                                                         <  1:1735959/125‑3 (MQ=255)
                                        aGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGc                                                                                                               >  2:564043/1‑150 (MQ=255)
                                           gccgcGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCAtgt                                                                                                                                                  <  1:996870/112‑1 (MQ=255)
                                           gccgcGGGTTCCGGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCGCAAGGGGTTGCCGGCAACGTTTCCGGCGGCCATCGTGCTGGTCCAGCAACTGGACCAGGTGTTTTCCGCAGGCATGGCCGAATGGCTca                                                                                                            <  1:2175350/150‑1 (MQ=255)
                                                     ggtggCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAgcagcgcagc                                                                                                  <  2:1363289/150‑1 (MQ=255)
                                                      gtggCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGgt                                                                                                                                         <  2:1256233/110‑1 (MQ=255)
                                                      gtggCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGgt                                                                                                                                         >  1:48220/1‑110 (MQ=255)
                                                      gtggCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGgt                                                                                                                                         <  2:48220/110‑1 (MQ=255)
                                                      gtggCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCgcagca                                                                                                 >  2:834153/1‑150 (MQ=255)
                                                           aaTCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGCTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCt                                                                                            >  2:436995/1‑150 (MQ=255)
                                                                 gTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGt                                                                                      <  1:834153/150‑1 (MQ=255)
                                                                  tCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTg                                                                                     >  1:964306/1‑150 (MQ=255)
                                                                        gcggGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTg                                                                               >  1:229491/1‑150 (MQ=255)
                                                                         cggGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGc                                                                                                                       >  2:1960596/1‑109 (MQ=255)
                                                                         cggGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGc                                                                                                                       <  1:1960596/109‑1 (MQ=255)
                                                                            gTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGccc                                                                           >  1:508740/1‑150 (MQ=255)
                                                                              ggCCCGGGAGCACGCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGt                                                                         >  2:480699/1‑150 (MQ=255)
                                                                                 ccGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGaa                                                                      >  2:648003/1‑150 (MQ=255)
                                                                                 ccGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCa                                                                                                                                            <  2:775658/80‑1 (MQ=255)
                                                                                  cGGCAGCACTCGAAGTGGTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAg                                                                     <  2:691372/150‑1 (MQ=255)
                                                                                  cGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAg                                                                     <  2:1981652/150‑1 (MQ=255)
                                                                                  cGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAg                                                                     <  2:691392/150‑1 (MQ=255)
                                                                                   atctgcaCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAgg                                                                    <  2:1977640/146‑1 (MQ=255)
                                                                                    gcagcaCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAANGGCTCAGCANCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGc                                                                   <  1:1048951/150‑1 (MQ=255)
                                                                                      agcaCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGa                                                                 >  1:963442/1‑150 (MQ=255)
                                                                                       gcaCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGAc                                                                                                                                              >  1:1756483/1‑72 (MQ=255)
                                                                                        caCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCg                                                                                                                                 <  2:556948/84‑1 (MQ=255)
                                                                                        caCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCg                                                                                                                                 >  1:556948/1‑84 (MQ=255)
                                                                                        caCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGt                                                                                      <  2:107448/127‑1 (MQ=255)
                                                                                        caCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAAc                                                               <  1:648003/150‑1 (MQ=255)
                                                                                         aCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAAcc                                                              <  1:915098/150‑1 (MQ=255)
                                                                                         aCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAAcc                                                              >  1:808656/1‑150 (MQ=255)
                                                                                          cTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGc                                                                                                                       >  1:2267776/1‑92 (MQ=255)
                                                                                          cTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGc                                                                                                                       <  2:2267776/92‑1 (MQ=255)
                                                                                          cTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAAccg                                                             >  2:2209634/1‑150 (MQ=255)
                                                                                           tCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAAccgc                                                            >  2:1856901/1‑150 (MQ=255)
                                                                                           tCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAAccgc                                                            >  1:1528296/1‑150 (MQ=255)
                                                                                             gAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGGCCGTGAAGGCGAAccgccg                                                          >  2:1587862/1‑150 (MQ=255)
                                                                                                gTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAg                                                       >  2:291630/1‑150 (MQ=255)
                                                                                                    tgctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAg                                                   >  2:2290829/1‑150 (MQ=255)
                                                                                                    tgctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAg                                                   >  1:1278062/1‑150 (MQ=255)
                                                                                                      ctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCANCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCACCCAggg                                                 >  1:1118094/1‑150 (MQ=255)
                                                                                                      ctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAggg                                                 <  1:2209634/150‑1 (MQ=255)
                                                                                                       tCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAg                                                                                                                              >  1:1653390/1‑72 (MQ=255)
                                                                                                       tCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAg                                                                                                                              <  2:1653390/72‑1 (MQ=255)
                                                                                                       tCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTca                                                                                                            >  1:1344110/1‑90 (MQ=255)
                                                                                                       tCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAggg                                                 >  1:388191/1‑149 (MQ=255)
                                                                                                       tCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAggg                                                 <  2:388191/149‑1 (MQ=255)
                                                                                                       tCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGACCAGCATGTGGGCCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCACCAGCGCAGCAGGCCTGGCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGc                                                >  1:1094581/1‑150 (MQ=255)
                                                                                                         aaGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGt                                                                         <  1:1454608/123‑1 (MQ=255)
                                                                                                          aGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGc                                                                                              <  1:2177316/101‑1 (MQ=255)
                                                                                                          aGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGc                                                                                              >  2:2177316/1‑101 (MQ=255)
                                                                                                          aGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAgg                                             <  1:1856901/150‑1 (MQ=255)
                                                                                                             gTTTGCCGGCAGGGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTAc                                          >  1:562884/1‑150 (MQ=255)
                                                                                                              tttGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACt                                         >  1:2056833/1‑150 (MQ=255)
                                                                                                                     gCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCa                                                    >  1:2022863/1‑132 (MQ=255)
                                                                                                                     gCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCa                                                    <  2:2022863/132‑1 (MQ=255)
                                                                                                                              ccGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTaccaac                         <  2:1528296/150‑1 (MQ=255)
                                                                                                                                 gcggcCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAaccac                      <  2:229491/150‑1 (MQ=255)
                                                                                                                                 gcggcCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAaccac                      <  2:1894561/150‑1 (MQ=255)
                                                                                                                                      cATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACAt                 >  1:1154761/1‑150 (MQ=255)
                                                                                                                                          gTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGc             <  2:1154761/150‑1 (MQ=255)
                                                                                                                                          gTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGc             >  2:2269801/1‑150 (MQ=255)
                                                                                                                                           tGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGcc            <  1:480699/150‑1 (MQ=255)
                                                                                                                                           accTGGTCCACATGCTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCg                                 <  2:1443912/127‑1 (MQ=255)
                                                                                                                                           accTGGTCCACATGCTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCg                                 >  1:1443912/3‑129 (MQ=255)
                                                                                                                                                gTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTa       >  2:20362/1‑150 (MQ=255)
                                                                                                                                                  ccAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACa     <  1:1386629/150‑1 (MQ=255)
                                                                                                                                                  ccAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACa     <  2:808656/150‑1 (MQ=255)
                                                                                                                                                   cAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTgcgc                                                                                  <  1:1049630/72‑1 (MQ=255)
                                                                                                                                                     gCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAGaa  >  1:864098/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCCGCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAGAA  >  NC_002947/1698703‑1699001

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: