Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 542,010 G→A 5.6% V83I (GTA→ATA) ‡ rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947542,0100GA 5.6% 676.9 / 18.1 215V83I (GTA→ATA) ‡rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (128/75);  new base A (7/5);  total (135/80)
Fisher's exact test for biased strand distribution p-value = 7.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 220 total aligned reads displayed.

GGTCGTTCGGTGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGA  >  NC_002947/541861‑542154
                                                                                                                                                     |                                                                                                                                                
    gTTCGGTGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAAt                                                                                                                                              >  1:3574106/1‑150 (MQ=255)
    gTTCGGTGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAAt                                                                                                                                              >  2:3061594/1‑150 (MQ=255)
        ggTGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGc                                                                                                                                          >  1:4320298/1‑150 (MQ=255)
          tGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCga                                                                                                                                        <  2:694063/150‑1 (MQ=255)
          tGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCga                                                                                                                                        <  2:637418/150‑1 (MQ=255)
          tGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCga                                                                                                                                        <  2:498210/150‑1 (MQ=255)
           gAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCgag                                                                                                                                       >  2:4723956/1‑150 (MQ=255)
               ttAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGt                                                                                                                                   >  2:2523413/1‑150 (MQ=255)
                   gaagCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTgcgg                                                                                                                               >  1:2094814/1‑150 (MQ=255)
                   gaagCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTgcgg                                                                                                                               >  2:4137140/1‑150 (MQ=255)
                     agCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTgcggcg                                                                                                                             >  1:433049/1‑150 (MQ=255)
                      gCCGGAGACGATCAACTATCGTACGTTCAAGCCTGAACGCGATGGTCTCTTCTGCGCGCGTATCTTTGGCCCGATCAAGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGt                                                                                                                            >  2:1004642/1‑150 (MQ=255)
                      gCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGt                                                                                                                            >  1:2809999/1‑150 (MQ=255)
                       ccGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGtt                                                                                                                           <  2:1076013/150‑1 (MQ=255)
                         ggAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGa                                                                                                                         >  2:389619/1‑150 (MQ=255)
                          gAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGaa                                                                                                                        <  1:4137140/150‑1 (MQ=255)
                          gAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGaa                                                                                                                        <  1:394584/150‑1 (MQ=255)
                          gAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGaa                                                                                                                        >  1:2452419/1‑150 (MQ=255)
                             aCCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGtt                                                                                                                     <  2:4167378/150‑1 (MQ=255)
                             aCCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGtt                                                                                                                     <  2:1681897/150‑1 (MQ=255)
                               gaTCAACTATCGTACGTTCAAGCCTGAACGCGATGGTCTCTTCTGCGCGCGTATCTTTGGCCCGATCAAGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGtcac                                                                                                                   <  2:2311651/149‑4 (MQ=255)
                               cATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGc                                                                                                                   <  1:1870162/150‑1 (MQ=255)
                               cATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGc                                                                                                                   >  2:4281109/1‑150 (MQ=255)
                               cATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGc                                                                                                                   >  1:3727140/1‑150 (MQ=255)
                                aTCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGcc                                                                                                                  >  2:4112756/1‑150 (MQ=255)
                                      tACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGCGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCg                                                                                                            >  2:2948801/1‑150 (MQ=255)
                                           catcagcAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGt                                                                                                       >  2:3350308/7‑150 (MQ=255)
                                                tCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtc                                                                                                  >  1:4302273/1‑150 (MQ=255)
                                                tCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtc                                                                                                  <  1:3073476/150‑1 (MQ=255)
                                                 cAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtcg                                                                                                 >  2:3291651/1‑150 (MQ=255)
                                                 cAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTtcgtcg                                                                                                 >  1:2420117/1‑150 (MQ=255)
                                                  aaGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTcgtcgt                                                                                                <  1:436264/150‑1 (MQ=255)
                                                  aaGCCTGAACGCGATGGTCTCTTCTGCGCGCGTATCTTTGGCCCGATCAAGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTcgtcgt                                                                                                >  1:2771896/1‑150 (MQ=255)
                                                    gCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGa                                                                                              <  1:3224855/150‑1 (MQ=255)
                                                    gCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGa                                                                                              >  1:1089079/1‑150 (MQ=255)
                                                    gCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGa                                                                                              <  1:3105031/150‑1 (MQ=255)
                                                      cTGAACGCGATGGTCTCTTCTGCGCGCGTATCTTTGGCCCGATCAAGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAgc                                                                                            >  2:3603291/1‑150 (MQ=255)
                                                       tGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAgcg                                                                                           >  1:3890132/1‑150 (MQ=255)
                                                       tGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAgcg                                                                                           >  2:2149021/1‑150 (MQ=255)
                                                         aGCGTGACGGCCTGATCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTAAGAGTGACTGAGCGGCAAGTACAAGCGCCTCAACCACCGCGGCGAAATCTGCTAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGATCGTCGTGAGCGCa                                                                                         >  2:319419/1‑150 (MQ=255)
                                                                     tGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCg                                                                             >  1:4848246/1‑150 (MQ=255)
                                                                      gTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCCCCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGa                                                                            >  1:2854351/1‑150 (MQ=255)
                                                                       ttCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACTAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAg                                                                           <  1:4969728/150‑1 (MQ=255)
                                                                       ttCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAg                                                                           >  1:43284/1‑150 (MQ=255)
                                                                         cTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCt                                                                         <  2:1916581/150‑1 (MQ=255)
                                                                          tGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTg                                                                        <  2:618695/150‑1 (MQ=255)
                                                                            cgcgcgtaTCTTTGGCCCGATCAAGGATTATGAATGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGc                                                                      <  2:707491/143‑1 (MQ=255)
                                                                             gcCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGcc                                                                     >  1:4767111/1‑150 (MQ=255)
                                                                              ccAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCt                                                                    >  1:4653644/1‑150 (MQ=255)
                                                                              ccAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCt                                                                    >  2:1258346/1‑150 (MQ=255)
                                                                              ccAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCt                                                                    >  2:3181185/1‑150 (MQ=255)
                                                                               cAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTc                                                                   >  1:4892234/1‑150 (MQ=255)
                                                                               cAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTc                                                                   >  1:224057/1‑150 (MQ=255)
                                                                                aaGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCg                                                                  >  2:1889494/1‑150 (MQ=255)
                                                                                 aGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGc                                                                 >  1:2395644/1‑150 (MQ=255)
                                                                                 aGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGATGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGc                                                                 <  2:373084/150‑1 (MQ=255)
                                                                                      tttGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGtt                                                            <  1:4214573/150‑1 (MQ=255)
                                                                                      tttGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGtt                                                            >  2:2471811/1‑150 (MQ=255)
                                                                                           cccAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCca                                                       >  2:4575718/1‑150 (MQ=255)
                                                                                           cccAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTAGCCca                                                       >  1:2944669/1‑150 (MQ=255)
                                                                                                  aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                >  2:4885154/1‑150 (MQ=255)
                                                                                                  aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                >  1:1449361/1‑150 (MQ=255)
                                                                                                  aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                >  2:4618259/1‑150 (MQ=255)
                                                                                                  aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                >  1:2190457/1‑150 (MQ=255)
                                                                                                  aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                <  2:1331379/150‑1 (MQ=255)
                                                                                                  aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGACCACATCTgg                                                >  2:4715061/1‑150 (MQ=255)
                                                                                                   aGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGt                                               >  2:2747295/1‑150 (MQ=255)
                                                                                                     gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                             <  1:4878692/150‑1 (MQ=255)
                                                                                                       cTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCt                                           <  2:3504269/150‑1 (MQ=255)
                                                                                                       cTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCt                                           >  1:740884/1‑150 (MQ=255)
                                                                                                        tACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTg                                          >  1:2341981/1‑150 (MQ=255)
                                                                                                         aCGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGa                                         <  2:2439120/150‑1 (MQ=255)
                                                                                                          cGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGaa                                        >  2:1691273/1‑150 (MQ=255)
                                                                                                           gaatGCTTGTGCGGCAAGTACAAGCGTATGAAATACAAAGGCATCATCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAg                                       >  2:1625166/4‑150 (MQ=255)
                                                                                                                ccTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCt                                  <  1:1570223/150‑1 (MQ=255)
                                                                                                                  tgtgCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCACATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGc                                >  1:1518651/1‑150 (MQ=255)
                                                                                                                   gtgCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGcc                               >  1:2470417/1‑150 (MQ=255)
                                                                                                                      cGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTc                            <  1:3035637/150‑1 (MQ=255)
                                                                                                                        gCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCAt                                                                                        <  2:2879627/88‑1 (MQ=255)
                                                                                                                        gCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCAt                                                                                        >  1:2879627/1‑88 (MQ=255)
                                                                                                                        gCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTccc                          <  2:798312/150‑1 (MQ=255)
                                                                                                                         cAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCg                         >  2:1971593/1‑150 (MQ=255)
                                                                                                                             tACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATc                     >  2:2219174/1‑150 (MQ=255)
                                                                                                                             tACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATc                     >  2:3918444/1‑150 (MQ=255)
                                                                                                                              aCAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCg                    <  1:673350/150‑1 (MQ=255)
                                                                                                                               cAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCgg                   <  1:555483/150‑1 (MQ=255)
                                                                                                                                aaGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGc                  <  1:4432747/150‑1 (MQ=255)
                                                                                                                                aaGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGc                  >  1:2283611/1‑150 (MQ=255)
                                                                                                                                aaGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGc                  >  2:4099049/1‑150 (MQ=255)
                                                                                                                                aaGCGACTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGc                  >  2:1944307/1‑150 (MQ=255)
                                                                                                                                    gcCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGc              >  1:1639655/1‑150 (MQ=255)
                                                                                                                                     ccTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCt             >  2:4475822/1‑150 (MQ=255)
                                                                                                                                      cTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTg            >  2:525826/1‑150 (MQ=255)
                                                                                                                                       tCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGa           <  1:4549246/150‑1 (MQ=255)
                                                                                                                                        cAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGAt          <  2:4320298/150‑1 (MQ=255)
                                                                                                                                         aaGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATg         <  1:3350308/150‑1 (MQ=255)
                                                                                                                                           gCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGa       >  1:4650259/1‑150 (MQ=255)
                                                                                                                                           gCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGa       >  2:3313412/1‑150 (MQ=255)
                                                                                                                                                gcggcgTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGa  >  2:394160/1‑150 (MQ=255)
                                                                                                                                                gcggcgTAATCTGCGAGAAGTGCAGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGa  >  2:1756539/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GGTCGTTCGGTGAAGTTAAGAAGCCGGAAACCATCAACTACCGTACGTTCAAGCCTGAGCGTGACGGCCTGTTCTGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGA  >  NC_002947/541861‑542154

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: