Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,341,095 T→A 5.6% intergenic (‑132/+111) PP_2935 ← / ← ybiT membrane protein/ABC‑F family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,341,0950TA 5.6% 799.5 / 27.1 252intergenic (‑132/+111)PP_2935/ybiTmembrane protein/ABC‑F family ATPase
Reads supporting (aligned to +/- strand):  ref base T (87/146);  new base A (2/12);  total (89/163)
Fisher's exact test for biased strand distribution p-value = 9.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

Only 100 of 254 total aligned reads displayed.

TGGTTGTTGCGCGGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTA  >  NC_002947/3340946‑3341244
                                                                                                                                                     |                                                                                                                                                     
gtcgctggccacgcgcatcatcgagctgtcggccgatggtgtggtggacttcagcggtacctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCtt                                                                                                                                                       <  2:8217747‑M1/14‑1 (MQ=255)
  gttgttGCGCGGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGc                                                                                                                                                     <  1:1232119/150‑1 (MQ=255)
  gttgttGCGCGGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGc                                                                                                                                                     <  1:6152837/150‑1 (MQ=255)
  gttgttGCGCGGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGc                                                                                                                                                     <  1:4901871/150‑1 (MQ=255)
  gttgttGCGCGGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGc                                                                                                                                                     <  1:4901802/150‑1 (MQ=255)
  gttgttGCGCGGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGc                                                                                                                                                     <  2:5367865/150‑1 (MQ=255)
   ttgttgCGCGGTCCTGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGcc                                                                                                                                                    <  2:9024596/150‑1 (MQ=255)
   gctggccacgcgcatcatcgcgctgtcggccgatggtgtggtggacttcagcggtacctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGcc                                                                                                                                                    >  2:8148743‑M1/134‑150 (MQ=255)
     tggccacgcgcatcatcgagctgtcggccgatggtgtggtggacttcagcggtacctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCg                                                                                                                                                  <  1:1393110‑M1/19‑1 (MQ=255)
     gttgCGCGGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCAAGCCCg                                                                                                                                                  <  2:4230907/150‑1 (MQ=255)
     gttgCGCGGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGAGGGCTTGcccc                                                                                                                                                  >  1:7298597/1‑149 (MQ=255)
        gcgcGGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGctc                                                                                                                                               <  2:6307550/150‑1 (MQ=255)
        gcgcGGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGctc                                                                                                                                               <  2:3839664/150‑1 (MQ=255)
        ccacgcgcatcatcgagctgtcggccgatggtgtggtggacttcagcggtacctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGctc                                                                                                                                               <  1:6812025‑M1/22‑1 (MQ=255)
         cgcgGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCAAGCCCGcgaa                                                                                                                                              <  1:5442625/150‑4 (MQ=255)
         cgcgGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCAAGCCCGcgaa                                                                                                                                              <  2:1164618/150‑4 (MQ=255)
         cgcgGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCAAGCCCGcgaa                                                                                                                                              <  1:779656/150‑4 (MQ=255)
           cgttCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTcc                                                                                                                                            <  2:4882216/147‑1 (MQ=255)
           cgGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCAAGCCCGcgaaga                                                                                                                                            <  1:7445608/150‑6 (MQ=255)
                aTGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCGCTCCCAAgg                                                                                                                                       >  1:8300011/1‑150 (MQ=255)
                 tGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCGCTACCAAGGc                                                                                                                                      >  1:568327/1‑150 (MQ=255)
                 tGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCGTGGCCGCGCTACCAAGGc                                                                                                                                      >  1:3404636/1‑150 (MQ=255)
                  gggCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGcc                                                                                                                                     <  2:6987944/150‑1 (MQ=255)
                   ggCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGccc                                                                                                                                    <  1:655502/150‑1 (MQ=255)
                      tctTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCcac                                                                                                                                 >  2:8523782/1‑150 (MQ=255)
                        tttCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCACGGCCCCACaa                                                                                                                               >  2:6812025/1‑150 (MQ=255)
                          tCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGc                                                                                                                             <  1:4537472/150‑1 (MQ=255)
                          tCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGc                                                                                                                             <  2:1853501/150‑1 (MQ=255)
                          tCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGc                                                                                                                             <  2:4287331/150‑1 (MQ=255)
                          tCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGc                                                                                                                             <  2:6453868/150‑1 (MQ=255)
                             cgggcgggcggTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTc                                                                                                                          <  1:90896/150‑1 (MQ=255)
                              gggcgggcggTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTcc                                                                                                                         <  1:4096308/150‑1 (MQ=255)
                                gcgggcggTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCaa                                                                                                                       <  1:7385774/150‑1 (MQ=255)
                                gcgggcggTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCaa                                                                                                                       <  1:100815/150‑1 (MQ=255)
                                   ggcggTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATAc                                                                                                                    <  1:2827345/150‑1 (MQ=255)
                                     cggTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCt                                                                                                                  <  2:7115304/150‑1 (MQ=255)
                                       gTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCtgt                                                                                                                <  1:1297368/150‑1 (MQ=255)
                                       gTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCtgt                                                                                                                <  1:436674/150‑1 (MQ=255)
                                        ttGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCtgtg                                                                                                               <  1:2741008/150‑1 (MQ=255)
                                         tGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGc                                                                                                              <  1:5425192/150‑1 (MQ=255)
                                               tGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTccc                                                                                                        <  2:3791510/150‑1 (MQ=255)
                                                  tCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCtgt                                                                                                     <  1:3978463/150‑1 (MQ=255)
                                                  tCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCtgt                                                                                                     <  1:2722553/150‑1 (MQ=255)
                                                   cagcggtacctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCtgtg                                                                                                    >  2:8645113‑M1/86‑150 (MQ=255)
                                                      atcatGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGa                                                                                                 <  2:5739264/150‑1 (MQ=255)
                                                         tacctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCg                                                                                              >  2:9089136‑M1/80‑150 (MQ=255)
                                                         tacctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCg                                                                                              >  1:7759379‑M1/80‑150 (MQ=255)
                                                          tGGCCGTTGCGCGAAGGAAAAGGCAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCgg                                                                                             <  1:993549/150‑1 (MQ=255)
                                                          tGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCgg                                                                                             <  1:8535572/150‑1 (MQ=255)
                                                          tGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCgg                                                                                             <  1:8535689/150‑1 (MQ=255)
                                                          tGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCgg                                                                                             <  1:2330111/150‑1 (MQ=255)
                                                          tGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCgg                                                                                             <  2:7799118/150‑1 (MQ=255)
                                                          tGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCgg                                                                                             <  2:1362584/150‑1 (MQ=255)
                                                          acctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCgg                                                                                             >  1:2553867‑M1/79‑150 (MQ=255)
                                                          acctacgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCgg                                                                                             >  1:4829806‑M1/79‑150 (MQ=255)
                                                           ggCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCggg                                                                                            <  1:7950914/150‑1 (MQ=255)
                                                           ggCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCggg                                                                                            <  1:2042120/150‑1 (MQ=255)
                                                            gCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGc                                                                                           >  1:5323019/1‑150 (MQ=255)
                                                              cGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCaa                                                                                         <  1:5003346/150‑1 (MQ=255)
                                                              cGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCaa                                                                                         <  2:2767785/150‑1 (MQ=255)
                                                              cGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCaa                                                                                         <  1:1862561/150‑1 (MQ=255)
                                                              acgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCaa                                                                                         >  2:474570‑M1/75‑150 (MQ=255)
                                                              acgatgactacctgcgtagccagggtgtgctggtctgattgatgggacggcctggccatctcggatggccgggcCTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCaa                                                                                         >  1:1455454‑M1/75‑150 (MQ=255)
                                                                ttGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGc                                                                                       <  1:3299429/150‑1 (MQ=255)
                                                                ttGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGc                                                                                       <  1:5997876/150‑1 (MQ=255)
                                                                            aaaGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCGTTCCCGCGCTACCAAGGCCCCACAAGCCTCCAAAACCAGTTCGGTGCCGGTGGGAGCGGGCACCCCCGCGCAGAGGc                                                                           >  1:7067673/1‑150 (MQ=255)
                                                                                  agcagGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCCGcc                                                                     >  2:5540230/1‑150 (MQ=255)
                                                                                       gCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCGTGCCCCCGCGCCCCACGCCCCCCAACCCGCCAAAACCACTGCGGTCCCGGGGGGCGCCGGCCCCCCCCCGAACACGCCCCGccccccc                                                                >  1:5606994/1‑145 (MQ=255)
                                                                                         tAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGTTGGGGATGTACTCTTCGGGGGCTTGCCCGCTCGCCCACGGCCCCACACGCGTCCAATACCTGTGCGGTCCCTGTGGGAGGGGGCAAGCCCGCGAAGAGGCCCGGCCATCCgag                                                              <  2:2557084/150‑1 (MQ=39)
                                                                                          aaCGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACCAGCGTCCAAAACCGGTGCGGTCCCGGTGGGAGCGGGCACGCCCGCGACGAGGCCCGGCCAGCCgcga                                                             >  1:8211160/1‑150 (MQ=255)
                                                                                          aaCGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCTGGCCATctcgga                                                             <  2:1861825/150‑6 (MQ=255)
                                                                                             gAATTATTCGTCGGGCCGTGCCGGTCGGTGGTATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCAAATACCTGTGCGGTCCCTGTGGTAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATgg                                                          <  2:222097/150‑1 (MQ=255)
                                                                                               attattCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCGGTGGGAGCGGGCCAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGcc                                                        >  2:6921699/1‑150 (MQ=255)
                                                                                                 tattCGTCGGGCGGTGCCGGTGGGTGGGTTGGGGATGTACTCTTCGCGGGCTGGCCCGCTCACCCCAGGCCCCACAAGCGTCCAATCCCTGTGCGATCCCTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAg                                                      <  1:3698103/150‑1 (MQ=255)
                                                                                                   ttCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCGGTGGGAGCGGGCAACCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGc                                                    >  2:4344710/1‑150 (MQ=255)
                                                                                                      gtcgGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATCCCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGt                                                 <  2:7934097/150‑1 (MQ=255)
                                                                                                      gtcgGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCGCTCCCCAGGCCCCACCAGCGTCCAATAACAGTGCGGGCCCGGTGGGAGCGGGCAAGCCCGCGCAGAGGCCCCGCCATCCCAGATGGCCccgccgg                                                 >  1:576209/1‑149 (MQ=255)
                                                                                                       tcgGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTGACCCAAGGCCCCACCAGCGTCCAATACCAGTGCGGGCCCGGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCCGCCATCCCAGAAGGGCCGGCCGTc                                                >  1:5268873/1‑150 (MQ=255)
                                                                                                          ggCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCa                                             <  2:5835850/150‑1 (MQ=255)
                                                                                                          ggCCGTGCCGGCCGGTGGTATGGGGATGTACTCTTCGCGGGCTTGCCCGCACTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCa                                             <  2:1806027/150‑1 (MQ=255)
                                                                                                           gCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCat                                            <  2:266308/150‑1 (MQ=255)
                                                                                                               tGCGGGTGGGTGGGTTGGGGATGTACTCTGCGCGGGCTTCCCCGCTCCCCCACGCCCACACACGCGTCCAATACGTGTGGGGGCCGGGGGGGCGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCatcaat                                        <  1:6554689/150‑1 (MQ=255)
                                                                                                                gCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCGCTCCCAAGGCCCCCCAAGCGTCCAAAAACAGTGCGGTCCCAGTGGGCGCGGGCCACCCCGCGAAGTGGCCCGGCCAACCCAGATTGCCCCGCCGTCCCctcaaaa                                       >  2:1564173/1‑148 (MQ=255)
                                                                                                                 ccGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTGGCCCGCTCACCCAAGGCCCCACAAGCGTCCAATCCCTGTGCGGGCCCTGGGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCatcaatca                                      <  1:8941944/150‑1 (MQ=255)
                                                                                                                 ccGGTCGGTGGGATGGGGATGTACTCGTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATCCCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCatcaatca                                      <  2:6003744/150‑1 (MQ=255)
                                                                                                                    gtcggtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGAc                                   >  2:5197471/1‑150 (MQ=255)
                                                                                                                       ggtggtATGGGGATGTAATCTTCGCGGGCTTGCCCGCTCTGCCAAGGCCCCACAAGCGTCCAATACCTGTGGGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGGCATCCGAGATGGCCAGGGCGTGCCATCAATCAGACCAg                                >  2:2171494/1‑150 (MQ=38)
                                                                                                                       ggtGGGATGGGGATCTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACCAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCCTCCCTCAGACCAg                                >  2:5217112/1‑150 (MQ=39)
                                                                                                                        gtGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGATGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGc                               <  2:4536216/150‑1 (MQ=255)
                                                                                                                           ggatggggatgTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGcaca                            >  2:8906415/1‑150 (MQ=255)
                                                                                                                               ggggatgTACTCTTCGGGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCGGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCt                        <  1:2055331/150‑1 (MQ=255)
                                                                                                                                 ggatgTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTgg                      <  2:7407159/150‑1 (MQ=255)
                                                                                                                                         ggtgcgcgGGCGGGCCCGCCCCCCCACGCCCCCACACGCGACCAACACCTGGGCGGTCCCGGGGGGCGCGGGCACCCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAg              <  1:8523782/146‑1 (MQ=255)
                                                                                                                                         aCTTCGCGGGCTTGCCCGCTCTCCCACGGCCCCACTTGCGTCCAATACCTGTGGGGTCCCTGTGGTAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGTCTACGcat              >  2:6574696‑M1/2‑149 (MQ=255)
                                                                                                                                         TCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAg              >  2:4574515‑M1/1‑150 (MQ=255)
                                                                                                                                         TCTTCGCGGGCTTGCCCCCTCTCCCAAGGCCCCACAATCGTCCAATACCTGTGCGGTCCCGGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGCTGGCCCCGCCGTCCCATCAAACAGACCAGCACACCCAGGCTAcgccg              <  2:5632294‑M1/150‑3 (MQ=255)
                                                                                                                                           TTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGt            <  2:7551105‑M1/150‑1 (MQ=255)
                                                                                                                                            TCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTa           >  2:4449680‑M1/1‑150 (MQ=255)
                                                                                                                                                  GCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATc     <  1:5013014‑M1/150‑1 (MQ=255)
                                                                                                                                                     TGCCCGCTCTCCCACGGCCCCACAAGCGTCCACTACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTa  >  2:6928656‑M1/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTTGTTGCGCGGTCATGGGCTCTTTCGCGGGCGGGCGGTTGGGCGTGGTCCCATCATGGCCGTTGCGCGAAGGAAAAGACAGCAGGCTAACGAATTATTCGTCGGGCCGTGCCGGTCGGTGGGATGGGGATGTACTCTTCGCGGGCTTGCCCGCTCTCCCAAGGCCCCACAAGCGTCCAATACCTGTGCGGTCCCTGTGGGAGCGGGCAAGCCCGCGAAGAGGCCCGGCCATCCGAGATGGCCAGGCCGTCCCATCAATCAGACCAGCACACCCTGGCTACGCAGGTAGTCATCGTA  >  NC_002947/3340946‑3341244

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: