Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 541,501 G→T 12.1% A1292A (GCG→GCT rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947541,5010GT12.1% 488.3 / 45.2 163A1292A (GCG→GCTrpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (70/73);  new base T (10/10);  total (80/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 170 total aligned reads displayed.

GTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTCCTTCAACGTGTTGAT  >  NC_002947/541353‑541650
                                                                                                                                                    |                                                                                                                                                     
gTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGgcgc                                                                                                                                                      >  2:5007926/1‑150 (MQ=255)
   acatgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTgg                                                                                                                                                   <  1:288851/150‑1 (MQ=255)
     atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                                                                                                                 <  1:2550726/150‑1 (MQ=255)
     atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                                                                                                                 <  1:4495470/150‑1 (MQ=255)
     atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                                                                                                                 <  2:1723846/150‑1 (MQ=255)
     atgtacatgCTCAAGCTGAACCACGTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGCTGGAAGTGTGGGCGCTGGaa                                                                                                                                                 <  2:3847936/150‑1 (MQ=255)
      tgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAg                                                                                                                                                >  1:1047128/1‑150 (MQ=255)
         acatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCAt                                                                                                                                             <  2:3081585/150‑1 (MQ=255)
           atgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATAc                                                                                                                                           >  2:3150908/1‑150 (MQ=255)
            tgttgaaGCTCCACCACCTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATatg                                                                                                                                          <  2:4704941/145‑3 (MQ=255)
                 aaGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGGTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGcgcg                                                                                                                                     <  2:3112184/150‑1 (MQ=255)
                  agctccaccaccTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCgg                                                                                                                                    >  1:1759016/7‑150 (MQ=255)
                  aGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCgg                                                                                                                                    <  2:2927816/150‑1 (MQ=255)
                   gCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGc                                                                                                                                   >  1:1792976/1‑150 (MQ=255)
                     tGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCAt                                                                                                                                 >  1:4035853/1‑150 (MQ=255)
                      gAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATa                                                                                                                                >  2:2837723/1‑150 (MQ=255)
                       caccaccTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTac                                                                                                                               <  1:4652525/149‑1 (MQ=255)
                       caccaccTGGTCGACGACAAGATCCATGCTCGTTAGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTac                                                                                                                               >  2:3786283/2‑150 (MQ=255)
                       aaCCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATac                                                                                                                               >  2:983391/1‑150 (MQ=255)
                           aCTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACAccc                                                                                                                           >  1:1688751/1‑150 (MQ=255)
                            cTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCt                                                                                                                          >  1:3263205/1‑150 (MQ=255)
                             ttGGTGGACGACAACATGCACGAGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTg                                                                                                                         >  2:2029224/1‑150 (MQ=255)
                              tggtcgACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGc                                                                                                                        <  1:5028228/150‑1 (MQ=255)
                              tggtcgACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGc                                                                                                                        <  1:2457769/150‑1 (MQ=255)
                                gtggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCaa                                                                                                                      <  1:2244274/150‑1 (MQ=255)
                                gtggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCaa                                                                                                                      <  1:4815179/150‑1 (MQ=255)
                                   gacgacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGACGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGaa                                                                                                                   >  2:2570233/1‑150 (MQ=255)
                                   gacgacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGGTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGaa                                                                                                                   <  2:4298892/150‑1 (MQ=255)
                                   gacgacAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGaa                                                                                                                   >  1:2840256/1‑150 (MQ=255)
                                      gacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTGGGGGAGATGGATGTGTGGGCGCTGGTAGCATACGGCGCGGCATAGACCCTGGAAGAAATg                                                                                                                >  2:4419005/1‑150 (MQ=255)
                                      gacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATg                                                                                                                >  2:5074579/1‑150 (MQ=255)
                                      gacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATg                                                                                                                <  1:417624/150‑1 (MQ=255)
                                            aTGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTcaca                                                                                                          <  2:137723/150‑1 (MQ=255)
                                            aTGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTcaca                                                                                                          <  2:1769069/150‑1 (MQ=255)
                                            aTGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTcaca                                                                                                          <  2:815915/150‑1 (MQ=255)
                                              gCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGGGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGt                                                                                                        >  1:720006/1‑150 (MQ=255)
                                                aCGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGa                                                                                                      >  2:2755234/1‑150 (MQ=255)
                                                 cgcgcgGTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGaa                                                                                                     <  1:4445926/150‑1 (MQ=255)
                                                  gcgcgTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAg                                                                                                    >  1:251692/1‑150 (MQ=255)
                                                   cgcgTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGt                                                                                                   <  2:2409437/150‑1 (MQ=255)
                                                        tCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgac                                                                                              >  1:2100529/1‑150 (MQ=255)
                                                          cACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacga                                                                                            >  1:3057273/1‑150 (MQ=255)
                                                           aCTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacgac                                                                                           >  2:2797303/1‑150 (MQ=255)
                                                           aCTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacgac                                                                                           <  1:4275730/150‑1 (MQ=255)
                                                            cTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGacgacg                                                                                          >  2:587947/1‑150 (MQ=255)
                                                               gTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTgg                                                                                                                                                   <  2:3052674/90‑1 (MQ=255)
                                                               gTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGa                                                                                       >  1:3066885/1‑150 (MQ=255)
                                                               gTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGAGGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGa                                                                                       >  1:1108134/1‑150 (MQ=255)
                                                                   cTACAGCCGGGTTACCCAGCAGCCGCTGGGTGGTACGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACgg                                                                                   <  2:2462606/150‑1 (MQ=255)
                                                                    tACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGc                                                                                  >  1:965683/1‑150 (MQ=255)
                                                                    tACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGc                                                                                  >  1:965712/1‑150 (MQ=255)
                                                                    tACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGc                                                                                  >  1:1781923/1‑150 (MQ=255)
                                                                     aCAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGcc                                                                                 >  2:935055/1‑150 (MQ=255)
                                                                        gCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTa                                                                              >  1:2198217/1‑150 (MQ=255)
                                                                        gCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTa                                                                              >  1:5118634/1‑150 (MQ=255)
                                                                           tcgtcaCCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCa                                                                           >  2:4034007/6‑150 (MQ=255)
                                                                             gtcaCCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAg                                                                         <  2:219636/147‑1 (MQ=255)
                                                                             gtcaCCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAg                                                                         <  2:4531014/147‑1 (MQ=255)
                                                                                  ccAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGTAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTa                                                                    <  1:58804/150‑1 (MQ=255)
                                                                                  ccAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAggtcta                                                                    >  2:2728499/1‑145 (MQ=255)
                                                                                  ccAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAggtcta                                                                    >  2:693311/1‑145 (MQ=255)
                                                                                  ccAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAggtcta                                                                    >  1:3431076/1‑145 (MQ=255)
                                                                                  ccAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAggtcta                                                                    >  1:3085067/1‑145 (MQ=255)
                                                                                    agcagcCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACa                                                                  <  1:2326284/150‑1 (MQ=255)
                                                                                     gcggctgCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACaa                                                                 <  1:3677657/144‑1 (MQ=255)
                                                                                        gcCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAGGTCTAcgaagc                                                              >  1:3788088/1‑145 (MQ=255)
                                                                                          cGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACa                                                            >  1:2304493/1‑150 (MQ=255)
                                                                                          cGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAGGTCTAcgaagcga                                                            <  1:5078235/150‑8 (MQ=255)
                                                                                               ggtggtAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTg                                                       <  2:3057273/150‑1 (MQ=255)
                                                                                                  ggtAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGAt                                                    <  2:1792976/150‑1 (MQ=255)
                                                                                                    tAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCg                                                                                                 <  1:43991/103‑1 (MQ=255)
                                                                                                    tAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCg                                                                                                 >  2:43991/1‑103 (MQ=255)
                                                                                                    tAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATgg                                                  <  1:5015831/150‑1 (MQ=255)
                                                                                                     aaGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGc                                                 <  2:4966894/150‑1 (MQ=255)
                                                                                                      aGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCg                                                >  1:433549/1‑150 (MQ=255)
                                                                                                            aGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCAcc                                          >  2:3547733/1‑150 (MQ=255)
                                                                                                                cGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTAt                                      >  1:2932146/1‑150 (MQ=255)
                                                                                                                      tCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGcc                                >  1:560483/1‑150 (MQ=255)
                                                                                                                        aGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCgg                              <  1:3150908/150‑1 (MQ=255)
                                                                                                                        aGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCgg                              >  2:3423702/1‑150 (MQ=255)
                                                                                                                            tttGGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCGTGGTAGAAATGCTCTCAGTGAAGTCGGACGTCGTGAACGGCCGTAGCAAGATGTACAATAACATGGTGGATGGCGATCACCGTATGGAGCCGGGCAt                          >  2:1181001/1‑150 (MQ=255)
                                                                                                                            tttCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTTGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCAt                          >  1:1002495/1‑150 (MQ=255)
                                                                                                                              tCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGc                        <  1:1776740/150‑1 (MQ=255)
                                                                                                                                ggggAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGccc                      <  1:2256052/150‑1 (MQ=255)
                                                                                                                                ggggAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGccc                      >  1:1367965/1‑150 (MQ=255)
                                                                                                                                ggggAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGccc                      >  1:3528677/1‑150 (MQ=255)
                                                                                                                                  ggAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGa                    <  2:1047128/150‑1 (MQ=255)
                                                                                                                                  ggAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGa                    <  1:4695903/150‑1 (MQ=255)
                                                                                                                                    agaTGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGt                  <  2:4253220/150‑1 (MQ=255)
                                                                                                                                     gaTGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTc                 <  2:2060472/150‑1 (MQ=255)
                                                                                                                                     gaTGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTc                 <  2:1872380/150‑1 (MQ=255)
                                                                                                                                      aTGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTcc                <  2:1688751/150‑1 (MQ=255)
                                                                                                                                      aTGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTcc                >  2:547582/1‑150 (MQ=255)
                                                                                                                                          aaGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCTATCACCGTATGGAGCCGGGCATGCCCGAGTCCTTCa            >  2:4154463/1‑150 (MQ=255)
                                                                                                                                           aGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTCCTTCaa           >  2:1872012/1‑150 (MQ=255)
                                                                                                                                               tgGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCGCAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTCCTTCAACgtg       <  1:3023082/150‑1 (MQ=255)
                                                                                                                                               tgGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGTCGGGCATGCCCGAGTCCTTCAACgtg       <  1:2832103/150‑1 (MQ=255)
                                                                                                                                                  gcgcTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCAGAGTCCTTCAAAttgttg    >  2:1213821/1‑150 (MQ=255)
                                                                                                                                                   cgcTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGAGCATGCCCGAGTCCTTCAACGTGTTGa   <  1:3115631/150‑1 (MQ=255)
                                                                                                                                                    gcTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTCCTTCAACGTGTTGAt  >  2:512260/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTCCTTCAACGTGTTGAT  >  NC_002947/541353‑541650

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: