Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 732,464 G→T 11.8% R667R (CGG→CGT clpB → chaperone protein ClpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947732,4640GT11.8% 739.9 / 66.8 242R667R (CGG→CGTclpBchaperone protein ClpB
Reads supporting (aligned to +/- strand):  ref base G (111/102);  new base T (11/18);  total (122/120)
Fisher's exact test for biased strand distribution p-value = 1.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

TGGCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGAC  >  NC_002947/732316‑732613
                                                                                                                                                    |                                                                                                                                                      
tGGCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTgcggc                                                                                                                                                       <  1:3526019/150‑1 (MQ=255)
 ggCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTgcggcg                                                                                                                                                      >  1:2432918/1‑150 (MQ=255)
    cGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCaa                                                                                                                                                   >  2:5825801/1‑150 (MQ=255)
    cGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCaa                                                                                                                                                   >  2:7498262/1‑150 (MQ=255)
      aGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGc                                                                                                                                                 <  2:1835938/150‑1 (MQ=255)
        ttCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCt                                                                                                                                               >  2:6068768/1‑150 (MQ=255)
           cTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTAc                                                                                                                                            >  1:8577412/1‑150 (MQ=255)
           cTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTAc                                                                                                                                            <  2:3019948/150‑1 (MQ=255)
             gTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTc                                                                                                                                          >  1:1897390/1‑150 (MQ=255)
                cGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCggt                                                                                                                                       <  1:4046656/150‑1 (MQ=255)
                cGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCggt                                                                                                                                       <  2:2640219/150‑1 (MQ=255)
                 gACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCggtg                                                                                                                                      <  2:2691074/150‑1 (MQ=255)
                    acCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGgtggtg                                                                                                                                   >  2:5527061/1‑150 (MQ=255)
                      cGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGtgct                                                                                                                                 >  1:6304960/1‑150 (MQ=255)
                       gAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGtgctg                                                                                                                                >  1:1441427/1‑150 (MQ=255)
                         agagGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTgctgct                                                                                                                              <  2:2264488/150‑1 (MQ=255)
                           agGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTgg                                                                                                                            <  2:8577412/150‑1 (MQ=255)
                            ggCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGa                                                                                                                           <  1:4691432/150‑1 (MQ=255)
                             gCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGAc                                                                                                                          >  2:4673659/1‑150 (MQ=255)
                                 tGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAgg                                                                                                                      <  2:8894103/150‑1 (MQ=255)
                                 tGGTGCGCAACGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAgg                                                                                                                      <  2:8347083/150‑1 (MQ=255)
                                   gTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTg                                                                                                                    <  2:5317150/150‑1 (MQ=255)
                                    tGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTgg                                                                                                                   >  2:7629523/1‑150 (MQ=255)
                                       acATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGaaa                                                                                                                <  2:8637448/149‑1 (MQ=255)
                                         aTCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGaaaaa                                                                                                              <  2:8154561/150‑1 (MQ=255)
                                          tCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAg                                                                                                             <  1:7629523/150‑1 (MQ=255)
                                               aTGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCAc                                                                                                        >  2:2279442/1‑150 (MQ=255)
                                                 gTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCAccc                                                                                                      <  1:4102836/150‑1 (MQ=255)
                                                  tCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCg                                                                                                     >  2:6715528/1‑150 (MQ=255)
                                                  tCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCg                                                                                                     >  1:357332/1‑150 (MQ=255)
                                                     gAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGAt                                                                                                  <  1:3362127/150‑1 (MQ=255)
                                                     gAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGAt                                                                                                  <  1:1076046/150‑1 (MQ=255)
                                                         tCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATgtgt                                                                                              >  2:2233501/1‑150 (MQ=255)
                                                           aTGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTc                                                                                                                                          <  1:361252/104‑1 (MQ=255)
                                                           aTGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTc                                                                                            <  1:874444/150‑1 (MQ=255)
                                                           aTGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTc                                                                                            >  1:4209830/1‑150 (MQ=255)
                                                             ggAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCaa                                                                                          <  1:8713414/150‑1 (MQ=255)
                                                                 aaGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTg                                                                                      >  1:2614323/1‑150 (MQ=255)
                                                                    cACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTg                                                                                   >  1:7212950/1‑150 (MQ=255)
                                                                        ccGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGc                                                                               <  2:5590934/150‑1 (MQ=255)
                                                                             gcgcgcCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTg                                                                          <  2:3750873/150‑1 (MQ=255)
                                                                                cgtcTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCtt                                                                       <  2:1091504/147‑2 (MQ=255)
                                                                                  tcTGATTGGTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGa                                                                     >  2:3416529/2‑150 (MQ=255)
                                                                                  ccTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGa                                                                     <  1:5388745/150‑1 (MQ=255)
                                                                                  ccTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGa                                                                     <  2:67909/150‑1 (MQ=255)
                                                                                   cTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGaa                                                                    <  2:7849356/150‑1 (MQ=255)
                                                                                   cTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGaa                                                                    <  1:2369024/150‑1 (MQ=255)
                                                                                    tGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAg                                                                   >  2:4894101/1‑150 (MQ=255)
                                                                                     gATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGa                                                                  >  1:905795/1‑150 (MQ=255)
                                                                                      aTCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGAc                                                                 <  1:5560604/150‑1 (MQ=255)
                                                                                       ttggTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACg                                                                <  1:7392211/148‑1 (MQ=255)
                                                                                        cGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACgg                                                               <  2:5338503/150‑1 (MQ=255)
                                                                                         ggagCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTAAACGTGCTGTTGCAGGTGCTGGAAGACGgc                                                              >  1:2783514/4‑150 (MQ=255)
                                                                                          gTGCACCTCCCGGTTATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACggtc                                                             <  2:5980784/150‑3 (MQ=255)
                                                                                          gTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgcc                                                             >  2:4048676/1‑150 (MQ=255)
                                                                                            gCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGgccgc                                                           <  2:1751409/150‑1 (MQ=255)
                                                                                                caccaGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGa                                                       >  1:8511071/1‑150 (MQ=255)
                                                                                                 tcccggttATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                                      >  2:2179133/8‑150 (MQ=255)
                                                                                                 tcccggttATGTCGGCTATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                                      >  1:508177/8‑150 (MQ=255)
                                                                                                 tcccggttATGTCGGCTATGAAGAGGGTGGCGTTCTGAAGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                                      >  2:7455129/8‑150 (MQ=255)
                                                                                                 accaGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGAc                                                      >  1:4592255/1‑150 (MQ=255)
                                                                                                 accaGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGAc                                                      <  1:2094474/150‑1 (MQ=255)
                                                                                                  ccaGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGAcc                                                     <  2:2963930/150‑1 (MQ=255)
                                                                                                   caGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCg                                                    >  1:8079580/1‑150 (MQ=255)
                                                                                                     ggCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGAc                                                  >  1:1260647/1‑150 (MQ=255)
                                                                                                      gCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACa                                                 >  1:2669523/1‑150 (MQ=255)
                                                                                                            tAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACg                                           >  2:8823311/1‑150 (MQ=255)
                                                                                                             aGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCTGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACgg                                          <  2:342868/150‑1 (MQ=255)
                                                                                                              gggTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACggg                                        >  2:4413209/1‑150 (MQ=255)
                                                                                                              ggctATGAAGAGGGTGGCGTTCTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGAc                                                      <  1:7825273/134‑1 (MQ=255)
                                                                                                                gTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACCGCCACGGgcg                                      >  2:3900902/1‑150 (MQ=255)
                                                                                                                gTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGgcg                                      >  2:4725646/1‑150 (MQ=255)
                                                                                                                  aTGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCa                                    <  1:6805359/150‑1 (MQ=255)
                                                                                                                   tGAAGAGGTCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCAc                                   >  1:8197970/1‑150 (MQ=255)
                                                                                                                   tGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCAc                                   >  1:1728502/1‑150 (MQ=255)
                                                                                                                     aaGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCg                                 >  2:8734498/1‑150 (MQ=255)
                                                                                                                       gagGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTg                               <  1:6434251/150‑1 (MQ=255)
                                                                                                                       gagGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTg                               >  1:3684047/1‑150 (MQ=255)
                                                                                                                        agGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGTGCGCCCGTggc                             <  1:1520511/150‑2 (MQ=255)
                                                                                                                         gggCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGa                             <  2:360506/150‑1 (MQ=255)
                                                                                                                          ggCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGTGCGCCCGTggct                            <  2:5647422/150‑3 (MQ=255)
                                                                                                                             ggTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTc                         >  2:1664258/1‑150 (MQ=255)
                                                                                                                              gcgttcTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTcc                        <  2:3800178/145‑1 (MQ=255)
                                                                                                                              gTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTcc                        <  1:5915628/150‑1 (MQ=255)
                                                                                                                               ttACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCg                       >  2:4451493/1‑150 (MQ=255)
                                                                                                                               cgttcTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTCCg                       <  2:2001487/146‑1 (MQ=255)
                                                                                                                                tACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGc                      <  1:4420150/150‑1 (MQ=255)
                                                                                                                                  tcTGACGGAAGCTGTGCGTCGCAGACCATATCAGGTGATCCTGTTCGATGAGATCGAAAAGGCGCATCCGGATGTGTTCAACGTGCTGCTGCAGGTTCTTGATGACGGTCGTCTGACCGACGGTCAGGGCCGTACGGTCGATTTCCGCaa                    >  1:5331113/2‑150 (MQ=255)
                                                                                                                                  ccTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCaa                    <  2:357332/150‑1 (MQ=255)
                                                                                                                                    tGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAaca                  <  1:6653072/150‑1 (MQ=255)
                                                                                                                                      aCCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACAcc                >  1:2244064/1‑150 (MQ=255)
                                                                                                                                        cGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgt              >  2:712285/1‑150 (MQ=255)
                                                                                                                                        cGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgt              >  2:4999142/1‑150 (MQ=255)
                                                                                                                                         gAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtg             >  2:7194251/1‑150 (MQ=255)
                                                                                                                                         gAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtg             <  2:7748346/150‑1 (MQ=255)
                                                                                                                                         gAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtg             >  1:2508804/1‑150 (MQ=255)
                                                                                                                                                tGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCTGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACcgtgatcgtga      <  1:4978725/150‑1 (MQ=255)
                                                                                                                                                   ggcgCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGtgatga   >  1:311931/1‑150 (MQ=255)
                                                                                                                                                   ggcgCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGtgatga   <  2:8671526/150‑1 (MQ=255)
                                                                                                                                                    gcgcAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGCCTTCCGCAACACCGTGATCGTGATGAc  <  2:1465243/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TGGCCGAGTTCCTGTTCGACACCGAAGAGGCGATGGTGCGCATCGACATGTCCGAATTCATGGAGAAGCACTCCGTCGCGCGCCTGATCGGTGCACCACCAGGCTATGTAGGGTATGAAGAGGGCGGTTACCTGACCGAGGCCGTGCGGCGCAAGCCTTACTCGGTGGTGCTGCTGGACGAGGTGGAAAAAGCCCACCCGGATGTGTTCAACGTGCTGTTGCAGGTGCTGGAAGACGGCCGCCTGACCGACAGCCACGGGCGCACCGTGGACTTCCGCAACACCGTGATCGTGATGAC  >  NC_002947/732316‑732613

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