Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,1320GA27.8% 98.8 / 14.3 55A16V (GCT→GTT) rarDputative transporter RarD
Reads supporting (aligned to +/- strand):  ref base G (8/31);  new base A (15/0);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 5.67e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCA  >  NC_000913/4003893‑4004372
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
tGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGc                                                                                                                                                                                                                                           <  1:235131/250‑1 (MQ=255)
    cACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                                                                                                                                                                                                                       <  2:113793/250‑1 (MQ=255)
     aCCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GGAGTAAAGCAACTTTAAGTACGCGCGCGCCATTCCCCCCATAAAATAACAGGCAAGAGAGAGTa                                                                                                                                                                                                                                      >  1:230578/1‑250 (MQ=255)
        ccccAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                                                                                                                                                                                                                                  <  2:43980/250‑1 (MQ=255)
             gATAAACAGTAGCCAGTTGCCACCAATCAGAACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCGGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACGTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGcccg                                                                                                                                                                                                                             <  2:345149/251‑2 (MQ=255)
                aaaCAGTAGCCAGTTGCCACCACTCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCCGGCGCGTTTCCCCCAATAAAATAAAAGGCCAAAGCGGGTAATAACCCCTGcc                                                                                                                                                                                                                          >  1:344408/1‑251 (MQ=255)
                 aaCAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcg                                                                                                                                                                                                                         <  1:176179/251‑1 (MQ=255)
                   cAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATTACGCGATGCGTCAGGATGTCAGCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATTCCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgc‑g                                                                                                                                                                                                                      <  1:58882/251‑1 (MQ=255)
                      tAGCCAGTTGCCCCCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑gtt                                                                                                                                                                                                                    <  1:205844/250‑1 (MQ=255)
                             ttGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GGAGTAAAACAACATTAAGTAAGCCTGCGCTATTCCCCAAAAAAAATAAACCGCAAAAACGAGGAATAAGCCCTGCCGC‑GGTTGtttttc                                                                                                                                                                                                            >  1:357668/1‑249 (MQ=255)
                                   ccaATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGCCTGACTCAGCGTTTTTAAATAGGACCACTGGCGGCCCATGCTCCTCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGG‑ACAGTAGTAAATCACCTTGGAGCCAGCTGGCGCTATACCCCACAAAAAATAAACGGCCAGGACGAGAAATACTCCCCTCCGC‑GTTTGtttttcctccct                                                                                                                                                                                                      >  1:278610/1‑243 (MQ=255)
                                            acgtcAGAGACTAGCAGCGTAAAACGTTTCTGTGGCGTCCGAATCAGCGTTTGTAAATAGGGCCACTGGCGGCAAATGCTCATCAGCACCACCAGAAAGAAAAACGACCAGAGCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCAGCAAGGGCGAGTAATGCGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATAc                                                                                                                                                                                             <  2:260/247‑1 (MQ=255)
                                            aCGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAAGAGGACCACTGGCGGCAAATGCTCAGCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATGCCCCAAATAAAATAAGCTGCAAGAGGTAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATAc                                                                                                                                                                                             <  1:29458/251‑1 (MQ=255)
                                                cAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCAt                                                                                                                                                                                         <  2:344408/251‑1 (MQ=255)
                                                           aGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCC‑GTAGTAAAACAAATTGAAGTACGCTGGCGCTCTACCCCCAATAAAATAAACGGCACGAGCGGGGAATAAGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGGATACTCATTTTTAATTGGt                                                                                                                                                                              >  1:114275/1‑251 (MQ=255)
                                                             cATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGAC‑GTAGGAAATCAACTTGAAGGACGCCGGCGCTATAACCCCAATAAAAAAAAAGGCAAGAGCCAGTAATACGCCCTCCCGC‑GGTTGTTTTTGATCCAACGGGAATAATCCTTTTTATTTGGtg                                                                                                                                                                             >  1:76902/1‑249 (MQ=255)
                                                             cATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTACATTAACTTGAAGAACGCTGGCGCTCTACCCCAAAAAAAAAAAACGGCCAGAGCGGGGAATACGCCCCTCCGC‑GGTTGTTTTGCATCCCTCGGGCAAACTCATTTTTAATTGGTa                                                                                                                                                                             >  1:14673/1‑250 (MQ=255)
                                                                    aTTTTCTGTGGAGGCTGAAGCTGCGTTTTTAAATAGGACCACTGACGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCATGATTTCATCGGTGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑Gt                                                                                                                                                                    <  2:10993/251‑1 (MQ=255)
                                                                          tgtgGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCCAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGCTGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAAATTGAAGTACGCTGGCGCTATACCCCAAATAAAAAAAACGGCAAGAGCGAGTAATACGCCCCGGCGC‑GGTTGTTTTTCATCCATCGGGAATACTCATTTTTTATTGGTAAAAGCA‑GTTTA‑CCt                                                                                                                                                             >  1:66139/1‑251 (MQ=255)
                                                                           gtgGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTg                                                                                                                                                            <  2:66139/251‑1 (MQ=255)
                                                                                gTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GGAGTAAATCAACTTTAAGTACGCTGGCGCTATACCCCAAAAAAAAAAAAAGGCAAGAGCGGGTAATAAGCCCCGCCGC‑GGTTTTTTTTCATCCATCGGGGATACTCATTTTTTATTTGTTACAACA‑GTTTA‑CCTGCtttt                                                                                                                                                       >  1:321780/1‑251 (MQ=255)
                                                                                     aaTCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTc                                                                                                                                                  <  2:335643/251‑1 (MQ=255)
                                                                                         aGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GAAGGAAATCAACTTGGAGTAAGCTGGCGGTATACCCCCAATAAAATAAAAGGCAAGAGGGAGGAATACGCCCTGCCCC‑GCTTGTTTTGcct                                                                                                                                                                                                          >  1:200405/1‑191 (MQ=255)
                                                                                         aGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGAGGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACGTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGcat                                                                                                                                                                                                          <  2:200405/193‑1 (MQ=255)
                                                                                                  aaaTAGGCCCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGCGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑TTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCa                                                                                                                                      <  2:133198/250‑1 (MQ=255)
                                                                                                  aaaTAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCat                                                                                                                                     <  2:192686/251‑1 (MQ=255)
                                                                                                            aCTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGGAAATCAACTTTAAGTAAGCTGGCGCTATACCCCCAATAAAAAAAACCGCAAAAGCGAGTAATACGCCCCGCCGC‑GTTTGTTTTTCATCCATCGGGAATACTCATTTTTTATTTGTAACAGCA‑GGTTT‑CCTTCTTTTTTTTCTTCAACCTACCCTTTAggtggcg                                                                                                                           >  1:7644/1‑251 (MQ=255)
                                                                                                               ggcggcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGAC‑GGAGGAAATCAAATTTAAGTAAGCCGGCGCTTTACCCCAAAAAAAAAAAACGGCAAAAGCGAGGAATACGCCCCCCCCCCGTTT‑TTTTGCATCCAAAGGGAATAATCCTTTTTTATTGGTAAAAGCA‑GATTA‑CCTGCTTTTTTGTTTTCTACCTACCAAATAggtggaggtg                                                                                                                        >  1:71165/1‑251 (MQ=255)
                                                                                                                     aaaTGCTCATCAGCCCCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAAGCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTG‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTg                                                                                                                  <  2:336/251‑1 (MQ=255)
                                                                                                                     aaaTGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GGAGTAAATCAAACTGACGGACGCTGGGGCCATTCCCCAAAAAAAATAAACGGCAAGAGGGGGTAATACGCCCTGCCCC‑GGTTGTTTTTCATCCATCGGGAATAATCATTTTTTTTTGGTAACAGCA‑GTGTATCCT‑CTTTTTTGTTTTCAAACTACCATCTATGTGGCGGTGGCCCTg                                                                                                                  >  2:116748/1‑251 (MQ=255)
                                                                                                                                  caccacCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GGAGTAAATCAAATTGAAGTAAGCCGGCGCCATACCCCAAATCAAATAAACGGCAAAAGCGAGGAATAAGCCCTGCCGC‑GCTTGTTTTTTCTCCCACGGGGATAATCATTTTTTATTGGTAAAAGCAAGTTT‑‑CCTGCTTTTTTTTCTTCAACCTCCCATATAGGAGGCGGCGGGACTGGGAAAATaaaact                                                                                                     >  1:280083/1‑248 (MQ=255)
                                                                                                                                             aGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATtg                                                                                          <  2:317041/251‑1 (MQ=255)
                                                                                                                                                   aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GGTTGTTTTGCATCCAGCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCAGATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAAt                                                                                    <  2:326126/251‑1 (MQ=255)
                                                                                                                                                    cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCC‑GTTGTAAATCCACTTGAAGCACGCTCGCGCACTACCCCCCATAAAATAAACGGCAAGAACGCGCAAAACGCCCCCCCGC‑GTTTGTTTTTCATCCATCGGGAATACTCCTTTTTTATTGGTAACAGCC‑GGTTA‑CCTGCTTTTTTGTCTTTAAACTCCCCTATAGGTGGCGGGGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAAtt                                                                                   >  1:25290/1‑251 (MQ=255)
                                                                                                                                                         aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GGTGTAAATCAACATGAAGTACGCTGGCGCCATACCCCCAATAAAATAAACGGCAAGAGGGAGGAATACGCCCCGCCGC‑GGTTGTTTTTCATCCCTCCGGGATACTCATTTTTTATTTGTTACAACA‑GTTTT‑CCTGCTTTTTTTTCTTCTAACTACCCTATAGGGGGGGGTGGGCCTGGCGATATAAAGCTGTTCTTCTTTTTGTAATTCAAcc                                                                              >  2:91182/1‑250 (MQ=255)
                                                                                                                                                                 cgATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAAGTCAACGCGggcgac                                                                      <  1:109727/251‑1 (MQ=255)
                                                                                                                                                                            gagttCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGACAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTAt                                                           <  2:58507/248‑1 (MQ=255)
                                                                                                                                                                                tCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGc                                                        <  2:135601/250‑1 (MQ=255)
                                                                                                                                                                                 cATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCt                                                      <  1:240475/251‑1 (MQ=255)
                                                                                                                                                                                            aC‑GTAGTAAATCAACTTGACGTACGATCGCGCTATAACCAAAATAAAATAAGCCGCAAGCGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCT‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCACCGCGGGCGACGGCGACTTTATTGCCTGCACGCAAc                                             >  2:210707/1‑249 (MQ=255)
                                                                                                                                                                                                    aaaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGctacta                                    <  1:91182/251‑1 (MQ=255)
                                                                                                                                                                                                        cAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCCTTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACttt                                                             >  1:274087/1‑222 (MQ=255)
                                                                                                                                                                                                        cAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACttt                                                             <  2:274087/222‑1 (MQ=255)
                                                                                                                                                                                                         aagttGAAGTTCGCTGGCGCTATACCCAAATTGAAATAAGCGGCAAGAGCGATTAATCCTCCCTGCCTC‑TTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑TTTTA‑CATGCTTTTATGTCTTCAACCTTCTATATAGGTAGCGGTGGCACTTGCAATAGAAAGCGTTTCTTCACTGTGTAATTCATCACGGTCGACGGCGACTGTATTGCCTGcccgtaccagggttctattagc                               <  2:245921/248‑22 (MQ=255)
                                                                                                                                                                                                                aGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTaaa                         >  1:127868/1‑250 (MQ=255)
                                                                                                                                                                                                                   aCGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGc                     >  1:34646/1‑251 (MQ=255)
                                                                                                                                                                                                                       tGGCGCTATAACCCCAATAAAATAAGAGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGc                 >  1:115086/1‑251 (MQ=255)
                                                                                                                                                                                                                         gcgcTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGTCTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGccc               <  1:184465/251‑1 (MQ=255)
                                                                                                                                                                                                                          cgcTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGcccc              <  1:177939/251‑1 (MQ=255)
                                                                                                                                                                                                                           gcTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTACTACGCCCCGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCcct             >  1:27182/1‑251 (MQ=255)
                                                                                                                                                                                                                           gcTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCcct             <  1:116748/251‑1 (MQ=255)
                                                                                                                                                                                                                                   ccAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGTCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGgc     <  1:174487/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aagaaaataaGCGGCAAGAGCGAGTAAAACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCa  <  2:371625/248‑1 (MQ=255)
                                                                                                                                                                                                                                         aaaataaGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGgcgc   >  1:180917/1‑247 (MQ=255)
                                                                                                                                                                                                                                         aaaataaGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGgcgc   >  1:390363/1‑247 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
TGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCAC‑GTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGC‑GTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAACAGCA‑GTTTA‑CCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGGCACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCA  >  NC_000913/4003893‑4004372

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