Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,541,800 A→T 5.6% D149V (GAC→GTC)  PP_3125 → Cro/CI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,541,8000AT 5.6% 1158.3 / 37.0 339D149V (GAC→GTC) PP_3125Cro/CI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (169/151);  new base T (8/11);  total (177/162)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 343 total aligned reads displayed.

GCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGTGAAACGCCA  >  NC_002947/3541651‑3541949
                                                                                                                                                     |                                                                                                                                                     
gcgcAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGa                                                                                                                                                       <  2:8744784/150‑1 (MQ=255)
    aGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGTCGcc                                                                                                                                                   <  1:29166/150‑1 (MQ=255)
      ccTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCt                                                                                                                                                 >  2:6713137/1‑150 (MQ=255)
      ccTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCt                                                                                                                                                 >  2:3608618/1‑150 (MQ=255)
       ctgctgCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGTCGCCCTg                                                                                                                                                <  1:2987185/150‑1 (MQ=255)
       ctgctgCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGTCGCCCTg                                                                                                                                                <  1:243515/150‑1 (MQ=255)
           tgcCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGTCGCCCTGctcc                                                                                                                                            <  2:726627/150‑4 (MQ=255)
           tgcCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCAtt                                                                                                                                            <  1:1073811/150‑1 (MQ=255)
            gcCGTTTCTGATTCAGCCGCCGTCCGAGTACAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGTCGCCCTGctcca                                                                                                                                           <  2:2502551/150‑5 (MQ=255)
             ccGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTc                                                                                                                                          <  1:6534208/150‑1 (MQ=255)
             ccGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTc                                                                                                                                          <  1:5827028/150‑1 (MQ=255)
              cGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCa                                                                                                                                         <  1:3429571/150‑1 (MQ=255)
               gTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCaa                                                                                                                                        <  1:4236025/150‑1 (MQ=255)
                tttCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAAc                                                                                                                                       <  2:6518084/150‑1 (MQ=255)
                tttCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAAc                                                                                                                                       >  2:8522331/1‑150 (MQ=255)
                  tCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGc                                                                                                                                     >  1:856293/1‑150 (MQ=255)
                  tCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGc                                                                                                                                     >  2:3540224/1‑150 (MQ=255)
                   cTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGcc                                                                                                                                    >  1:2621999/1‑150 (MQ=255)
                    tGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGccc                                                                                                                                   >  1:3777703/1‑150 (MQ=255)
                    tGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGccc                                                                                                                                   >  1:3179288/1‑150 (MQ=255)
                           gccgccGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGAcacc                                                                                                                            <  1:7329012/150‑1 (MQ=255)
                           gccgccGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGAcacc                                                                                                                            >  1:6683820/1‑150 (MQ=255)
                           gccgccGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGAcacc                                                                                                                            >  2:4239565/1‑150 (MQ=255)
                                        ttCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTgcgc                                                                                                               >  1:1754608/1‑150 (MQ=255)
                                         tCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCt                                                                                                              >  2:7773196/1‑150 (MQ=255)
                                          cAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTc                                                                                                             <  2:5615390/150‑1 (MQ=255)
                                           aGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCg                                                                                                            >  1:6800084/1‑150 (MQ=255)
                                             cGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGt                                                                                                          <  2:6772978/150‑1 (MQ=255)
                                                 ccTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGgt                                                                                                      >  1:3255349/1‑150 (MQ=255)
                                                 ccTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGgt                                                                                                      <  2:4592274/150‑1 (MQ=255)
                                                 ccTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGgt                                                                                                      <  2:6063786/150‑1 (MQ=255)
                                                   tACGTTCAAGGAACACCAGGGCGAAGAGGTCCTGTTTGTTCACACAGGCCAGGTGGACGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGgtgt                                                                                                    <  2:3027735/150‑1 (MQ=255)
                                                        tCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGc                                                                                               >  1:5999265/1‑150 (MQ=255)
                                                          aaGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCg                                                                                             >  2:8002853/1‑150 (MQ=255)
                                                          aaGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGAAGCCg                                                                                             >  1:2801096/1‑150 (MQ=255)
                                                           aGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCgg                                                                                            <  1:3033538/150‑1 (MQ=255)
                                                           aGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCgg                                                                                            >  2:2031963/1‑150 (MQ=255)
                                                            ggAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGgc                                                                                           >  1:7365605/1‑150 (MQ=255)
                                                            ggAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTAGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGgc                                                                                           >  2:1496211/1‑150 (MQ=255)
                                                                cacCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAg                                                                                       <  2:7795453/150‑1 (MQ=255)
                                                                 acCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGt                                                                                      >  2:2198692/1‑150 (MQ=255)
                                                                      ggCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTg                                                                                 <  2:2805665/150‑1 (MQ=255)
                                                                      ggCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTg                                                                                 <  1:1533943/150‑1 (MQ=255)
                                                                          aaGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttg                                                                             <  2:1977462/150‑1 (MQ=255)
                                                                            gagTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgtt                                                                           <  2:7095350/150‑1 (MQ=255)
                                                                            gagTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgtt                                                                           >  1:5028515/1‑150 (MQ=255)
                                                                             agTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgttg                                                                          <  1:1931490/150‑1 (MQ=255)
                                                                             agTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgttg                                                                          <  1:3481596/150‑1 (MQ=255)
                                                                               ttCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgttgtt                                                                        <  1:5644556/150‑1 (MQ=255)
                                                                               ttCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGgttgttgtt                                                                        <  1:8356119/150‑1 (MQ=255)
                                                                                tCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTc                                                                       <  2:5700105/150‑1 (MQ=255)
                                                                                  cTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTcac                                                                     >  2:694665/1‑150 (MQ=255)
                                                                                  cTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTcac                                                                     >  2:6343733/1‑150 (MQ=255)
                                                                                  cTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTcac                                                                     >  1:2949526/1‑150 (MQ=255)
                                                                                  cTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTcac                                                                     <  1:6474938/150‑1 (MQ=255)
                                                                                   tgtttgttCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTcacc                                                                    >  2:7572300/1‑150 (MQ=255)
                                                                                    gtttgttCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTcacca                                                                   <  1:1674815/150‑1 (MQ=255)
                                                                                       tgttCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGc                                                                <  1:5115061/150‑1 (MQ=255)
                                                                                          tCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGa                                                             <  2:8441876/150‑1 (MQ=255)
                                                                                           cacacaGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGAc                                                            <  1:3884528/150‑1 (MQ=255)
                                                                                           cacacaGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGAc                                                            >  1:4479264/1‑150 (MQ=255)
                                                                                               caGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAAt                                                        <  1:4966062/150‑1 (MQ=255)
                                                                                               caGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAAt                                                        <  2:7474063/150‑1 (MQ=255)
                                                                                                 ggCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGa                                                      >  2:3210715/1‑150 (MQ=255)
                                                                                                      ggTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCAACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAgcg                                                 >  1:6889297/1‑150 (MQ=255)
                                                                                                       gTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAgcgc                                                >  2:582185/1‑150 (MQ=255)
                                                                                                        tGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCa                                               <  1:7773196/150‑1 (MQ=255)
                                                                                                         ggAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCaa                                              >  1:8141203/1‑150 (MQ=255)
                                                                                                         ggAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCaa                                              >  1:4931243/1‑150 (MQ=255)
                                                                                                         ggAAGAGCATATCATGAGTGAGCGTGTGTATCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGCTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCaa                                              >  2:2335632/1‑150 (MQ=255)
                                                                                                          gAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAg                                             >  1:6267861/1‑150 (MQ=255)
                                                                                                           aaGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAgg                                            <  1:2407539/150‑1 (MQ=255)
                                                                                                            aGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAggg                                           <  1:7007399/150‑1 (MQ=255)
                                                                                                             gTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAgggg                                          >  1:8269929/1‑150 (MQ=255)
                                                                                                                  tttCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAAc                                     >  2:3797046/1‑150 (MQ=255)
                                                                                                                    tCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTa                                   >  1:1003430/1‑150 (MQ=255)
                                                                                                                     cATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTAt                                  >  1:5465530/1‑150 (MQ=255)
                                                                                                                     cATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTAt                                  <  2:8383630/150‑1 (MQ=255)
                                                                                                                     cATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCACACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTAt                                  >  1:8432871/1‑150 (MQ=255)
                                                                                                                      aTGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATg                                 >  2:5117418/1‑150 (MQ=255)
                                                                                                                      aTGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATg                                 >  1:536054/1‑150 (MQ=255)
                                                                                                                      aTGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATg                                 >  1:7927611/1‑150 (MQ=255)
                                                                                                                           tgagCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGccccc                            >  2:106792/1‑150 (MQ=255)
                                                                                                                           tgagCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGccccc                            <  1:272263/150‑1 (MQ=255)
                                                                                                                              gCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCcac                         >  1:2306777/1‑150 (MQ=255)
                                                                                                                               cGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCcaca                        >  2:8176660/1‑150 (MQ=255)
                                                                                                                               cGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCcaca                        >  2:4499946/1‑150 (MQ=255)
                                                                                                                                gggTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCacac                       <  1:7571335/150‑1 (MQ=255)
                                                                                                                                 ggTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACg                      >  2:1316892/1‑150 (MQ=255)
                                                                                                                                    gtTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTtc                   >  2:5361459/1‑150 (MQ=255)
                                                                                                                                       gCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTgg                                                                                <  2:3895553/86‑1 (MQ=255)
                                                                                                                                       gCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCGCCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCa                >  1:6131036/1‑150 (MQ=255)
                                                                                                                                        gaaatgCAGGGCGTCGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAg               >  1:1257688/6‑150 (MQ=255)
                                                                                                                                         tGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAgg              >  2:4959318/1‑150 (MQ=255)
                                                                                                                                         tGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAgg              >  2:6627977/1‑150 (MQ=255)
                                                                                                                                           atgCAGGGCGTCGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCg            >  2:6173965/3‑150 (MQ=255)
                                                                                                                                            tgCAGGGCGTCGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGt           >  1:5741215/2‑150 (MQ=255)
                                                                                                                                              cAGGGCGACGCCCTGCATTTCAACGCCCATACACCGCACCGTCTGCGCTCGGTATGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGTGa         <  1:394158/150‑1 (MQ=255)
                                                                                                                                                 ggCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGTGAAAc      <  2:1136070/150‑1 (MQ=255)
                                                                                                                                                 ggCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGTGAAAc      >  2:409513/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGCAGCCTGCTGCCGTTTCTGATTCAGCCGCCGTCCGAGTTCAGCGACCCTACGTTCAAGGAACACCAGGGCGAAGAGTTCCTGTTTGTTCACACAGGCCAGGTGGAAGTGGATTTCATGAGTGAGCGGGTGTTGCTGGAGCAGGGCGACGCCCTGCATTTCAACGCCCAGACACCGCACCGGCTGCGCTCGGTAGGTGTGACGCCGGCGCAGTTACTGGTTGTTGTTCACCACGCAGACGAATGAAAGCGCAAGGGGGCAACTATGCCCCCCACACGTTCTCAGGCGTGAAACGCCA  >  NC_002947/3541651‑3541949

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: