Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 541,492 A→G 13.5% E1289E (GAA→GAG 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,4920AG13.5% 435.6 / 48.6 174E1289E (GAA→GAGrpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (75/72);  new base G (12/11);  total (91/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 = 2.13e-01

Only 100 of 178 total aligned reads displayed.

GTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTCCTTCAA  >  NC_002947/541343‑541641
                                                                                                                                                     |                                                                                                                                                     
gtgACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGaa                                                                                                                                                       <  1:3760505/150‑1 (MQ=255)
gtgACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGaa                                                                                                                                                       >  2:2285071/1‑150 (MQ=255)
    ccGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAgtgt                                                                                                                                                   <  1:699890/150‑1 (MQ=255)
     cGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGtgtg                                                                                                                                                  >  1:3653853/1‑150 (MQ=255)
     cGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGtgtg                                                                                                                                                  >  2:2256052/1‑150 (MQ=255)
      gTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTgg                                                                                                                                                 >  1:2060472/1‑150 (MQ=255)
          gTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGgcgc                                                                                                                                             >  2:5007926/1‑150 (MQ=255)
           ttACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCt                                                                                                                                            <  2:387088/150‑1 (MQ=255)
             acatgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTgg                                                                                                                                          <  1:288851/150‑1 (MQ=255)
               atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                                                                                                        <  1:4495470/150‑1 (MQ=255)
               atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                                                                                                        <  2:3465683/150‑1 (MQ=255)
               atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                                                                                                        <  2:1149519/150‑1 (MQ=255)
               atgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGaa                                                                                                                                        <  2:1583652/150‑1 (MQ=255)
                tgtacatgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAg                                                                                                                                       >  1:1047128/1‑150 (MQ=255)
                     atgCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATAc                                                                                                                                  >  2:3150908/1‑150 (MQ=255)
                       gCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACgg                                                                                                                                >  2:831023/1‑150 (MQ=255)
                            agctccaccaccTGGTCGACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCgg                                                                                                                           >  1:1759016/7‑150 (MQ=255)
                            aGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCgg                                                                                                                           <  2:2927816/150‑1 (MQ=255)
                              cTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCa                                                                                                                         <  2:4825639/150‑1 (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:2832103/1‑150 (MQ=255)
                                 aaCCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATac                                                                                                                      >  2:983391/1‑150 (MQ=255)
                                        tggtggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGc                                                                                                               <  1:2285071/150‑1 (MQ=255)
                                        tggtcgACGACAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGc                                                                                                               <  1:2457769/150‑1 (MQ=255)
                                          gtggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCaa                                                                                                             <  1:4815179/150‑1 (MQ=255)
                                          gtggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCaa                                                                                                             >  1:1729486/1‑150 (MQ=255)
                                          gtggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCaa                                                                                                             <  1:2244274/150‑1 (MQ=255)
                                          gtggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCaa                                                                                                             <  2:2929148/150‑1 (MQ=255)
                                          gtggACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCaa                                                                                                             <  1:3956370/150‑1 (MQ=255)
                                          gtggACGACAAGATGCACGCGCGGTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCaa                                                                                                             <  2:2657645/150‑1 (MQ=255)
                                             gacgacAAGATCCATGCTCGTTCGATCGGACCTTACTCGCTCGTCACCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGaa                                                                                                          >  1:2840256/1‑150 (MQ=255)
                                              acgacgAGATGCACGCGCGGTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGaaa                                                                                                         <  2:3351055/150‑1 (MQ=255)
                                                gacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTGGGGGAGATGGATGTGTGGGCGCTGGTAGCATACGGCGCGGCATAGACCCTGGAAGAAATg                                                                                                       >  2:4419005/1‑150 (MQ=255)
                                                gacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATg                                                                                                       <  1:417624/150‑1 (MQ=255)
                                                gacAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATg                                                                                                       >  1:216758/1‑150 (MQ=255)
                                                   aaGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTc                                                                                                    >  2:1021289/1‑150 (MQ=255)
                                                      aTGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTcaca                                                                                                 <  1:584151/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:3015501/1‑150 (MQ=255)
                                                            gcgcgTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAg                                                                                           >  2:5015831/1‑150 (MQ=255)
                                                            gcgcgTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAg                                                                                           >  1:251692/1‑150 (MQ=255)
                                                             cgcgTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGt                                                                                          <  2:2409437/150‑1 (MQ=255)
                                                              gcgGTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTc                                                                                         <  2:5096689/150‑1 (MQ=255)
                                                                gTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCgg                                                                                       >  1:3269723/1‑150 (MQ=255)
                                                                  tCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgac                                                                                     >  1:2100529/1‑150 (MQ=255)
                                                                    cACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacga                                                                                   >  1:3057273/1‑150 (MQ=255)
                                                                     aCTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGgacgac                                                                                  <  1:4275730/150‑1 (MQ=255)
                                                                         gTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTgg                                                                                                                                          <  2:3052674/90‑1 (MQ=255)
                                                                         gTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGa                                                                              >  1:3066885/1‑150 (MQ=255)
                                                                         gTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGAGGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGa                                                                              >  1:1108134/1‑150 (MQ=255)
                                                                          ttCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGaa                                                                             <  2:4333/150‑1 (MQ=255)
                                                                             cTACAGCCGGGTTACCCAGCAGCCGCTGGGTGGTACGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACgg                                                                          <  2:2462606/150‑1 (MQ=255)
                                                                              tACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGc                                                                         >  1:1786036/1‑150 (MQ=255)
                                                                              tACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGc                                                                         >  1:965712/1‑150 (MQ=255)
                                                                              tACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGc                                                                         >  1:4794442/1‑150 (MQ=255)
                                                                               aCAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGcc                                                                        >  2:935055/1‑150 (MQ=255)
                                                                               aCAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGcc                                                                        >  2:966703/1‑150 (MQ=255)
                                                                                  gCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGtc                                                                     >  1:2778958/1‑149 (MQ=255)
                                                                                     tcgtcaCCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCa                                                                  >  2:4034007/6‑150 (MQ=255)
                                                                                     tcgtcaCCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCa                                                                  >  2:1780318/6‑150 (MQ=255)
                                                                                       gtcaCCCAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAg                                                                <  2:219636/147‑1 (MQ=255)
                                                                                         caCCCAGCGGCTGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGAt                                                              >  2:1851300/2‑150 (MQ=255)
                                                                                          aCCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATg                                                             >  1:1833671/1‑150 (MQ=255)
                                                                                           cccAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGgcg                                                                                                                                              <  2:4277542/68‑1 (MQ=255)
                                                                                            ccAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAggtcta                                                           >  1:3085067/1‑145 (MQ=255)
                                                                                            ccAGCAGCCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAggtcta                                                           >  1:3431076/1‑145 (MQ=255)
                                                                                             cagcagCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTAc                                                          >  2:4519882/1‑150 (MQ=255)
                                                                                             cagcagCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTAc                                                          >  2:3099034/1‑150 (MQ=255)
                                                                                              agcagcCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACa                                                         <  1:2326284/150‑1 (MQ=255)
                                                                                                  gcCGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAGGTCTAcgaagc                                                     >  1:3788088/1‑145 (MQ=255)
                                                                                                    cGCTCGGTGGTAAGGCGCAGTTCGGTGGTCAGCGCTTCGGCGAAATGGAGGTCTGGGCTCTGGAAGCATATGGTGCGGCTTACACGCTGCAGGAAATGCTTACAGTCAAGTCGGACGACGTTGCAGGCCGTACCAAGGTCTAcgaagcga                                                   <  1:5078175/150‑8 (MQ=255)
                                                                                                              tAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATgg                                         <  1:5015831/150‑1 (MQ=255)
                                                                                                              tAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATgg                                         >  1:3482693/1‑150 (MQ=255)
                                                                                                                aGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCg                                       >  1:433549/1‑150 (MQ=255)
                                                                                                                      aGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGCAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCAcc                                 <  2:4618335/150‑1 (MQ=255)
                                                                                                                       gTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCg                                >  2:3062808/1‑150 (MQ=255)
                                                                                                                             tggtCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGa                          >  1:3193006/1‑150 (MQ=255)
                                                                                                                                tCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGcc                       >  1:560483/1‑150 (MQ=255)
                                                                                                                                 cAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCg                      <  2:1223030/150‑1 (MQ=255)
                                                                                                                                     gTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCa                  <  1:948561/150‑1 (MQ=255)
                                                                                                                                      tttCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTTGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCAt                 >  1:1002495/1‑150 (MQ=255)
                                                                                                                                      tttCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCAt                 >  2:617937/1‑150 (MQ=255)
                                                                                                                                        tCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGc               <  1:1776740/150‑1 (MQ=255)
                                                                                                                                          ggggAGATTGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGccc             <  1:4419005/150‑1 (MQ=255)
                                                                                                                                          ggggAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGccc             <  2:1764660/150‑1 (MQ=255)
                                                                                                                                          ggggAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGccc             >  1:1367965/1‑150 (MQ=255)
                                                                                                                                          ggggAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGccc             >  2:2733544/1‑150 (MQ=255)
                                                                                                                                          ggggAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGccc             <  1:2256052/150‑1 (MQ=255)
                                                                                                                                            ggAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGa           <  1:5100385/150‑1 (MQ=255)
                                                                                                                                             gagaTGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAg          <  2:1974292/150‑1 (MQ=255)
                                                                                                                                             gagaTGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAg          <  2:547765/150‑1 (MQ=255)
                                                                                                                                               gaTGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTc        <  2:1872380/150‑1 (MQ=255)
                                                                                                                                               gaTGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTc        <  2:2060428/150‑1 (MQ=255)
                                                                                                                                                aTGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTcc       <  2:1688751/150‑1 (MQ=255)
                                                                                                                                                    aaGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCTATCACCGTATGGAGCCGGGCATGCCCGAGTCCTTCa   >  2:4154463/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGACCGTTGGTTACATGTACATGCTCAAGCTGAACCACTTGGTGGACGACAAGATGCACGCGCGTTCCACTGGTTCCTACAGCCTGGTTACCCAGCAGCCGCTGGGTGGTAAGGCGCAGTTCGGTGGTCAGCGTTTCGGGGAGATGGAAGTGTGGGCGCTGGAAGCATACGGCGCGGCATACACCCTGCAAGAAATGCTCACAGTGAAGTCGGACGACGTGAACGGCCGTACCAAGATGTACAAGAACATCGTGGATGGCGATCACCGTATGGAGCCGGGCATGCCCGAGTCCTTCAA  >  NC_002947/541343‑541641

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: