Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,698,856:1 +C 6.3% coding (562/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,8561.C 6.3% 460.4 / 23.6 111coding (562/1014 nt)cheBCchemotaxis response regulator protein‑glutamate methylesterase
Reads supporting (aligned to +/- strand):  ref base . (51/53);  new base C (4/3);  total (55/56)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

Only 100 of 114 total aligned reads displayed.

GGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAGAATGGC  >  NC_002947/1698710‑1699005
                                                                                                                                                      |                                                                                                                                                     
 gTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGa                                                                                                                                                   >  2:2868825/1‑150 (MQ=255)
  tCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGAc                                                                                                                                                  >  2:1377968/1‑150 (MQ=255)
    ggTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGAATGTGGACCa                                                                                                                                                <  2:2423016/150‑1 (MQ=255)
      tGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAgg                                                                                                                                              <  1:2013729/150‑1 (MQ=255)
      tGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAgg                                                                                                                                              <  2:2573280/150‑1 (MQ=255)
           gCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGACATCGTGCTGGTCCAGCATGTGGACCAGgtgtct                                                                                                                                         <  2:897038/150‑3 (MQ=255)
             agcGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGTTATTGGGGCTTCGGTGGGTG‑‑GCCCGGCGGCGCTCGAAGTGCTGCTGAAGTGTTTGCCGGCAGCGGTGCCGGCGGCCAGCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGc                                                                                                                                       <  2:346465/148‑1 (MQ=255)
             aCCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGc                                                                                                                                       >  1:1110660/1‑150 (MQ=255)
             aCCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGc                                                                                                                                       <  2:678157/150‑1 (MQ=255)
               cGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCAGTTGAAGTGCTGCTCCAGGGTTTGCCGGGAGCGTTTTCGGCGGCAATCGTGCTGGTCCAGCATGTGGACCAGGAGTTTGCCg                                                                                                                                     >  1:1449472/1‑150 (MQ=255)
                gCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGc                                                                                                                                    >  1:2393932/1‑150 (MQ=255)
                 cTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCATTCGAAGTGCTGCTCAAGGGTTTGCCGTCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCa                                                                                                                                   <  1:2675003/150‑1 (MQ=255)
                  tGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAg                                                                                                                                  <  1:287924/150‑1 (MQ=255)
                     gcgAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCa                                                                                                                               >  2:119559/1‑150 (MQ=255)
                     gcgAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCa                                                                                                                               >  1:1129687/1‑150 (MQ=255)
                             tCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTGGAGGCCGACAG‑‑GTCGGAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGaa                                                                                                                       <  2:2384698/150‑1 (MQ=255)
                               ccAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATg                                                                                                                     >  2:2268458/1‑150 (MQ=255)
                                cAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGgt                                                                                                                                             <  1:467298/125‑1 (MQ=255)
                                cAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGgt                                                                                                                                             >  2:467298/1‑125 (MQ=255)
                                    gccgcGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTca                                                                                                                <  1:2035417/150‑1 (MQ=255)
                                       gcgGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGgtt                                                                                                                                            >  1:2744141/1‑118 (MQ=255)
                                         gggTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGGACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCagcagc                                                                                                           >  2:352639/1‑150 (MQ=255)
                                         gggTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCagcagc                                                                                                           >  1:1043412/1‑150 (MQ=255)
                                               gtggCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCgccgca                                                                                                     >  1:2792316/1‑150 (MQ=255)
                                               gtggCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCgcagca                                                                                                     >  1:2600600/1‑150 (MQ=255)
                                                tggCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTgg                                                                                                                                                    <  1:2876695/102‑1 (MQ=255)
                                                tggCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTgg                                                                                                                                                    >  2:2876695/1‑102 (MQ=255)
                                                    aaTCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATgg                                                                                                                    <  1:2695795/130‑1 (MQ=255)
                                                    aaTCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATgg                                                                                                                    >  2:2695795/1‑130 (MQ=255)
                                                    aaTCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCt                                                                                                >  2:2443855/1‑150 (MQ=255)
                                                                ggcggGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCTAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCt                                                                                    >  1:2021045/1‑150 (MQ=255)
                                                                          gCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGCCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTg                                                                                   >  2:2052055/1‑143 (MQ=255)
                                                                            ccGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCCGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGaa                                                                          >  2:687087/1‑150 (MQ=255)
                                                                            ccGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGaa                                                                          <  2:2600600/150‑1 (MQ=255)
                                                                             gagcagcaCTTCGAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAg                                                                         >  1:214859/3‑150 (MQ=255)
                                                                             cGGCAGCACTCGAAGTGCTGCTCAAGGGTTTTCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGAATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAg                                                                         <  1:1370002/150‑1 (MQ=255)
                                                                             cGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAg                                                                         >  2:2112808/1‑150 (MQ=255)
                                                                             cGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAg                                                                         <  2:2065793/150‑1 (MQ=255)
                                                                                 agcaCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGa                                                                     <  1:449400/150‑1 (MQ=255)
                                                                                  gcaCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGaa                                                                    >  2:617508/1‑150 (MQ=255)
                                                                                       tCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAAccgc                                                                <  2:2694674/150‑1 (MQ=255)
                                                                                        cGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAAccgcc                                                               <  2:2310343/150‑1 (MQ=255)
                                                                                           aGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTTGCCTGTGAAGGCGAACCGCCGCa                                                            <  1:621071/150‑1 (MQ=255)
                                                                                                tgctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAgg                                                                                                   <  1:1588144/106‑1 (MQ=255)
                                                                                                tgctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAgg                                                                                                   >  2:1588144/1‑106 (MQ=255)
                                                                                                tgctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCATGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAg                                                       >  1:1879264/1‑150 (MQ=255)
                                                                                                tgctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCa                                                                                                                                                <  1:2550883/61‑1 (MQ=37)
                                                                                                  ctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAggg                                                     >  1:97667/1‑150 (MQ=255)
                                                                                                   tCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTTGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGc                                                    >  1:2917708/1‑150 (MQ=255)
                                                                                                   tCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATgg                                                                                                                    >  1:1043592/1‑86 (MQ=255)
                                                                                                   tCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATgg                                                                                                                    <  2:1043592/86‑1 (MQ=255)
                                                                                                   tCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCa                                                                                                                                                <  2:2015028/58‑1 (MQ=37)
                                                                                                   tCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCa                                                                                                                                                >  1:2015028/1‑58 (MQ=37)
                                                                                                    cAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCTAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCa                                                   <  2:2021045/150‑1 (MQ=255)
                                                                                                     aaGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGAc                                                                                                                                                  <  2:1150119/54‑1 (MQ=255)
                                                                                                     aaGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGAc                                                                                                                                                  >  1:1150119/1‑54 (MQ=255)
                                                                                                        ggTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTa                                               <  1:1432813/150‑1 (MQ=255)
                                                                                                             gCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCt                                          <  1:2112808/150‑1 (MQ=255)
                                                                                                              ccGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAAt                                                                                                                      >  1:3003817/1‑73 (MQ=255)
                                                                                                              ccGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAAt                                                                                                                      <  2:3003817/73‑1 (MQ=255)
                                                                                                                ggCAGCGTTTCCGGCGGGCATGGGGCCGGTCCACCATGTGGACCAGGTGTGTCCCGCAGGAATGGCCGAATGGCTCAGCAGCGCAGCCGGCTTGCCGGTGCGCTTGGCCCGTGAAGGCGACCCGCCGCAGCCAGGGCAGGTACTCCTGGc                                       <  2:439861/150‑1 (MQ=255)
                                                                                                                ggCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGc                                       <  1:3058980/150‑1 (MQ=255)
                                                                                                                 gCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCGGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAg                                                  >  1:714868/1‑138 (MQ=255)
                                                                                                                 gCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAg                                                  <  2:714868/138‑1 (MQ=255)
                                                                                                                  cAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCg                                     <  2:1129687/150‑1 (MQ=255)
                                                                                                                       tttCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGaa                                                                          <  1:1065219/108‑1 (MQ=255)
                                                                                                                       tttCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGaa                                                                          >  2:1065219/1‑108 (MQ=255)
                                                                                                                       tttCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTacc                                <  2:2708700/150‑1 (MQ=255)
                                                                                                                         tCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTaccaa                              <  1:119559/150‑1 (MQ=255)
                                                                                                                            ggcggcCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTaccaacca                           >  2:1671981/1‑150 (MQ=255)
                                                                                                                               ggcCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGt                                   >  2:911111/1‑139 (MQ=255)
                                                                                                                               ggcCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGt                                   <  1:911111/139‑1 (MQ=255)
                                                                                                                               ggcCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAAccacca                        >  1:152447/1‑150 (MQ=255)
                                                                                                                                gcCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAg                                                  >  1:2143519/1‑123 (MQ=255)
                                                                                                                                gcCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAg                                                  <  2:2143519/123‑1 (MQ=255)
                                                                                                                                 ccATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACcaca                      <  1:1970705/150‑1 (MQ=255)
                                                                                                                                  cATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACAt                     <  2:214859/150‑1 (MQ=255)
                                                                                                                                   aTCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACAtt                    >  2:144734/1‑150 (MQ=255)
                                                                                                                                     cGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCg                  >  1:2416205/1‑150 (MQ=255)
                                                                                                                                       tGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGcc                >  1:2448714/1‑150 (MQ=255)
                                                                                                                                       accTGGTCCACATGCTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGAGCAGGTACTCCTgg                                        >  2:2641632/3‑126 (MQ=255)
                                                                                                                                       accTGGTCCACATGCTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGAGCAGGTACTCCTgg                                        <  1:2641632/124‑1 (MQ=255)
                                                                                                                                        gCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCt               <  1:1397218/150‑1 (MQ=255)
                                                                                                                                        ccTGGTCCAGCATGTGGACCAGGTGTTTGCGGCAGGCATGGCCGAATGGCTCAGCAGCGCGGCAGGCCTGCCAGTGCGCCTTGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCATGTACTCCTGGCCTGTACCCACCACCACATTCGCCt               <  2:1449472/149‑1 (MQ=255)
                                                                                                                                         cTGGTCCACATGCTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTg              >  1:1122840/1‑150 (MQ=255)
                                                                                                                                          tGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGc             >  1:1962682/1‑150 (MQ=255)
                                                                                                                                            gTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGc                                                                       >  2:128104/1‑90 (MQ=255)
                                                                                                                                            gTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAgg                                                 >  2:2220335/1‑112 (MQ=255)
                                                                                                                                            gTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAgg                                                 <  1:2220335/112‑1 (MQ=255)
                                                                                                                                              ccAGCATGTGGACCAGGTGTTTGCAGCAAGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACa         >  1:2094766/1‑150 (MQ=255)
                                                                                                                                               cAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAg        >  2:2743781/1‑150 (MQ=255)
                                                                                                                                               cAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAg        >  1:612242/1‑150 (MQ=255)
                                                                                                                                               cAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCAGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTa                                               >  1:2496188/1‑111 (MQ=255)
                                                                                                                                               cAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCAGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTa                                               <  2:2496188/111‑1 (MQ=255)
                                                                                                                                                aGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGc                                       <  1:75156/118‑1 (MQ=255)
                                                                                                                                                aGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGc                                       >  2:75156/1‑118 (MQ=255)
                                                                                                                                                   aTGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCGGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCAGCAGCCATGGCAGGTGCTCCGGGCCGGTACCAACCACCACATTCGCCTGCTACAGAATg    <  1:60918/150‑1 (MQ=255)
                                                                                                                                                    tgtgGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAGAATgg   <  2:1962682/150‑1 (MQ=255)
                                                                                                                                                     gtgGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCTCCACATTCGCCTGCTACAGAATGGc  >  2:1086262/1‑150 (MQ=255)
                                                                                                                                                     gtgGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAGAATGGc  <  1:1377968/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAGAATGGC  >  NC_002947/1698710‑1699005

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: