Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 681,295 G→A 5.5% intergenic (+32/+37) PP_0583 → / ← PP_0584 AraC family transcriptional regulator/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947681,2950GA 5.5% 521.3 / 12.1 164intergenic (+32/+37)PP_0583/PP_0584AraC family transcriptional regulator/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base G (82/73);  new base A (7/2);  total (89/75)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 174 total aligned reads displayed.

GGAGGCGGGGGAGAATGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCC  >  NC_002947/681146‑681442
                                                                                                                                                     |                                                                                                                                                     
ggAGGCGGGGGAGAATGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCggg                                                                                                                                                       <  1:3117498/150‑1 (MQ=255)
ggAGGCGGGGGAGAATGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCggg                                                                                                                                                       <  1:2003734/150‑1 (MQ=255)
ggAGGCGGGGGAGAATGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCggg                                                                                                                                                       <  1:5882058/150‑1 (MQ=255)
 gAGGCGGGGGAGAATGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCgggg                                                                                                                                                      <  2:4251753/150‑1 (MQ=255)
 gAGGCGGGGGAGAATGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCgggg                                                                                                                                                      <  2:8035109/150‑1 (MQ=255)
  aGGCGGGGGAGAATGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGt                                                                                                                                                     >  2:6828377/1‑150 (MQ=255)
       ggggAGAATGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTAtac                                                                                                                                                >  1:5087118/1‑148 (MQ=255)
       ggggAGAATGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATcc                                                                                                                                                >  1:2780110/1‑150 (MQ=255)
       ggggAGAATGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATcc                                                                                                                                                <  2:6716969/150‑1 (MQ=255)
          gagaATGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGGGCGTTGGGGGGGttttggggg                                                                                                                                             >  1:1528444/1‑143 (MQ=255)
             aaTGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCTCCTCGGCCTACATTGCGGCATTCTAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGaaa                                                                                                                                          <  2:4039845/150‑1 (MQ=255)
             aaTGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGaaa                                                                                                                                          <  2:2379694/150‑1 (MQ=255)
             aaTGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGaaa                                                                                                                                          <  1:8113280/150‑1 (MQ=255)
             aaTGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGaaa                                                                                                                                          <  1:132387/150‑1 (MQ=255)
              aTGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAg                                                                                                                                         <  2:6963542/150‑1 (MQ=255)
              aTGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAg                                                                                                                                         <  2:6963572/150‑1 (MQ=255)
                gtgACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAggg                                                                                                                                       <  1:2991946/150‑1 (MQ=255)
                gtgACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAggg                                                                                                                                       <  1:4473342/150‑1 (MQ=255)
                 tgACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGc                                                                                                                                      <  2:4709666/150‑1 (MQ=255)
                 tgACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGc                                                                                                                                      <  1:1181396/150‑1 (MQ=255)
                 tgACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGATAACCCCGCGAAAGGGc                                                                                                                                      <  1:34568/150‑1 (MQ=39)
                   aCCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGGCCGTTCGCGGGGGTATCCGCGAAAGGGCCg                                                                                                                                    >  1:8052564/1‑150 (MQ=255)
                   aCCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCAGCGAAAGGGCCg                                                                                                                                    >  2:5748742/1‑150 (MQ=255)
                     cGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGc                                                                                                                                 >  2:5657096/1‑150 (MQ=255)
                      gAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGca                                                                                                                                >  2:5275047/1‑150 (MQ=255)
                      gAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGca                                                                                                                                >  2:1966709/1‑150 (MQ=255)
                        ggcggcGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATGGCGGCATTCAAGCAGCTATTGGGCGCGACACGGGGGGAATTGACGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGGGGGGTTATCCGCGAAAGGGCCGGcaca                                                                                                                              <  2:8260923/150‑1 (MQ=255)
                         gcggcgCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATACGCGAAAGGGCCGGCACAg                                                                                                                             >  2:3280246/1‑150 (MQ=255)
                          cggcgCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCATTAGACGGGTTATACGCGAAAGGGCCGGCACAgg                                                                                                                            >  1:4530367/1‑150 (MQ=255)
                           ggcgCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGc                                                                                                                           >  2:3672166/1‑150 (MQ=255)
                                       ggggTATGACTCCACCTCGGCCTACTTTGCGGAATTCAAGCAGCTATTGGGCGCGACCCGGGGGGAATTGAAGCTGTTTTTTTGGGTGTCCTGTGCGGGCCCTTGCGGGGGGTTAGCCGCGAACGGGCCGGCACAGGCCATGCAtcatca                                                                                                               <  2:3204744/149‑1 (MQ=255)
                                         ggTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCCTGCATCATCAgg                                                                                                             >  2:3059633/1‑150 (MQ=255)
                                         ggTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATACGCGAAAGGGCCGGCACAGGCCCTGCATCATCAgg                                                                                                             >  2:1846803/1‑150 (MQ=255)
                                          gTATGACTCCACCTCGGCCTACCTTGCGGCAGTCACGCAGCTTTTTGGCGCGCCCCGGGGGGCATTGAAGCTCTGAGTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGt                                                                                                            <  2:9071018/150‑1 (MQ=255)
                                              gACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCCGGGTTTTACGCGAAAGGGCCCGCACAGGCCATGCAACATCAGGTGccg                                                                                                        >  1:746437/1‑148 (MQ=255)
                                                cTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCATTCGCGGGGTTATCCGCGAAAGGGCCGGCACACGCCATGCATCATAAGGTGCGGaa                                                                                                      >  2:2649170/1‑150 (MQ=255)
                                                 tCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTGGCGGGGTTGTCCGCGAAAGGGCCGGCACAGGCCCAGCAGCATAAGGTGCGGAAc                                                                                                     >  2:30404/1‑150 (MQ=255)
                                                 cccACCTCGGCCTACCTTGCGGCTTTCACGCAGCTTTTTGGCGCGACACGGGGGGAATTGATGCTCTGATTTTGGGTTTCCTGTGCCGGCCCTTTCGGGGGGTTATCCGCGAACGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAAc                                                                                                     <  2:2673441/149‑1 (MQ=255)
                                                   cACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCATTGGCGGGGTTTTACGCGAAAAGGCCGGCCAACGCCAAGCATAAACACGTGCGGAAccc                                                                                                   >  1:2295711/1‑149 (MQ=255)
                                                      cTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGTACCTgcg                                                                                                >  2:3449434/1‑150 (MQ=255)
                                                      cTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCGTTAGCGGGGGTTTGCGCGCAAGGGCCGGCACAGGCCCAGCAACATCAGGGGCGGAACCAgcg                                                                                                >  1:1730392/1‑150 (MQ=255)
                                                                    cGGCATTCACGCCGCTATTGGGCGCGCCACGGGGGTAATGTACGCTCTTATTTGGGGTTACCTGTGCCGCCCCTGCCGCGGGTTTACCCGCTAACGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCa                                                                                  <  1:8044615/150‑1 (MQ=255)
                                                                     ggCATTCACGCAGCTATTGGGCGCGACACGGGGGTAATTGAAGCTCTGATTTTGGGTTTCCTGGGCGGGCCCTTTCGCGGGGTTATCCGCGAACGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAg                                                                                 <  1:4379890/150‑1 (MQ=255)
                                                                     ggCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGGCCTTTCGGGGGGTTATGCGCGTTAGGGGCGGGACAGGGCATGGATGATCATGTGCGGAACGTGCGGACCAGCGTATGCAg                                                                                 >  2:5460112/1‑150 (MQ=255)
                                                                     ggCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAg                                                                                 <  2:7511904/150‑1 (MQ=255)
                                                                     ggCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCATTCGCGGGGTTATCCGCGAAATGGCCGGCACCGGCCAAGCAACATCAGGTGCGGAACCTGCGCCCCAGCGTCTCCAg                                                                                 >  1:5418756/1‑150 (MQ=255)
                                                                      gCATTCAAGCAGCTTTTGGGCGCGACCCGGGGGGAATTGACGCTCTGATTTCGGGTTTCCTGTGCCGGCCCTTGCGGGGGGTTATCCGCGAACGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGc                                                                                <  1:1426287/150‑1 (MQ=255)
                                                                         ttCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGGTATACGCGAAATGGCCGGCACAGGCCATGCATCAACAGGGGCGGCAACAGCCCACCCGCGTATACAGcaca                                                                            >  1:4723465/1‑150 (MQ=255)
                                                                         ttCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCATTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGGGCGGAACCTGCGCACCAGCGTATCCAGCTCa                                                                            >  1:6073019/1‑150 (MQ=255)
                                                                         ttCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCATTCGCGGGGTTATCCGCGAAAGGGCCGGCACACGCCCGGCATCCACCCGTGCGGCACCGGCCCCCCACCGGATCCCGCCCa                                                                            >  1:7339370/1‑150 (MQ=255)
                                                                            aaGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATACAGCTCATTg                                                                         >  2:5896724/1‑150 (MQ=255)
                                                                               cagcTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACACGCCATGCATCATAAGGTGCGGAACCAGCCCACCACCGTAACCAGCTCATTGGAc                                                                      >  1:7453196/1‑150 (MQ=255)
                                                                                agcTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATACGCGAAAGGGCCGGCACAGGCCATGCATCAACAGGTGCGGAACCTGCGCACCAGCGTATACAGCTCAATTGACa                                                                     >  1:3759822/1‑150 (MQ=255)
                                                                                  ctttttGGCGCGACCCGGGGGTATTTGATGGTCTGTTTTTGGGTTTCGTGGGCGGGCCGTTTCGGGGGGTTATCCGCTAACGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGc                                                                   <  1:2041260/147‑1 (MQ=255)
                                                                                  cTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCGTTAGCCGGGGTATACGCGAAACGGCCGGCACAGGGCCAGCCTCAAAAGGGGCGGAACCAGCGCCCCAGCGGATCCACCGCATTTGACAGc                                                                   >  1:1897290/1‑150 (MQ=255)
                                                                                   tATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCAGCGCACCAGCGTATCCAGCTCATTGGACAGcc                                                                  >  2:8042643/1‑150 (MQ=255)
                                                                                   tATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCATTGGCGGGGGTGTACGCGAAAGGGCCCGCACACGCCCATCATCAAAAGGTGCGGAACCAGCGCACCAGCGGATCCAGCGCCATTGGaagcc                                                                   >  1:213484/1‑146 (MQ=255)
                                                                                     tttGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTAGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCAAGCATCATAAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTTGACAGCCCg                                                                >  2:8867632/1‑150 (MQ=255)
                                                                                       tGGCGCGCCCCCGGGGGATTTTAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGGGGGTTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGc                                                              <  1:6787553/150‑1 (MQ=255)
                                                                                       tGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTGCGGGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGc                                                              <  2:6381949/150‑1 (MQ=255)
                                                                                        ggCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGcc                                                             <  2:5226871/150‑1 (MQ=255)
                                                                                        ggCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGcc                                                             >  2:8150274/1‑150 (MQ=255)
                                                                                        ggCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCGTTGGCCGGGGGATGCGCGAAAGGGCCGGCCCCCGGCCATCCTCAACAGGTGCGGAACCGGCGCCCCACCGTATACAGCGCATTTGAAACCCCGGcc                                                             >  1:5801432/1‑150 (MQ=255)
                                                                                         gcgcGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCa                                                            <  2:7118713/150‑1 (MQ=255)
                                                                                           gcgACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGGACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAgg                                                          >  2:6695392/1‑150 (MQ=255)
                                                                                           gcgACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTAGCGGGGTTATACGCGAAAGGGCCCGCACAGGCCATGCATCATAAGGTGCGGAACCTGCGCACCAGCGGATCCAGCTCATTGGAAAGCCCCGCCacg                                                          >  1:719071/1‑148 (MQ=255)
                                                                                                acCGGGGTAATTGAAGCTCTGTTTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCtctc                                                     <  1:5298233/150‑1 (MQ=255)
                                                                                                   gggggAATTTAAGCTCTGTTTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGa                                                  <  2:6339462/150‑1 (MQ=255)
                                                                                                   gggggAATTGAAGCTCTGTTTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTCTCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGa                                                  <  2:2941220/150‑1 (MQ=255)
                                                                                                   gggggAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGa                                                  <  1:7567173/150‑1 (MQ=255)
                                                                                                   gggggAATTGAAGATCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCTAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGa                                                  <  1:7567566/150‑1 (MQ=255)
                                                                                                    gggtAATTGACGCGCTTTTTTTGGGTTGCGTGTGCGGGCCGTTCGGGGGGGTTATCCGCTAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAg                                                 <  2:1544362/150‑1 (MQ=255)
                                                                                                       tatttGAAGCTCTTTTTTTGGTTTGCGTGGGCGGGCCCTTTGGGGGGTTTCCCCGCTAACGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTcc                                              <  1:8563505/147‑1 (MQ=255)
                                                                                                         aTTGACGCTCTTTTTTTGGTTTTCCTGTGCGGGCCCTTTCGGGGGGTTATCCGCGAACGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTg                                            <  1:3613169/150‑1 (MQ=255)
                                                                                                         aTTGAAGCTCTGTTTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTg                                            <  2:7200299/150‑1 (MQ=255)
                                                                                                         aTTGAAGCTCTGATTTTCGGGTTCCTGTGCCGGCCCTTTCGCGGGGTTACCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTg                                            <  1:6084975/150‑1 (MQ=255)
                                                                                                               gCTCTGATTTTCGGTTTCCTGTGCCGGCCCATTCCCGGGGTTATCCGCGAAAGGGCCCGCACAGGCCATGCATCCTCACGTGCGGAACCAGCGCACCAGCGTATCCAGCACAATGGGCCGCCCCGCCACGCTCAAGGAGTCCAGgcgcgc                                      >  1:3727052/1‑150 (MQ=255)
                                                                                                                 tctGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGGTATCCGCGAAAGGGCCGGCCCACGCCATGCATCATCAGGTGCGGCACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGcgccg                                    >  2:4880473/1‑150 (MQ=255)
                                                                                                                     aTTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCAGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTg                                >  2:624335/1‑150 (MQ=255)
                                                                                                                      ttttCGGTTTCCTGTGCCGGCCCTTTCGGGGGGTTATCCGCTAACGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTgg                               <  1:6615601/150‑1 (MQ=255)
                                                                                                                        ttCGGTTGCCTGTGCCGGCCCTTTCGGGGGGTTAGCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGc                             <  2:3090359/150‑1 (MQ=255)
                                                                                                                          cGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCCTCATCCGGTGCGGAACCTGCGCACCAGCGGATCCAGCTCCTTGGACAGCCCGGCCCGGCTCTCGGCGTCCGGGCGCGCCGCCGGGGCCa                           >  2:8616092/1‑150 (MQ=255)
                                                                                                                                 cTGTGCCGGCCATTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGc                    <  1:8154371/150‑1 (MQ=255)
                                                                                                                                  tgtgCCGGCCCTTTAGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGcc                   >  2:737894/1‑150 (MQ=255)
                                                                                                                                    tgCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCACTCGGAGTCCTTGCGCGCCGCCTGGGCCAGCTCGGccac                 >  2:5974108/1‑150 (MQ=255)
                                                                                                                                    tgCCGGCCCTTTCGCGGATAACCCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTc                                                                             <  1:2950718/90‑1 (MQ=21)
                                                                                                                                     gCCGGCCCTTTCGCGGATAACCCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGccacc                >  2:5269825/1‑150 (MQ=38)
                                                                                                                                      ccGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGccacca               >  2:481580/1‑150 (MQ=255)
                                                                                                                                      ccGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGccacca               >  2:110122/1‑150 (MQ=255)
                                                                                                                                       cGGCCCTTTCGCGGATAACCCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGACACCaa              >  2:147480/1‑150 (MQ=38)
                                                                                                                                       cGGCCCTTTCGCGGATAACCCCGCGAAAGGGCCGGCACAGGCCAAGCAACAACAGGTGCGGAACCAGCGCACCACCGTATCCAACTCATTGGACACCCCCGCCAGGCACACCGAGTCCAGGCCCGCCCCCAGGGCCAGCACCGCCACCaa              >  1:7154620/1‑150 (MQ=16)
                                                                                                                                           ccTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTgc          >  1:6893247/1‑150 (MQ=255)
                                                                                                                                            cTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTgcg         >  2:7300221/1‑150 (MQ=255)
                                                                                                                                            cTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTgcg         >  2:8515989/1‑150 (MQ=255)
                                                                                                                                            cTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTgcg         >  1:3992067/1‑150 (MQ=255)
                                                                                                                                             tttCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTgcgc        >  2:4744473/1‑150 (MQ=255)
                                                                                                                                             tttCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTgcgc        >  1:1178938/1‑150 (MQ=255)
                                                                                                                                              ttCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGTCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCa       >  2:7978892/1‑150 (MQ=255)
                                                                                                                                                  cgGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATccc   >  2:5639694/1‑150 (MQ=255)
                                                                                                                                                   gggtTTATCCGCTAACGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATcccc  <  1:2600543/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGAGGCGGGGGAGAATGTGACCGAGGCGGCGCTGGGGTGTGGGTATGACTCCACCTCGGCCTACATTGCGGCATTCAAGCAGCTATTTGGCGCGACACCGGGGGAATTGAAGCTCTGATTTTCGGTTTCCTGTGCCGGCCCTTTCGCGGGGTTATCCGCGAAAGGGCCGGCACAGGCCATGCATCATCAGGTGCGGAACCTGCGCACCAGCGTATCCAGCTCATTGGACAGCCCGGCCAGGCTCTCGGAGTCCTGGCGCGCCGCCTGGGCCAGCTCGGCCACCAACTGCGCATCCCC  >  NC_002947/681146‑681442

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