Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,698,852 Δ1 bp 6.3% 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. 6.3% 351.7 / 17.4 111coding (558/1014 nt)cheBCchemotaxis response regulator protein‑glutamate methylesterase
Reads supporting (aligned to +/- strand):  ref base G (52/51);  new base . (4/3);  total (56/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 111 total aligned reads displayed.

GCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAG  >  NC_002947/1698706‑1698999
                                                                                                                                                     |                                                                                                                                                    
gCAAGGTCGGTGGCTGCACCGCTGCGCGAGGATTGCGATTTGCTGGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCAt                                                                                                                                                   >  1:2122653/1‑150 (MQ=255)
     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                                                                                                                                        <  2:2573280/150‑1 (MQ=255)
          tGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAgg                                                                                                                                        <  1:2013729/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)
                     cTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCATTCGAAGTGCTGCTCAAGGGTTTGCCGTCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCa                                                                                                                             <  1:2675003/150‑1 (MQ=255)
                      tGCGTGAGGCTTCCCAGCGCCGCGTGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAg                                                                                                                            <  2:2181898/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)
                                  cccAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAAt                                                                                                                <  2:1242242/150‑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)
                                           gcgGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGgtg                                                                                                                                      <  2:2744141/119‑1 (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                                                                                                                                              >  2:2876695/1‑102 (MQ=255)
                                                    tggCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTgg                                                                                                                                              <  1:2876695/102‑1 (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)
                                                            gggTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCg                                                                                      <  1:2443855/150‑1 (MQ=255)
                                                                    ggcggGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCTAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCt                                                                              >  1:2021045/1‑150 (MQ=255)
                                                                     gcggGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTg                                                                             <  1:875605/150‑1 (MQ=255)
                                                                              gCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTg                                                                             <  1:2052055/143‑1 (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:2065793/150‑1 (MQ=255)
                                                                                 cGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAg                                                                   >  2:2112808/1‑150 (MQ=255)
                                                                                     agcaCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGa                                                               <  1:449400/150‑1 (MQ=255)
                                                                                      gcaCTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGaa                                                              >  2:617508/1‑150 (MQ=255)
                                                                                            cGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAAccgcc                                                         <  2:2310343/150‑1 (MQ=255)
                                                                                               aGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTTGCCTGTGAAGGCGAACCGCCGCa                                                      <  1:621071/150‑1 (MQ=255)
                                                                                                    tgctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAgg                                                                                             >  2:1588144/1‑106 (MQ=255)
                                                                                                    tgctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAgg                                                                                             <  1:1588144/106‑1 (MQ=255)
                                                                                                    tgctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCATGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAg                                                 >  1:1879264/1‑150 (MQ=255)
                                                                                                    tgctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCa                                                                                                                                          <  1:2550883/61‑1 (MQ=37)
                                                                                                    tgctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCACATGCTGGACCa                                                                                                                                          >  2:2550883/1‑61 (MQ=37)
                                                                                                      ctCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAggg                                               >  2:60918/1‑150 (MQ=255)
                                                                                                      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)
                                                                                                        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)
                                                                                                                    ggCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGc                                 <  1:3058980/150‑1 (MQ=255)
                                                                                                                     gCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCGGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAg                                            >  1:714868/1‑138 (MQ=255)
                                                                                                                     gCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAg                                            <  2:714868/138‑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)
                                                                                                                              ccGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTaccaac                       <  1:352639/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                                  <  1:2641632/124‑1 (MQ=255)
                                                                                                                                           accTGGTCCACATGCTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGAGCAGGTACTCCTgg                                  >  2:2641632/3‑126 (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)
                                                                                                                                                gTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGc                                                                 <  1:128104/90‑1 (MQ=255)
                                                                                                                                                gTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAgg                                           >  2:2220335/1‑112 (MQ=255)
                                                                                                                                                gTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAgg                                           <  1:2220335/112‑1 (MQ=255)
                                                                                                                                                  ccAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACa   <  1:1604136/150‑1 (MQ=255)
                                                                                                                                                  ccAGCATGTGGACCAGGTGTTTGCAGCAAGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACa   >  1:2094766/1‑150 (MQ=255)
                                                                                                                                                   cAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAg  >  1:612242/1‑150 (MQ=255)
                                                                                                                                                   cAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAg  >  2:2743781/1‑150 (MQ=255)
                                                                                                                                                   cAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCAGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTa                                         <  2:2496188/111‑1 (MQ=255)
                                                                                                                                                   cAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCAGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTa                                         >  1:2496188/1‑111 (MQ=255)
                                                                                                                                                    aGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGc                                 <  1:75156/118‑1 (MQ=255)
                                                                                                                                                    aGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGc                                 >  2:75156/1‑118 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAAGGTCGGTGGCTGCACCGCTGCGCGAGGCTTCCCAGCGCCGCGGGTTGGTGGCAATCGGGTCTTCGGCGGGTG‑‑GCCCGGCAGCACTCGAAGTGCTGCTCAAGGGTTTGCCGGCAGCGTTTCCGGCGGCCATCGTGCTGGTCCAGCATGTGGACCAGGTGTTTGCCGCAGGCATGGCCGAATGGCTCAGCAGCGCAGCAGGCCTGCCGGTGCGCCTGGCCCGTGAAGGCGAACCGCCGCAGCCAGGGCAGGTACTCCTGGCCGGTACCAACCACCACATTCGCCTGCTACAG  >  NC_002947/1698706‑1698999

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