Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,543,692 T→A 10.1% intergenic (+33/+30) PP_3128 → / ← galE exopolysaccharide biosynthesis/transport protein/UDP‑glucose 4‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,543,6920TA10.1% 409.9 / 29.2 138intergenic (+33/+30)PP_3128/galEexopolysaccharide biosynthesis/transport protein/UDP‑glucose 4‑epimerase
Reads supporting (aligned to +/- strand):  ref base T (71/52);  new base A (7/7);  total (79/59)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 148 total aligned reads displayed.

GCCAACGGAATGAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGCCTTGCTCGCATCCGCCACCAG  >  NC_002947/3543543‑3543841
                                                                                                                                                     |                                                                                                                                                     
gCCAACGGAATGAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCtt                                                                                                                                                       <  2:92354/150‑1 (MQ=255)
gCCAACGGAATGAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCtt                                                                                                                                                       >  1:954632/1‑150 (MQ=255)
gCCAACGGAATGAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCtt                                                                                                                                                       <  1:1621286/150‑1 (MQ=255)
  cAACGGAATGAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTtgt                                                                                                                                                     >  1:450801/1‑150 (MQ=255)
  cAACGGAATGAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTtgt                                                                                                                                                     <  2:4130290/150‑1 (MQ=255)
  cAACGGAATGAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTtgt                                                                                                                                                     >  1:2979245/1‑150 (MQ=255)
    aCGGAATGAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTtgt                                                                                                                                                     >  2:2697426/1‑148 (MQ=255)
    aCGGAATGAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTtgt                                                                                                                                                     <  1:2697426/148‑1 (MQ=255)
        aaTGAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTgcg                                                                                                                                               >  2:5107799/1‑150 (MQ=255)
         aTGAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTgcgc                                                                                                                                              <  2:5269376/150‑1 (MQ=255)
          tGAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGcc                                                                                                                                             <  2:902042/150‑1 (MQ=255)
           gAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCt                                                                                                                                            >  1:2882954/1‑150 (MQ=255)
               aacaaGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGaaa                                                                                                                                        >  2:1359624/1‑150 (MQ=255)
               aacaaGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGaaa                                                                                                                                        >  2:1359815/1‑150 (MQ=255)
               aacaaGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGaaa                                                                                                                                        >  2:4608050/1‑150 (MQ=255)
               aacaaGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGaaa                                                                                                                                        <  2:2577420/150‑1 (MQ=255)
               aacaaGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGaaa                                                                                                                                        <  1:1955027/150‑1 (MQ=255)
                 caaGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGc                                                                                                                                      >  1:1525504/1‑150 (MQ=255)
                  aagaagGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCt                                                                                                                                     >  1:4824432/1‑150 (MQ=255)
                      agGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTTAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGa                                                                                                                                 <  2:3025862/150‑1 (MQ=255)
                       ggAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAg                                                                                                                                >  1:1620680/1‑150 (MQ=255)
                          aaTCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGcc                                                                                                                             >  1:5117671/1‑150 (MQ=255)
                          aaTCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGcc                                                                                                                             <  2:4896793/150‑1 (MQ=255)
                            tCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCaa                                                                                                                           >  2:1732698/1‑150 (MQ=255)
                              gAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACt                                                                                                                         >  1:1510601/1‑150 (MQ=255)
                                      ccTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCt                                                                                                                 <  1:4666696/150‑1 (MQ=255)
                                       cTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTTCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTc                                                                                                                <  1:3480169/150‑1 (MQ=255)
                                       cTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTc                                                                                                                >  1:2331412/1‑150 (MQ=255)
                                        tGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTcc                                                                                                               >  2:446780/1‑150 (MQ=255)
                                        tGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCACAAGGCTGCGCCTGaaa                                                                                                                                        <  2:622164/125‑1 (MQ=39)
                                        tGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCACAAGGCTGCGCCTGaaa                                                                                                                                        >  1:622164/1‑125 (MQ=39)
                                          cATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAg                                                                                                             <  1:3765541/150‑1 (MQ=255)
                                           aTCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCACAAGGCTGCGCCTGa                                                                                                                                          <  1:3877500/120‑1 (MQ=39)
                                              aCCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGAt                                                                                                         >  2:1898746/1‑150 (MQ=255)
                                              aCCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGAt                                                                                                         >  2:2627350/1‑150 (MQ=255)
                                                 tGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACt                                                                                                      <  2:2449963/150‑1 (MQ=255)
                                                  gggTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAgc                                                                                                                     <  1:4240825/134‑1 (MQ=255)
                                                  gggTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAgc                                                                                                                     >  2:4240825/1‑134 (MQ=255)
                                                   ggTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGc                                                                                                    <  2:5186624/150‑1 (MQ=255)
                                                     tACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCACAAGGCTGCGCCTGaaa                                                                                                                                        >  1:3890849/1‑112 (MQ=38)
                                                     tACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCACAAGGCTGCGCCTGaaa                                                                                                                                        <  2:3890849/112‑1 (MQ=38)
                                                     tACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCACAAGGCTGCGCCTGaaa                                                                                                                                        <  2:3890845/112‑1 (MQ=38)
                                                        ggTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTCTGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTc                                                                                               >  2:3522341/1‑150 (MQ=255)
                                                        ggTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTc                                                                                               >  2:3630086/1‑150 (MQ=255)
                                                        ggTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTc                                                                                               >  1:3318715/1‑150 (MQ=255)
                                                         gTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTcc                                                                                              >  1:182286/1‑150 (MQ=255)
                                                         gTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTcc                                                                                              >  1:1661/1‑150 (MQ=255)
                                                          tGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCACCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTccc                                                                                             >  2:4807404/1‑150 (MQ=255)
                                                            aaCCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCAc                                                                                           <  2:2253903/150‑1 (MQ=255)
                                                            aaCCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCAc                                                                                           <  1:3527646/150‑1 (MQ=255)
                                                             aCCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACt                                                                                          <  1:3216314/150‑1 (MQ=255)
                                                                  aTCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCAc                                                                                     >  2:3862564/1‑150 (MQ=255)
                                                                     gggACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTCCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCa                                                                                  >  1:4453848/1‑150 (MQ=255)
                                                                      ggACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCAt                                                                                 >  2:3602620/1‑150 (MQ=255)
                                                                         cAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGcc                                                                              <  1:4838723/150‑1 (MQ=255)
                                                                            tgtTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGca                                                                           >  1:3390638/1‑150 (MQ=255)
                                                                                gAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAAt                                                                       >  2:865982/1‑150 (MQ=255)
                                                                                  aCAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTgcgc                                                                                                                                              <  1:2786889/77‑1 (MQ=255)
                                                                                  aCAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTgcgc                                                                                                                                              >  2:2786889/1‑77 (MQ=255)
                                                                                    aaGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGc                                                                   >  1:5094939/1‑150 (MQ=255)
                                                                                      ggCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCACAAGGCTGCGCCTGa                                                                                                                                          <  2:1905695/77‑1 (MQ=25)
                                                                                       gCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCt                                                                     <  2:3062805/145‑1 (MQ=255)
                                                                                       gCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTcc                                                                >  1:811189/1‑150 (MQ=255)
                                                                                         cgacgaACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAg                                                              >  1:1943790/1‑150 (MQ=255)
                                                                                          gacgaACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGt                                                             >  1:258511/1‑150 (MQ=255)
                                                                                          gacgaACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGt                                                             >  1:4075726/1‑150 (MQ=255)
                                                                                             gaACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACt                                                                                                                         >  2:2815747/1‑87 (MQ=255)
                                                                                               aCAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGc                                                        >  2:115896/1‑150 (MQ=255)
                                                                                               aCAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGc                                                        >  1:3503128/1‑150 (MQ=255)
                                                                                                    cACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATg                                                   >  1:542765/1‑150 (MQ=255)
                                                                                                    cACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATg                                                   <  2:1283469/150‑1 (MQ=255)
                                                                                                      cGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCa                                                 <  1:4001958/150‑1 (MQ=255)
                                                                                                        cACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAgg                                               <  1:4658925/150‑1 (MQ=255)
                                                                                                         aCCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAggg                                              <  1:3543705/150‑1 (MQ=255)
                                                                                                          ccTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGgc                                             <  2:4148204/150‑1 (MQ=255)
                                                                                                           cTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGgcg                                            >  2:4840676/1‑150 (MQ=255)
                                                                                                            tATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGgcgc                                           >  2:3445348/1‑150 (MQ=255)
                                                                                                                 cTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGccagcc                                      >  2:2950882/1‑150 (MQ=255)
                                                                                                                 cTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGccagcc                                      >  1:2363706/1‑150 (MQ=255)
                                                                                                                  tGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCAc                                                                                           <  2:5048416/96‑1 (MQ=255)
                                                                                                                  tGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCAc                                                                                           >  1:5048416/1‑96 (MQ=255)
                                                                                                                   gAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGccagccca                                    <  2:542765/150‑1 (MQ=255)
                                                                                                                   gAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGccagccca                                    <  1:5191583/150‑1 (MQ=255)
                                                                                                                     gCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGccagcccagc                                  <  2:12590/150‑1 (MQ=255)
                                                                                                                        ttCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCAcc                               >  2:1500479/1‑150 (MQ=255)
                                                                                                                                 tACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGcc                      <  2:4588273/150‑1 (MQ=255)
                                                                                                                                  tagCTTTCAGGCGCAGCACAAGGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTTCATGGCCt                     >  1:4538589/3‑150 (MQ=39)
                                                                                                                                  aCGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGCCt                     >  1:2525056/1‑150 (MQ=255)
                                                                                                                                      tttCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCCGCACCTGCATGGCCTAGcg                 >  1:4050909/1‑149 (MQ=255)
                                                                                                                                         cAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGCCTTGCTCGc              <  1:5131805/150‑1 (MQ=255)
                                                                                                                                          aGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGCCTTGCTCGCa             >  1:3129968/1‑150 (MQ=255)
                                                                                                                                           ggCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGCCTTGCTCGCAt            >  2:1102586/1‑150 (MQ=255)
                                                                                                                                            gcgcAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCTGGGCGCCAGCCCAGCACCTGCATGGCCTTGCTCGCagc           >  2:4768398/1‑148 (MQ=255)
                                                                                                                                            gcgcAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGCCTTGCTCGCATc           >  1:133030/1‑150 (MQ=255)
                                                                                                                                            gcgcAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGCCTTGCTCGCATc           <  2:954632/150‑1 (MQ=255)
                                                                                                                                            gcgcAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGCCTTGCTCGCATc           >  1:3321982/1‑150 (MQ=255)
                                                                                                                                             cgcAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGCCTTGCTCGCATcc          <  2:258511/150‑1 (MQ=255)
                                                                                                                                             cgcAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGCCTTGCTCGCATcc          <  1:1359815/150‑1 (MQ=255)
                                                                                                                                             cgcAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGCCTTGCTCGCATcc          <  1:1359624/150‑1 (MQ=255)
                                                                                                                                                     tgtgCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGCCTTGCTCGCATCCGctgcaac  <  2:2887429/150‑7 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAACGGAATGAACAACAAGAAGGAAATCGAGGATTGCCTGCATCACCTGGGTACGGTGAACCTGATCGGGACAGTGTTGAACAAGGCCGACGAACAGCCACGCACCTATTACTGAGCCTTCTGGCAATACGCTTTCAGGCGCAGCCTTGTGCTGCGCCTGAAAGCTATGAGGCCAACTCAGCGCCTCCAGGGATACTGCAACTCCCACTGCCACGCATGCCGCACAATCTGCTCCAGTGAGGCGAATGCAGGGCGCCAGCCCAGCACCTGCATGGCCTTGCTCGCATCCGCCACCAG  >  NC_002947/3543543‑3543841

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: