Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,491,722 T→A 7.4% intergenic (‑102/+3) yhfG ← / ← ppiA putative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,491,7220TA 7.4% 1336.6 / 12.4 440intergenic (‑102/+3)yhfG/ppiAputative antitoxin for Fic/peptidyl‑prolyl cis‑trans isomerase A (rotamase A)
Reads supporting (aligned to +/- strand):  ref base T (142/247);  new base A (45/2);  total (191/249)
Fisher's exact test for biased strand distribution p-value = 4.19e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 441 total aligned reads displayed.

TGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGAAATC  >  NC_000913/3491574‑3491871
                                                                                                                                                    |                                                                                                                                                      
tGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATc                                                                                                                                                       >  2:4701625/1‑150 (MQ=255)
 gAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCa                                                                                                                                                      <  2:2189517/150‑1 (MQ=255)
 caaccTCCCACACACGGGACATTTGCTTATCGCTGACTTTCTTCACTTGATGCCTCCCAGAACCTCTCCAATTCCCACAACTATAAGCCGCAAATCCGCGTTACGCCGCGAGAGCAGAAGCGCGGGCAAGATTCCCCGCGCGAGAAAACa                                                                                                                                                      <  2:8164019/146‑1 (MQ=255)
  aaGCTCCCAGAGACGGGACTTTTGCTTATCGTTGTGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCAAAAAGTATAAGCGGAAATTCCGGGTTACGCCGGGAGAGCATAAGCGGGGGCAAGATTGCCCGCGCGAGAAATCAt                                                                                                                                                     <  1:6839368/150‑1 (MQ=255)
      tCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACg                                                                                                                                                 <  2:3994789/150‑1 (MQ=255)
          agagACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGATTGCCCGCGCGAGAAAGCATTAAGGCAg                                                                                                                                             >  1:625423/1‑150 (MQ=255)
                gggACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGGGACATTTCCCGCGCGGGCAAACAATAAGGCCGGACatt                                                                                                                                       >  2:7113950/1‑150 (MQ=255)
                gggACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGTTTTCCCGCGCGCGGAAAACTTAAGGCAGGACgtt                                                                                                                                       >  1:7092341/1‑150 (MQ=255)
                gggACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGATTGCCCGCGCGCGAAATAAGTACCGCAGGACttt                                                                                                                                       >  1:436694/1‑150 (MQ=255)
                gggACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCAATACGGCAGGACttt                                                                                                                                       >  2:3875678/1‑150 (MQ=255)
                 ggACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTa                                                                                                                                      <  1:7026584/150‑1 (MQ=255)
                   aCTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCGGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGc                                                                                                                                    <  1:7927094/150‑1 (MQ=255)
                    cTTTTGCTTTTCGTTGAGTTTCTTCAGTTGTTGCCTCCCTTAACGTGTCATTTTGCCACAATTATAAGCGGCAAATGCGGGTTTCGCGGGGAGATCATAAGCGCGGGCAATATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCg                                                                                                                                   <  1:5522004/150‑1 (MQ=255)
                     ttttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACtttttcgg                                                                                                                                  >  1:1506099/1‑145 (MQ=255)
                     ttttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAACAATGCCCCCGCGCGAAAACAATAAGGCACGACATTAACgg                                                                                                                                  >  1:3743292/1‑150 (MQ=255)
                       ttGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGCTTGCCCGCGCGAGAAAACAATACCGCCGGACTTTAGCGGaa                                                                                                                                >  1:6068521/1‑150 (MQ=255)
                             aTCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATa                                                                                                                          <  1:6523076/150‑1 (MQ=255)
                              tCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATaa                                                                                                                         <  2:3246583/150‑1 (MQ=255)
                                  tGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTAc                                                                                                                    <  1:7685128/150‑1 (MQ=255)
                                        tcttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATAAATAAGGCAGGACTTTAGCGGAAAGGATAAATACCGGttt                                                                                                              >  2:7722326/1‑150 (MQ=255)
                                          ttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGATTGGCCGCGCGAGGAAACAATACGGCAGGACTTTATCGGAAAGGATAACAAACCGGtttt                                                                                                             >  1:7220879/1‑150 (MQ=255)
                                          ttcACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTATCGGAAAGGATAAATACCGGttttt                                                                                                            >  2:9229628/1‑149 (MQ=255)
                                            cACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGAc                                                                                                          <  2:6446329/150‑1 (MQ=255)
                                                 gATGCCTCCCTGAACGTGTCCATTTGCCACACGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGcac                                                                                                     <  1:4649910/150‑1 (MQ=255)
                                                     ccTCCCTGGACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCCCGCGTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGGTAACTACCGGTTTTGACGGCACATTc                                                                                                 >  2:9318955/1‑150 (MQ=255)
                                                      cTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCCACATTTCCCGCGCCCGGAAACAATACGGCAGGCCTTTTGCGGGAACGATAAAGACACGTTTTGTCGGCCCAATCt                                                                                                >  2:341681/1‑150 (MQ=255)
                                                      cTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCCGATTGCCCGCGCCGGAAAACATTACGGCAGGCCGTTAGCGGAAAGGATAACTACCGGGTTTGACGGCCCAttag                                                                                                >  2:5037185/1‑148 (MQ=255)
                                                      cTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCt                                                                                                >  2:6006304/1‑150 (MQ=255)
                                                       tCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCCCCGTGCCCCCGCGCCCAAACAATACGGCAGGCCGGTAGCGGGAAGGCTAACAACCCCGGTTGGCGGcacagactg                                                                                               >  1:1124533/1‑145 (MQ=255)
                                                       tCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTg                                                                                               <  2:5681415/150‑1 (MQ=255)
                                                         ccTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGGTACGCCGGGAGAGCAGAAGCGCGGGCAACATTGCCCGCGCGAGCAATCAATACGGCCGGACGTTAGCGGCAACGCTAACTACCCGTTTTGACGGCACCGTCTggg                                                                                             >  2:8481606/1‑149 (MQ=255)
                                                           tGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTAc                                                                                           <  1:4302025/150‑1 (MQ=255)
                                                             aaCGTGTCCATTCGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGTGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACgg                                                                                         <  1:8181400/150‑1 (MQ=255)
                                                               cGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGAc                                                                                       <  2:301350/150‑1 (MQ=255)
                                                                  gtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCaa                                                                                    <  2:7088154/150‑1 (MQ=255)
                                                                  gtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCaa                                                                                    <  2:1406222/150‑1 (MQ=255)
                                                                  gtCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAACATTGCCCGCGCGAGAAATAATTAAGGCAGGACTTTAGCGGAAAGGATAACTAACCGTTTTTAAGGCACATTATGGTAAGGACCaa                                                                                    >  2:5168214/1‑150 (MQ=255)
                                                                     cATTTGCCACAAGTATAAGCGGCGAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGt                                                                                 <  2:1806129/150‑1 (MQ=255)
                                                                       tttGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCa                                                                               <  1:2398916/150‑1 (MQ=255)
                                                                        ttGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCAt                                                                              <  2:620423/150‑1 (MQ=255)
                                                                          gCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGGGAGATTTGGCGCGCGGGAAATAATTAAGTCAGGAGTTTTTGGGTAAGGATAAGTAACGGGTTTTAGGGCACATTATTGTTCGGTCCAACGTCATGa                                                                            >  1:3124696/1‑150 (MQ=255)
                                                                            cacaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGt                                                                          <  2:1032150/150‑1 (MQ=255)
                                                                            cacaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGt                                                                          <  1:4022107/150‑1 (MQ=255)
                                                                             acaAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCCAGATTGCCCCCCCGAGAAATAAATAAGGCCGGACGATATCGGAAAGGCGAAAAACCCGGTTTTAAGGCCCAATATGGTACGGACCAACGTCAAGAGTc                                                                         >  2:2560184/1‑150 (MQ=255)
                                                                                    taAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGTCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCt                                                                  >  2:4278600/1‑150 (MQ=255)
                                                                                    taAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCt                                                                  <  2:1492659/150‑1 (MQ=255)
                                                                                    taAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCt                                                                  <  2:6942440/150‑1 (MQ=255)
                                                                                     aaGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTg                                                                 <  2:4502362/150‑1 (MQ=255)
                                                                                     aaGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTg                                                                 <  2:5852026/150‑1 (MQ=255)
                                                                                      aGCGGCAAATCCGGGTTACGCCGTGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGTAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTgg                                                                <  2:2255265/150‑1 (MQ=255)
                                                                                       gCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCCAGATTTCCCCCGCGAGAAAAAATTACGGCAGGACTTTATCGGAAAGGATAACAACCCGTTTTTACGGCACATAATGGTAAGGGCCAACCTCATGAGTCGGCACCTggg                                                               >  2:3511670/1‑150 (MQ=255)
                                                                                         ggCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGGGGGCAAGAGTGGCCCCGCGAGAAAAAATTAAGGCAGGAGTTTATCGGAAAGGAAAACTACCCGTTTTTTAGGCAAATACTTGTACGGAACAACCTAATGAGTCGGCAACAGGGaa                                                             >  2:8081040/1‑150 (MQ=255)
                                                                                            aaaTCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATc                                                          <  1:8836253/150‑1 (MQ=255)
                                                                                             aaTCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCt                                                         <  1:7377986/150‑1 (MQ=255)
                                                                                                ccGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGACATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTAt                                                      <  2:5433715/150‑1 (MQ=255)
                                                                                                 cGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATc                                                     <  2:4668901/150‑1 (MQ=255)
                                                                                                      tACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCaa                                                <  1:626263/150‑1 (MQ=255)
                                                                                                       aCGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAAc                                               <  2:7374477/150‑1 (MQ=255)
                                                                                                         gCCGGGAGAGCAGAAGCGCGGGGCAGATTGCCCGCGCGCGAAAAAAATTCGGCAGGACGTTAGCGGAAAGGATAACGAACCGGGTTTTCCGCACAATCGGGGACGGCCCAACCGAATGAGGCGGCACCGGGGCCATCGTATCGGCCACgg                                             >  2:8740229/1‑149 (MQ=255)
                                                                                                          ccGGGAGAGCAGAAGCGCGGGCCACCATGCCCGCCCCCGAAAAAAATACAGCAGGACGTTTTCGGAAAGGATAAAAAAACGTTTTTACCGCACAATCTGGTCCGGACCACACTCAAGAGTAGGCGCCGGGGCAAACGTATAGGCAAccgc                                            >  2:2711992/1‑147 (MQ=255)
                                                                                                             ggAGAGCAGAAGCGCGGGCCAGAATGCCCGCGCCAGAAATAAATAAGGCAGGACATTAGCGGAAACGATAACTACCCGTGTTGACGGCACAATCTGGTACGGACCCACGTCAAGAGTCGGCCCCGGGGGAAACTTATCGGCAACGGCCAt                                         >  2:2899511/1‑150 (MQ=255)
                                                                                                             ggAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGGCCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCAt                                         <  2:1724571/150‑1 (MQ=255)
                                                                                                             ggAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCAt                                         <  2:9031293/150‑1 (MQ=255)
                                                                                                             ggAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCAt                                         <  2:4971858/150‑1 (MQ=255)
                                                                                                             ggAGAGCAGAAGCGCGGGCAAGATGGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCAt                                         <  1:830082/150‑1 (MQ=255)
                                                                                                                gagCATAAGCGCGGGCAAGTTTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGcc                                      <  1:2548298/150‑1 (MQ=255)
                                                                                                                 agCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCt                                     <  2:2050983/150‑1 (MQ=255)
                                                                                                                  gCAGACGCGCGGTCAAGATTGCCCGCGCTAGAAATCATTACGGCAGGACTTTAGCGTAATTTATAACTACCGTTTTTGACGGCACATTCTGGTACGGACAAACGTCATTAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCtt                                    <  1:3312941/150‑1 (MQ=255)
                                                                                                                  gCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCtt                                    <  1:7962224/150‑1 (MQ=255)
                                                                                                                  gCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCtt                                    <  1:535844/150‑1 (MQ=255)
                                                                                                                   cAGAAGCGCGGGCAAGATTTCCCCCGCGAGAAATAAATAAGGCAGGACTTTAGCCGAAAGGATAAAAAACGGTTTTTACGGCACAATATGGTACGGACCAACGTCAAGAGTCGGCACCTGGGAAATATTATAGGCAACGTCCATGCCttt                                   >  2:7636082/1‑150 (MQ=255)
                                                                                                                   cAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCttt                                   <  2:9161351/150‑1 (MQ=255)
                                                                                                                     gAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTca                                 <  2:6669665/150‑1 (MQ=255)
                                                                                                                     gAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACAACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTca                                 >  2:7932643/1‑150 (MQ=255)
                                                                                                                        gcgcGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTcacca                              >  2:6223868/1‑150 (MQ=255)
                                                                                                                             ggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTa                         >  2:1989524/1‑150 (MQ=255)
                                                                                                                             ggCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTa                         <  2:4817379/150‑1 (MQ=255)
                                                                                                                              gCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAc                        >  1:5786604/1‑150 (MQ=255)
                                                                                                                               cAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTAcc                       <  2:6571633/150‑1 (MQ=255)
                                                                                                                                aaGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCa                      <  2:6772152/150‑1 (MQ=255)
                                                                                                                                  gATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCaaa                    >  2:1729599/1‑150 (MQ=255)
                                                                                                                                  gATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTATCACCACTTTACCaaa                    >  1:6180820/1‑150 (MQ=255)
                                                                                                                                   aTTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAAt                   >  2:117393/1‑150 (MQ=255)
                                                                                                                                   aTTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAAt                   >  1:4823689/1‑150 (MQ=255)
                                                                                                                                     tGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATAc                 >  2:3781408/1‑150 (MQ=255)
                                                                                                                                       cccGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcg               >  2:7495196/1‑150 (MQ=255)
                                                                                                                                       cccGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcg               <  1:5116736/150‑1 (MQ=255)
                                                                                                                                       cccGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcg               >  1:8609263/1‑150 (MQ=255)
                                                                                                                                        ccGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcgc              <  1:2752748/150‑1 (MQ=255)
                                                                                                                                        ccGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcgc              <  1:5090717/150‑1 (MQ=255)
                                                                                                                                        ccGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACcgc              <  2:5596681/150‑1 (MQ=255)
                                                                                                                                            gcgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTaa          >  1:1716854/1‑150 (MQ=255)
                                                                                                                                            gcgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTaa          >  2:1304795/1‑150 (MQ=255)
                                                                                                                                            gcgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTaa          <  2:5315019/150‑1 (MQ=255)
                                                                                                                                             cgAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAAc         <  2:6177481/150‑1 (MQ=255)
                                                                                                                                               agaAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGTTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCg       >  1:9207495/1‑150 (MQ=255)
                                                                                                                                                gaAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGa      >  1:5021841/1‑150 (MQ=255)
                                                                                                                                                gaAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGa      >  1:6460093/1‑150 (MQ=255)
                                                                                                                                                 aaaTCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGaa     <  1:4478111/150‑1 (MQ=255)
                                                                                                                                                   aTCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGAAAt   >  2:835033/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TGAAGCTCCCAGAGACGGGACTTTTGCTTATCGGTGAGTTTCTTCACTTGATGCCTCCCTGAACGTGTCCATTTGCCACAAGTATAAGCGGCAAATCCGGGTTACGCCGGGAGAGCAGAAGCGCGGGCAAGATTGCCCGCGCGAGAAATCATTACGGCAGGACTTTAGCGGAAAGGATAACTACCGGTTTTGACGGCACATTCTGGTACGGACCAACGTCATGAGTCGGCACCTGGGAAATCTTATCGGCAACGTCCATGCCTTTCACCACTTTACCAAATACCGCGTAACCGAAATC  >  NC_000913/3491574‑3491871

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: