Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913501,5570AC32.8% 144.0 / 10.4 61intergenic (+128/+5)gsk/ybaLinosine/guanosine kinase/putative transporter YbaL
Reads supporting (aligned to +/- strand):  ref base A (37/4);  new base C (0/20);  total (37/24)
Fisher's exact test for biased strand distribution p-value = 1.70e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGC  >  NC_000913/501309‑501795
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
gcAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAaccacc                                                                                                                                                                                                                                               <  1:457023/251‑1 (MQ=255)
  aGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGt                                                                                                                                                                                                                                            <  1:723910/251‑1 (MQ=255)
          aaaTATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTAcccc                                                                                                                                                                                                                                     <  1:79187/250‑1 (MQ=255)
           aaTATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGt                                                                                                                                                                                                                                   <  1:408888/251‑1 (MQ=255)
                gCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcacc                                                                                                                                                                                                                               >  1:303189/1‑250 (MQ=255)
                    aCCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTc                                                                                                                                                                                                                          >  2:181920/1‑251 (MQ=255)
                           gAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAGGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccgcc                                                                                                                                                                                                                   <  1:193698/251‑1 (MQ=255)
                                aTCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATAGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTgg                                                                                                                                                                                                              >  2:86352/1‑251 (MQ=255)
                                                   ttccccTCGTTAACCGCGCGGTTTGCGGGGGCGTAAAGAGAGCGTGAAAGGGTCTAACGGGGATCTTTAATTTATGGTGGTTCGGAGGGCCAGAAAGGCCTTTCCCGCCGCTTCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCTTCTTTTCTTGCCTACTTTTTTTTCATATTTGACACCCGGCACCCCCCGTTTCCCCCGTCACCCCCACACCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                                                                                                                            <  2:571678/246‑1 (MQ=255)
                                                     ccccTCGTTTACCGCGCGGTTTCCCGGGGGGTGAAGACGGCCTGGAGGTGTCTTACTGGGTTCTTTAAGTTATCGTCGGTTCGTAGGCCAGATAGGGCTTTCCCGCCGCTTCGGGCATTTGGCTCTGGATGCCTGATGCGCCGCTGGCGCTTCTTACCTTGCCTACATTTTTTTCATATTTTACACCCGGAACCCCCCGTTTACCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCAt                                                                                                                                                                                         <  1:740707/249‑1 (MQ=255)
                                                              tAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGc                                                                                                                                                                                >  2:211631/1‑251 (MQ=255)
                                                                        cTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                                                                                                      >  1:661893/1‑251 (MQ=255)
                                                                             cGGAGCGTGAATACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATAGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCCACCACCGTTTACCCCGTCACAACCTCACCCGCCGGTGGCGTTTCCCGCAGTTCCAGCATCGTACGGGCGATTTCAcgcccgcc                                                                                                                                                                 >  2:458907/1‑251 (MQ=255)
                                                                             cGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGcc                                                                                                                                                                 >  1:350829/1‑251 (MQ=255)
                                                                                   gTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCAc                                                                                                                                                           >  2:431725/1‑251 (MQ=255)
                                                                                    agaagaccgctggaaagatttttACGGGGATGGTAAATTTATCGGCGTTTCGGAGGCCAGAAAGGGCGTTCCCGCCGCATCTGGCTTTGGCCTCTCTATCCCTGATGCGACGCGGGCGCGTCTTATCATGCCTACATTTTTTTCATTTTTTACACCCGGCACCCCCCGTTTCCCCCGTCACCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                                                                                                          <  2:189730/230‑1 (MQ=255)
                                                                                        gACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCt                                                                                                                                                      >  2:168472/1‑251 (MQ=255)
                                                                                            agcaggAAGAGTTTTACGGGGATCGTTAAGTTATCGTCGTTTCGTAGGCCAGAAAAGGCGTTCCCGCCGCATCTGGTTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCAACCCCCGTTTACCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                                                                                                                  <  1:43024/247‑1 (MQ=255)
                                                                                                            gtGGGATGGTTAATTTATCGTGGGTTGGTGGGCAAGATAAGGCGTTCCCGCCGCATCTGGCTTTGGGCTCTCGTGGCTTGATGCGACGCGGGCGGTTTTTATCATGCCTAAATATTTTTCATATTTTCCACCCGGCAACCCCCTTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                                                                                                                   <  1:755151/249‑1 (MQ=255)
                                                                                                              gggATCGTTACGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                                                                                                                >  2:598062/1‑251 (MQ=255)
                                                                                                              gggATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGat                                                                                                                                >  1:231691/1‑251 (MQ=255)
                                                                                                                    tttAAGTTAGCTTCTGTTCGTAGGCCAGATAAGGCTTCCCCCCCGCTTCTGGCTTTTGTCTCTCGATGCCTGATGCGACGCTGGCGCGTTTTATCAGCCCTACATTTTTTTCATTTTTTACACCCGCCAACCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                                                                                                                           <  1:257713/249‑1 (MQ=255)
                                                                                                                              cgcgggTTCGTGGGCCAGAAAGGGCGTTCCCGCCGCATCTGGATTTTGCCTCCCGATGCCTGATGCGCCGCGGGCGCGTTTTATCATCCCTACATTTTTTTCATTTTTTACATCCGGCAACCCCCGTTTACCCCGTCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCg                                                                                                                <  1:698925/247‑1 (MQ=255)
                                                                                                                               gtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGt                                                                                                               >  2:21576/1‑251 (MQ=255)
                                                                                                                               gtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGACATACGCCACTTCATCGt                                                                                                               >  2:149016/1‑251 (MQ=255)
                                                                                                                                tcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc                                                                                                              >  1:614850/1‑251 (MQ=255)
                                                                                                                                  ggTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt                                                                                                            >  2:13776/1‑251 (MQ=255)
                                                                                                                                    ttCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATa                                                                                                           >  1:76887/1‑250 (MQ=255)
                                                                                                                                      cGTAGCCCAGATAGGGCGTTCCCGCGGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCTTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCCCCGTTTCCCCCGCCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAAt                                                                                                         <  1:13776/250‑1 (MQ=255)
                                                                                                                                         cggCCAGATAAGGCTTTCCCGCCCCACCTCGCTTTTGGCTCCTGAGGCCTGATGCGCCGCTGGCGCGTCTTATCTCCCCTAATTTTTTTCCATATTTGACACCCGGCATCCCCCGTTTACCCCGCCACCACCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgc                                                                                                     <  1:680791/250‑1 (MQ=255)
                                                                                                                                                          ttccccccgcATCTGGCATTTGCCTCTCGTTGCCTGAGGCGACGCGGGCGCGTCTTATCAGGCCTACATTTTTTTCATATTTAACACCCGCAACCCCCCGTTTCCCCCGCCACCACCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc                                                                                        <  1:598062/243‑1 (MQ=255)
                                                                                                                                                               gccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTcc                                                                                                                                                                                                       >  1:661327/1‑133 (MQ=255)
                                                                                                                                                                ccccATCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCCCCTTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATc                                                                                <  1:571985/248‑1 (MQ=255)
                                                                                                                                                                      cTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGAtt                                                                          >  1:392301/1‑251 (MQ=255)
                                                                                                                                                                        ggCATTTGGCTCTCGATGCCTGAGGCGCCGCTGGCGCGTTTTATCATCCCTACTTTTTTTTCATATTTTACATCCGGCACCCCCCTTTCCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGAtttt                                                                        <  1:289091/251‑1 (MQ=255)
                                                                                                                                                                         gCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGAttttt                                                                       >  2:442453/1‑251 (MQ=255)
                                                                                                                                                                           aTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATgg                                                                                                       >  1:294539/1‑217 (MQ=255)
                                                                                                                                                                            tttGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTcgc                                                                    >  1:789181/1‑251 (MQ=255)
                                                                                                                                                                               ggCTCTCGATGCCTGATGCGGCGCTGGCGCGTTTTATCTTGCCTACATTTTTTTCATATTTTCCACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTcgcgc                                                                  <  1:211631/250‑1 (MQ=255)
                                                                                                                                                                                 ctctCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCggg                                                               >  1:69691/1‑251 (MQ=255)
                                                                                                                                                                                   nctcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGc                                                              >  2:344356/2‑251 (MQ=255)
                                                                                                                                                                                       gATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGc                                                              >  1:636176/1‑247 (MQ=255)
                                                                                                                                                                                          cccTGATGCGACGCTGGCGCTTGTTATCATGCCTACTTTTTTTTAATTTTTCCCGCCCGCCACCCCCCGTTTCCCCCGCCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGcn                                                       <  2:350829/250‑2 (MQ=255)
                                                                                                                                                                                              gATGCGCCGCTGGCGCCTCTTCTCTTGCCTACTTTTTTTACATATTTTACATCCGCCAACCCCCGTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCcac                                                     <  2:76887/249‑1 (MQ=255)
                                                                                                                                                                                               atgccccGCTGGCGTGTTTTATCATCCCTACTTTTTTTTCATATTTTCCCCCCGCCACCCCCCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAAt                                                  <  2:47537/245‑1 (MQ=255)
                                                                                                                                                                                                 gCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAAtct                                                >  1:433514/1‑251 (MQ=255)
                                                                                                                                                                                                   gACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAAtctc                                               >  2:703752/1‑250 (MQ=255)
                                                                                                                                                                                                     cGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCAcc                                            >  2:37465/1‑251 (MQ=255)
                                                                                                                                                                                                      gCTGGCGCGTCTTACCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCAccc                                           >  1:488314/1‑251 (MQ=255)
                                                                                                                                                                                                       cTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCg                                          >  2:787114/1‑251 (MQ=255)
                                                                                                                                                                                                          gcgcGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAAtc                                                 >  2:735941/1‑241 (MQ=255)
                                                                                                                                                                                                            gcttcttctCATGCCTATATATTTTTAATATTTACCACCCGACACCCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCa                                     <  2:69691/248‑1 (MQ=255)
                                                                                                                                                                                                              gTCTTATCATGACTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGAGGAGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATACTCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCATCCGCTTCAt                                    >  1:148550/1‑250 (MQ=255)
                                                                                                                                                                                                                   aTCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGt                              >  1:308338/1‑251 (MQ=255)
                                                                                                                                                                                                                         ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGa                         >  1:401330/1‑250 (MQ=255)
                                                                                                                                                                                                                             catttttttCATATTTTACACCCGGCAACCCCCGTTTACCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGt                    <  1:168472/247‑1 (MQ=255)
                                                                                                                                                                                                                                atTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTc                   >  1:182773/1‑249 (MQ=255)
                                                                                                                                                                                                                                   tttCATTTTTTCCACCCGGCACCCCCCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGa               <  1:431725/250‑1 (MQ=255)
                                                                                                                                                                                                                                       catattttccccccGGCACCCCCCGTTTCCCCCGCCCCCCCCCCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCa            <  1:37465/238‑1 (MQ=255)
                                                                                                                                                                                                                                         tatTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGcc         >  1:161182/1‑251 (MQ=255)
                                                                                                                                                                                                                                          atTTTACATCCGTCAACCACCGTTTACTACGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGCCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGCATCGTCCGGATCAGcc         >  2:397505/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 aTCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGc  >  1:700535/1‑250 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
GCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGC  >  NC_000913/501309‑501795

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: