Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,698,414 G→A 6.8% Q40Q (CAG→CAA 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,4140GA 6.8% 320.4 / 8.6 111Q40Q (CAG→CAAcheBCchemotaxis response regulator protein‑glutamate methylesterase
Reads supporting (aligned to +/- strand):  ref base G (49/53);  new base A (5/3);  total (54/57)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 113 total aligned reads displayed.

TGTTGTGGAACTGATCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAACGTCCACCGGG  >  NC_002947/1698265‑1698562
                                                                                                                                                     |                                                                                                                                                     
tgttgtGGAACTGATCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAg                                                                                                                                                       >  1:673492/1‑150 (MQ=255)
tgttgtGGAACTGATCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAg                                                                                                                                                       >  2:1197270/1‑150 (MQ=255)
  ttgtGGAACTGATCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCg                                                                                                                                                     <  1:911944/150‑1 (MQ=255)
    gtgGAACTGATCGGAGTTGCTCAGGGATGTAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGtt                                                                                                                                                   <  1:1819577/150‑1 (MQ=255)
    gtgGAACTGATCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGtt                                                                                                                                                   <  1:1197270/150‑1 (MQ=255)
     tgGAACTGATCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTg                                                                                                                                                  >  2:2439381/1‑150 (MQ=255)
       gAACTGATCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGGTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGGGCCGAGGCGGTGCAGCGTTGCa                                                                                                                                                >  2:277109/1‑150 (MQ=255)
           tGATCGGAGTTGCTCAGGGATGAAGCGGGNCATCGTCAATGACATGCCCNTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAANGCTGCACCgc                                                                                                                                            >  1:2275670/1‑149 (MQ=255)
           tGATCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAACGCTGCACCgc                                                                                                                                            <  1:1747665/150‑2 (MQ=255)
            gATCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAg                                                                                                                                           <  1:105708/150‑1 (MQ=255)
             aTCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGc                                                                                                                                          >  2:599152/1‑150 (MQ=255)
             aTCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGc                                                                                                                                          >  2:1393211/1‑150 (MQ=255)
                 gAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGtt                                                                                                                                      >  2:1403349/1‑150 (MQ=255)
                   gTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTAc                                                                                                                                    <  1:1256527/150‑1 (MQ=255)
                   gTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTAc                                                                                                                                    <  1:2320382/150‑1 (MQ=255)
                        tCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGAc                                                                                                                              >  1:48792/1‑150 (MQ=255)
                        tCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGAc                                                                                                                              <  1:249997/150‑1 (MQ=255)
                        tCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGAc                                                                                                                              >  1:601800/1‑150 (MQ=255)
                        tCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGCCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTAGCGGAc                                                                                                                              >  1:649714/1‑150 (MQ=255)
                          aGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCt                                                                                                                            <  1:99137/150‑1 (MQ=255)
                             gATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGAt                                                                                                                         <  2:1903903/150‑1 (MQ=255)
                              aTGAAAATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGagc                                                                                                                        <  2:1009168/150‑3 (MQ=255)
                                 aaGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTc                                                                                                                     <  1:1955372/150‑1 (MQ=255)
                                   gATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCAt                                                                                                                   <  1:736425/150‑1 (MQ=255)
                                   gATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCAt                                                                                                                   <  2:2449009/150‑1 (MQ=255)
                                        ccATCGTCAATGACATGCCCCCGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGAACGCACACCAGGTCGTGTGGGGGGCCAGCAACGGTGCCGATGCGGTGCAGCGTTGCACTGAGTAGTTACCGGACCTGATCCTCATGGAcc                                                                                                              <  1:716665/150‑1 (MQ=255)
                                           tCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCTGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCATCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGa                                                                                                           <  2:2190056/150‑1 (MQ=255)
                                           tCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGa                                                                                                           <  1:2439381/150‑1 (MQ=255)
                                             gTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATc                                                                                                         >  1:1952659/1‑150 (MQ=255)
                                                 aTGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGc                                                                                                     >  1:527141/1‑150 (MQ=255)
                                                 aTGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAACGCTGCACCgc                                                                                                                                            <  1:2212173/112‑2 (MQ=255)
                                                 aTGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAACGCTGCACCgc                                                                                                                                            >  2:2212173/1‑111 (MQ=255)
                                                    aCATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAg                                                                                                                                        >  1:718806/1‑113 (MQ=255)
                                                    aCATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAg                                                                                                                                        <  2:718806/113‑1 (MQ=255)
                                                    aCATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCAt                                                                                                                   <  1:1600870/133‑1 (MQ=255)
                                                    aCATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCgg                                                                                                  <  2:2447993/150‑1 (MQ=255)
                                                     cATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGt                                                                                                 >  2:2134850/1‑150 (MQ=255)
                                                      aTGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTg                                                                                                >  1:1492320/1‑150 (MQ=255)
                                                        gCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGAt                                                                                              <  2:2207733/150‑1 (MQ=255)
                                                         ccccTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATg                                                                                             <  2:673492/150‑1 (MQ=255)
                                                          cccTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTCCCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATgg                                                                                            >  1:2161381/1‑150 (MQ=255)
                                                           ccTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATg                                                                                                                  >  1:1810773/1‑127 (MQ=255)
                                                           ccTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATg                                                                                                                  <  2:1810773/127‑1 (MQ=255)
                                                           ccTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGa                                                                                           >  2:1540389/1‑150 (MQ=255)
                                                           ccTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTCCTCATGCCGATGaagca                                                                                           >  2:1093320/1‑146 (MQ=255)
                                                            cTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGAc                                                                                          >  1:1806776/1‑150 (MQ=255)
                                                             tGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACg                                                                                         >  2:2307203/1‑150 (MQ=255)
                                                               gCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTg                                                                                                            >  1:1019399/1‑129 (MQ=255)
                                                               gCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTg                                                                                                            <  2:1019399/129‑1 (MQ=255)
                                                                  gTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTc                                                                                    <  1:599152/150‑1 (MQ=255)
                                                                  gTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGGGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTc                                                                                    <  2:102377/150‑1 (MQ=255)
                                                                   tCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCg                                                                                   >  1:1427765/1‑150 (MQ=255)
                                                                     gAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAg                                                                                 <  1:2173708/150‑1 (MQ=255)
                                                                     gAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAg                                                                                 <  1:1093320/150‑1 (MQ=255)
                                                                     gAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAg                                                                                 <  2:601800/150‑1 (MQ=255)
                                                                     gAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAg                                                                                 <  2:675677/150‑1 (MQ=255)
                                                                           ttGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAAcc                                                                           >  1:1039020/1‑150 (MQ=255)
                                                                              cgcCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCcgg                                                                        <  2:1039020/150‑1 (MQ=255)
                                                                               gcCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGtt                                                                                                                                                   >  2:2370617/1‑75 (MQ=255)
                                                                               gcCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGtt                                                                                                                                                   <  1:2370617/75‑1 (MQ=255)
                                                                               gcCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCcggc                                                                       >  1:88270/1‑150 (MQ=255)
                                                                               gcCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCcggc                                                                       >  2:2003946/1‑150 (MQ=255)
                                                                               gcCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCcggc                                                                       >  2:1364045/1‑150 (MQ=255)
                                                                                      ggTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATgg                                                                                                                 <  2:1680692/101‑1 (MQ=255)
                                                                                      ggTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATgg                                                                                                                 >  1:1680692/1‑101 (MQ=255)
                                                                                        tCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGAc                                                                                                                              <  2:565143/86‑1 (MQ=255)
                                                                                        tCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCcggc                                                                       >  2:391009/1‑141 (MQ=255)
                                                                                        tCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCcggc                                                                       <  1:391009/141‑1 (MQ=255)
                                                                                        tCGCCCTGGAGCCCGCACACCAGGTCGGGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGAc                                                                                                                              >  1:565143/1‑86 (MQ=255)
                                                                                         cGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCAAGCAGCGACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGc                                                             <  1:290480/150‑1 (MQ=255)
                                                                                                  gCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAAcccc                                                    >  2:2419538/1‑150 (MQ=255)
                                                                                                  gCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGGGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAAcccc                                                    <  1:277109/150‑1 (MQ=255)
                                                                                                    ccGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGt                                                  >  1:1855172/1‑150 (MQ=255)
                                                                                                     cGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGgt                                                                                       <  1:231769/112‑1 (MQ=255)
                                                                                                     cGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGgt                                                                                       >  2:231769/1‑112 (MQ=255)
                                                                                                     cGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTCCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGgt                                                                                       >  2:1957484/1‑112 (MQ=255)
                                                                                                     cGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTCCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGgt                                                                                       <  1:1957484/112‑1 (MQ=255)
                                                                                                       cacaCCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATg                                                                                             <  1:1965634/104‑1 (MQ=255)
                                                                                                        acacCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTgcgc                                              >  1:943942/1‑150 (MQ=255)
                                                                                                                  gtgtGGGTGGCCAGCAACGGTGCCGCGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATc                                    <  1:2210515/150‑1 (MQ=255)
                                                                                                                  gtgtGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCgg                                                                                                  <  1:124041/88‑1 (MQ=255)
                                                                                                                  gtgtGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCgg                                                                                                  >  2:124041/1‑88 (MQ=255)
                                                                                                                  gtgtGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATc                                    <  2:2013723/150‑1 (MQ=255)
                                                                                                                      gggtggCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCGTCGTCa                                >  1:320326/1‑150 (MQ=255)
                                                                                                                       ggtggCCAGCAACGGTGCCGAGGCGGTGCAGAGTAGAAACAAGCAGTACCCGCATTGTGATTTCATGGACCTGGGCATGGCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCAc                               >  1:2125851/1‑150 (MQ=255)
                                                                                                                        gtggCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACg                              >  2:1837397/1‑150 (MQ=255)
                                                                                                                         tggCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACgg                             >  2:1934256/1‑150 (MQ=255)
                                                                                                                                         cGAGGCGGTGCAGCGTTGCTCCGCGCAGTTACCGGACCTGATCCTCATGGATCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGaa             <  2:1098756/150‑1 (MQ=255)
                                                                                                                                          gAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAAc            >  1:296459/1‑150 (MQ=255)
                                                                                                                                          gAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATCACGGTGATGGAAGTTGCTTAGGCGACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAAc            >  1:884965/1‑150 (MQ=255)
                                                                                                                                             tcGGTGCAACGCTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAACGTc         >  2:1763826/2‑150 (MQ=255)
                                                                                                                                             tcGGTGCAACGCTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCATAAACAGAACGTc         >  2:1996546/2‑150 (MQ=255)
                                                                                                                                             gCGGTGCAGCGTTGCACAGCGCACTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAACGTc         <  1:154091/150‑1 (MQ=255)
                                                                                                                                                gTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAACGTCCAc      >  2:1231327/1‑150 (MQ=255)
                                                                                                                                                   cAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACCGAACGTCCACCgg   <  1:2133589/150‑1 (MQ=255)
                                                                                                                                                   cAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAACGTCCACCgg   <  1:2003946/150‑1 (MQ=255)
                                                                                                                                                   cAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAACGTCCACCgg   <  1:1364045/150‑1 (MQ=255)
                                                                                                                                                    aGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGc                                                             <  1:2067143/91‑1 (MQ=255)
                                                                                                                                                    aGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGc                                                             >  2:2067143/1‑91 (MQ=255)
                                                                                                                                                    aGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAACGTCCACCggg  >  2:1289683/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTTGTGGAACTGATCGGAGTTGCTCAGGGATGAAGATCGCCATCGTCAATGACATGCCCCTGGCCGTCGAAGCGTTGCGCCGGGCGGTCGCCCTGGAGCCCGCACACCAGGTCGTGTGGGTGGCCAGCAACGGTGCCGAGGCGGTGCAGCGTTGCACCGAGCAGTTACCGGACCTGATCCTCATGGACCTGATCATGCCGGTGATGGACGGTGTCGAGGCAACCCGGCGGATCATGGCCGAAACCCCGTGCGCCATTGTCATCGTCACGGTCGACCGTAAACAGAACGTCCACCGGG  >  NC_002947/1698265‑1698562

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: