Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,306,729 C→G 7.4% intergenic (+17/+15) cysS → / ← PP_2906 cysteine‑‑tRNA ligase/sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,306,7290CG 7.4% 471.1 / 20.2 145intergenic (+17/+15)cysS/PP_2906cysteine‑‑tRNA ligase/sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (72/62);  new base G (4/7);  total (76/69)
Fisher's exact test for biased strand distribution p-value = 3.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

GATCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCGGC  >  NC_002947/3306580‑3306878
                                                                                                                                                     |                                                                                                                                                     
gATCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCAc                                                                                                                                                       >  1:842211/1‑150 (MQ=255)
  tCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCg                                                                                                                                                     >  1:1723610/1‑150 (MQ=255)
    cgcacgcgcCTTCAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCa                                                                                                                                                   >  2:4790745/5‑150 (MQ=255)
    cAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCa                                                                                                                                                   >  1:1742182/1‑150 (MQ=255)
     aGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGCAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAg                                                                                                                                                  >  1:2185472/1‑150 (MQ=255)
      ggCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGc                                                                                                                                                 <  2:891846/150‑1 (MQ=255)
          cgcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGCCCCa                                                                                                                                             <  2:877070/150‑1 (MQ=255)
           gcCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAg                                                                                                                                            >  1:3198773/1‑150 (MQ=255)
            ccTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCg                                                                                                                                                     <  2:5154460/140‑1 (MQ=255)
            ccTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGc                                                                                                                                           <  2:650678/150‑1 (MQ=255)
             cTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCGGTGCAGCCCCAgtt                                                                                                                                          <  1:4381853/150‑3 (MQ=255)
             cTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCt                                                                                                                                          >  2:5028420/1‑150 (MQ=255)
               gAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAGCTGGGGCTGCGGTGCAGCCCCAgttac                                                                                                                                        <  1:3972690/150‑5 (MQ=255)
                  ggCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTAt                                                                                                                                     <  1:673737/150‑1 (MQ=255)
                   gcgcGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTAtt                                                                                                                                    <  1:5028420/150‑1 (MQ=255)
                      cgTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGc                                                                                                                                 <  2:2690480/150‑1 (MQ=255)
                           ggACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGa                                                                                                                            >  1:2437609/1‑150 (MQ=255)
                           ggACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGa                                                                                                                            >  2:2816959/1‑150 (MQ=255)
                           ggACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGa                                                                                                                            <  1:2236253/150‑1 (MQ=255)
                           ggACAAGAACTGGGCTGAGTCTGACCGCAGTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGa                                                                                                                            <  1:2566231/150‑1 (MQ=255)
                            gACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGAc                                                                                                                           >  2:3069744/1‑150 (MQ=255)
                              cAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGAccc                                                                                                                         <  1:2502798/150‑1 (MQ=255)
                               aaGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCg                                                                                                                        >  1:2758764/1‑150 (MQ=255)
                               aaGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCg                                                                                                                        >  2:234262/1‑150 (MQ=255)
                               aaGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCg                                                                                                                        <  2:1547018/150‑1 (MQ=255)
                                     tGGGCTGAGTCTGACCGCATGCGCGATCAGCTGACCGCCATGGGTGGAGTGCTGGAGGACAGCAAGGGCACGACTACCGGGCGTCTGGCTGACTGAAGGAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCg                                                                                                                  <  2:3244903/150‑1 (MQ=255)
                                      gggCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGCGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCATTTATTTGCTTAGACCCGCAAGCgg                                                                                                                 <  2:2189144/150‑1 (MQ=255)
                                         cTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCa                                                                                                              <  1:5165204/150‑1 (MQ=255)
                                          tGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAg                                                                                                             >  2:1357698/1‑150 (MQ=255)
                                              tCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCa                                                                                                         >  1:3806603/1‑150 (MQ=255)
                                                tGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAAt                                                                                                       >  1:4439473/1‑150 (MQ=255)
                                                    cGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCa                                                                                                   >  2:315161/1‑150 (MQ=255)
                                                      cATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGc                                                                                                 <  2:4233820/150‑1 (MQ=255)
                                                       aTTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGcc                                                                                                >  2:1088148/1‑150 (MQ=255)
                                                         tCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTg                                                                                              >  2:3715840/1‑150 (MQ=255)
                                                           gcgATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGc                                                                                            <  2:4583265/150‑1 (MQ=255)
                                                           gcgATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGc                                                                                            <  2:5248329/150‑1 (MQ=255)
                                                            cgATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGcc                                                                                           >  1:3080027/1‑150 (MQ=255)
                                                               tCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCccca                                                                                        <  1:5373114/150‑1 (MQ=255)
                                                                  gCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTTACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCcccaccc                                                                                     <  1:1939476/150‑1 (MQ=255)
                                                                     gACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATa                                                                                                      <  2:1828363/130‑1 (MQ=255)
                                                                     gACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATa                                                                                                      >  1:1828363/1‑130 (MQ=255)
                                                                      aCCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACt                                                                                 >  2:1504071/1‑150 (MQ=255)
                                                                       ccGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTc                                                                                <  2:4319675/150‑1 (MQ=255)
                                                                        cGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCg                                                                               <  2:2794921/150‑1 (MQ=255)
                                                                        cGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCg                                                                               <  1:4825332/150‑1 (MQ=255)
                                                                        cGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCg                                                                               <  2:2794942/150‑1 (MQ=255)
                                                                         gCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGc                                                                              <  1:411379/150‑1 (MQ=255)
                                                                           cATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCtt                                                                            <  2:2383138/150‑1 (MQ=255)
                                                                           cATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCtt                                                                            <  1:3744578/150‑1 (MQ=255)
                                                                           cATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCtt                                                                            >  2:5264485/1‑150 (MQ=255)
                                                                            aTGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTc                                                                           >  2:185113/1‑150 (MQ=255)
                                                                                    aGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAtc                                                                   >  1:4182558/1‑150 (MQ=255)
                                                                                       gCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCACAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTcc                                                                >  1:3573124/1‑150 (MQ=255)
                                                                                         tGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAg                                                              >  2:1173161/1‑150 (MQ=255)
                                                                                          ggaggaCAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGc                                                             <  1:3472515/150‑1 (MQ=255)
                                                                                                cAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCATTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGaaa                                                       <  1:3138108/150‑1 (MQ=255)
                                                                                                 aGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAAt                                                                                                       >  1:4665382/1‑101 (MQ=255)
                                                                                                 aGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAAt                                                                                                       <  2:4665382/101‑1 (MQ=255)
                                                                                                 aGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAAc                                                      >  1:2718930/1‑150 (MQ=255)
                                                                                                  gCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAAcc                                                     >  1:4135122/1‑150 (MQ=255)
                                                                                                   cAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCa                                                    >  2:4917349/1‑150 (MQ=255)
                                                                                                      gggCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGc                                                 <  1:1799559/150‑1 (MQ=255)
                                                                                                       ggCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCa                                                >  1:3518365/1‑150 (MQ=255)
                                                                                                            gACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTcg                                           >  2:5381140/1‑150 (MQ=255)
                                                                                                            gACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTcg                                           >  2:5381181/1‑150 (MQ=255)
                                                                                                            gACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTcg                                           >  2:2658502/1‑150 (MQ=255)
                                                                                                                 ccTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGAt                                      >  1:4010430/1‑150 (MQ=255)
                                                                                                                  cTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGc                                                                                                                                 <  1:3123041/58‑1 (MQ=255)
                                                                                                                  cTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGc                                                                                                                                 >  2:3123041/1‑58 (MQ=255)
                                                                                                                  cTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATc                                     >  2:4030305/1‑150 (MQ=255)
                                                                                                                  cTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATc                                     >  2:3434350/1‑150 (MQ=255)
                                                                                                                    ggCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTc                                                                                >  1:1270343/1‑105 (MQ=255)
                                                                                                                    ggCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTc                                                                                <  2:1270343/105‑1 (MQ=255)
                                                                                                                     gCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATcgcc                                  >  2:3172705/1‑150 (MQ=255)
                                                                                                                           ggCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACg                            >  1:366914/1‑150 (MQ=255)
                                                                                                                              tgactgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATc                         <  1:1900905/150‑1 (MQ=255)
                                                                                                                              tgactgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATc                         <  1:4917349/150‑1 (MQ=255)
                                                                                                                              tgactgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATc                         <  1:3434350/150‑1 (MQ=255)
                                                                                                                                actgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGa                       <  2:3045300/150‑1 (MQ=255)
                                                                                                                                  tgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAt                                                                    <  1:3092376/103‑1 (MQ=255)
                                                                                                                                  tgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGAt                          <  2:2878385/145‑1 (MQ=255)
                                                                                                                                  tgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACa                     >  2:1206018/1‑150 (MQ=255)
                                                                                                                                  tgaAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACa                     >  2:1726609/1‑150 (MQ=255)
                                                                                                                                  tgaAGTAACTGGGGCGGCACCGCAGCCCCAGCGTTTAGTTGCTTTGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCCCCCACCTCGCTTCCCCCAATCTCCGGCCCGAAACCCTGCAGGGCGACGCTCGCCGCCACGAGCGACa                     >  2:1726868/1‑150 (MQ=255)
                                                                                                                                    aaGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAAc                   >  2:1523365/1‑150 (MQ=255)
                                                                                                                                     aGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAAcc                  <  1:241680/150‑1 (MQ=255)
                                                                                                                                      gTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAAccc                 <  1:4993452/150‑1 (MQ=255)
                                                                                                                                      gTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAAccc                 >  2:3199085/1‑150 (MQ=255)
                                                                                                                                        agcTGGGGCTGCGGTGCAGCCCCAGTTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATc                         >  1:5425694/3‑140 (MQ=255)
                                                                                                                                        agcTGGGGCTGCGGTGCAGCCCCAGTTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATc                         >  1:5425975/3‑140 (MQ=255)
                                                                                                                                        agcTGGGGCTGCGGTGCAGCCCCAGTTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATc                         <  2:5425975/138‑1 (MQ=255)
                                                                                                                                           tGGGGCTGCGGTGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAAGCTCCAGCCCGAAACCATGCAGGTCGACGAGCGCCGCCCCGAGCGACAACCCCAGcc            >  2:4432260/1‑150 (MQ=255)
                                                                                                                                           tGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATc                         <  1:4470449/137‑1 (MQ=255)
                                                                                                                                            ggggCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGccc           >  2:1106666/1‑150 (MQ=255)
                                                                                                                                             gggCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCg          >  2:254413/1‑150 (MQ=255)
                                                                                                                                               gCTGCACCGCAGCCCCAGCTTTTCTTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGACACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCCGCCCGaa        <  2:2545860/150‑1 (MQ=255)
                                                                                                                                               gCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGaa        <  2:2545814/150‑1 (MQ=255)
                                                                                                                                                 tGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGc                                                 <  1:1114761/107‑1 (MQ=255)
                                                                                                                                                   cACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCg    >  1:1060018/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATCCAGGCGCGCCTGAAGGCGCGTGCGGACAAGAACTGGGCTGAGTCTGACCGCATTCGCGATCAGCTGACCGCCATGGGTGTAGTGCTGGAGGACAGCAAGGGCACGACTACCTGGCGTCTGGCTGACTGAAGTAACTGGGGCTGCACCGCAGCCCCAGCTTTTATTTGCTTAGACCCGCAAGCGGGCAGTGCAATACCAGCCTGGCCCCACCCAACTCGCTTCCCCCAATCTCCAGCCCGAAACCATGCAGGTCGACGATCGCCGCCACGATCGACAACCCCAGCCCGAAGCCGGC  >  NC_002947/3306580‑3306878

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: