Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,1320GA30.4% 81.3 / 11.8 46A16V (GCT→GTT) rarDputative transporter RarD
Reads supporting (aligned to +/- strand):  ref base G (6/26);  new base A (14/0);  total (20/26)
Fisher's exact test for biased strand distribution p-value = 1.62e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAAC  >  NC_000913/4003897‑4004382
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                            
cACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCCCGGAGTAAAACAAATTGAAGTAAGCTGGCGCTATACCCCAAATAAAATAAAAGGCAAGAGGGAGt                                                                                                                                                                                                                                                >  1:228303/1‑250 (MQ=255)
cACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGt                                                                                                                                                                                                                                                <  1:263733/250‑1 (MQ=255)
 aCCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTTAAGTAAGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAAAGCGAgga                                                                                                                                                                                                                                               >  1:306781/1‑248 (MQ=255)
  ccGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAAt                                                                                                                                                                                                                                             <  2:344623/251‑1 (MQ=255)
    gCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                                                                                                                                                                                                                                           <  1:488275/251‑1 (MQ=255)
         gATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGACGTAGTAAAACAAATTGAAGTACGCTGGCGCTATACCCCAAAAAAAAAAAACGGCAAGAGCGAGTAATACGcccc                                                                                                                                                                                                                                      >  1:150809/1‑250 (MQ=255)
            aaaCAGTAGCCAGTTGCGACCAATCAGAACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACGGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCTTCGGCGGGCAAGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGcc                                                                                                                                                                                                                                   <  2:565431/251‑1 (MQ=255)
            aaaCAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATTCCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGcc                                                                                                                                                                                                                                   <  2:319919/251‑1 (MQ=255)
             aaCAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAAACAACTTGAAATACACTGGCGCCATAACCCCAATAAAAAAAAAGGCAAGAGCGAgaaaaacaccccgccc                                                                                                                                                                                                                                  >  1:329966/1‑236 (MQ=255)
                   aGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgttt                                                                                                                                                                                                                             <  1:275949/250‑1 (MQ=255)
                     ccAGGTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGAGCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgtt                                                                                                                                                                                                                          <  2:563859/251‑1 (MQ=255)
                           gCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTAAGCTGGCGCTATACCCCAAATAAAAAAAACGGCAAGAGCGAGTAATACGCCCCGCCGCGTTTTTTTTGcat                                                                                                                                                                                                                    >  1:237049/1‑251 (MQ=255)
                                    cAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGa                                                                                                                                                                                                          <  2:150809/251‑1 (MQ=255)
                                        aCGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGACGGAGTAAATCAACTTTTAGTACGCTGGCCCTATACCCCCAATAAAATAAACGGCCAAAGCGAGGAATACGCCCCGCCGCGGTTGTTTTTcattccttgggaaaac                                                                                                                                                                                                       >  2:399469/1‑238 (MQ=255)
                                        aCGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTTAAGTACGCCGGCGCTTTACCCCCAATAAAAAAAACGGCCAGAGCGAGTAAAACGCCCCGCCGCGTTTGTTTTTCATCCCTC‑GGGGATAc                                                                                                                                                                                                      >  1:395626/1‑251 (MQ=255)
                                        aCGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATAc                                                                                                                                                                                                      <  1:146886/251‑1 (MQ=255)
                                          ggCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGcn                                                                                                                                                                                                                                   <  2:270472/221‑2 (MQ=255)
                                                          aTAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAAc                                                                                                                                                                                    <  2:421862/251‑1 (MQ=255)
                                                              aaaTTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCCAATAAAATAAACGGCAAGAGCGAGTAATACGCCCCGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTTATTGGTAAaagc                                                                                                                                                                                 >  1:317158/1‑250 (MQ=255)
                                                               aaTTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGACGCAGTAAATCAACTTGAAAGACGCTGGCGGTATACCCCCAATAAAAAAAAAGGCAAAAGCGAGTAATACGCCCCGCCGCGTTTGTTTTTCCTCCATC‑GGGAATACTCCTTTTTTATTTGTAAcaac                                                                                                                                                                                 >  1:309/1‑247 (MQ=255)
                                                                  tttCTGGGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAAGGGTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATGCGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGtt                                                                                                                                                                             <  2:163687/250‑1 (MQ=255)
                                                                          gCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCtt                                                                                                                                                                    <  1:384507/251‑1 (MQ=255)
                                                                                    cAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTc                                                                                                                                                          <  1:272919/251‑1 (MQ=255)
                                                                                             tAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCt                                                                                                                                                                                                                                      <  1:192522/167‑1 (MQ=255)
                                                                                                   ggACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAAACAACTTGAAGTACGCTGGCGCTATACCCCCAAAAAAATAAACGGCAAGAGCGAGGAATAAGCCCTGCCGCGGTTGTTTTTCATCCCTC‑GGGAATACTCCTTTTTAATTGGTAACAGCAGTTTACCTGCTTTTTTGTCTTCAACCTACCATATAgg                                                                                                                                           >  1:145004/1‑251 (MQ=255)
                                                                                                     aCCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTTTAGTAAGCCGGCGCTATACCCCAAATAAAAAAAACGGCAAGAGCGAGTAATAAGCCCCGCCGCGGTTTTTTTTCATCCATC‑CGGAATACTCCTTTTTAATTGGTAACAACAGTTTACCTGCTTTTTTGTCTTCAACCTACCATATAggtg                                                                                                                                         >  1:569817/1‑251 (MQ=255)
                                                                                                           ggcggcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAggtggcggtg                                                                                                                                   <  1:319593/251‑1 (MQ=255)
                                                                                                               gcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGGAAATCAACTTGAAGTACGCTGGCGCTAAACCCCCAATAAAAAAAACGGCCAGAGCGAGTAATACGCCCCGCCGCGTTTGTTTTTTATCCCTC‑GGGGATACTCATTTTTTATTGGTAACAGCAGTTTTCCTGCTTTTATTTCTTCAACCTACCCTATAGGGGGCGGGGG‑CAc                                                                                                                              >  1:472861/1‑251 (MQ=255)
                                                                                                                              caccacCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCt                                                                                                               <  2:329966/251‑1 (MQ=255)
                                                                                                                                ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTCGGAAAACAACTTGAAATAAGCTGGCGCTATACCCCAAATAAAAAAAACGGCAAGAGCGAGTAATAAGCCCCCCCGCGGTTGTTTTTCATCCATC‑GGGAAAACTCATTTTTAATTTGTAAAAGCAGTTTTCCCGCTTTTTTTTCTTCAACCTACCAAATAGGGGGGGGGGGGCAC‑GGCAATATAAAgacgt                                                                                                             >  2:318758/1‑247 (MQ=255)
                                                                                                                                ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGGAAATCAACTTGAAAGACGCTGGCGCTATACCCCAAATAAAAAAAGCGGCAAGAGCGAGGAAAAAGCCCTGCCGCGGTTTTTTTTCATCCATCCGGGAA‑ACTCCTTTTTTATTGGTAACAACAGTTTACCCGCTTTTATGTCTTTAACCTACCATATAGGGGGCGGGGG‑CACTGGGAATATAAAACTGt                                                                                                             >  1:312593/1‑251 (MQ=255)
                                                                                                                                ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGGAAATCAAATTGAAGTACGCTGGCGCTATACCCCAAAAAAAATAAACGGCAAGAGCGGGTAATACGCCCTGCCGCGTTTGTTTTGCATCCCTC‑GGGGATACTCATTTTTAATTTGTAACAGCAGTTTACCTGCTTTTTTTTTTTCAACCTACCCTATAGGTGGCGGGGG‑CACTGGCAATATAAAGCTGt                                                                                                             >  1:465186/1‑251 (MQ=255)
                                                                                                                                    caTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGttctt                                                                                                         <  1:75151/251‑1 (MQ=255)
                                                                                                                                                         cACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCgg                                                                                     <  2:329357/250‑1 (MQ=255)
                                                                                                                                                             cgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGggcgac                                                                                <  2:317158/251‑1 (MQ=255)
                                                                                                                                                                   gTCAGGATTTCATCGGCGGGCACGCCGTAAATCAACTTGAAAGAAGCTGGCGCTATACCCCCAATAAAATAAAAGGCAAGAGCGAGGAATACGCCCTGCCGCGGTTGTTTTGCATCCATC‑GGGGAAAATCATTTTTAATTGGGAACAGCAGTTTAACTGCTTTTATGTCTTCAAACTACCATATAGGTGGCGGGGG‑CACTGGCAATATAAAACTGTTCTTCATTGTGTAATTCAACGCGggagccggcgac                                                                          >  1:423611/1‑251 (MQ=255)
                                                                                                                                                                                      gCACGTAGTAAATCAACTTGAAGTACGCTAGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAAc                                                       >  1:214473/1‑251 (MQ=255)
                                                                                                                                                                                            agtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGct                                                  >  1:237134/1‑250 (MQ=255)
                                                                                                                                                                                                         gAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGggcgacggc                                                                             >  2:182249/1‑210 (MQ=255)
                                                                                                                                                                                                         gAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGggcgacggc                                                                             <  1:182249/210‑1 (MQ=255)
                                                                                                                                                                                                         gAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTaa                                    <  1:243078/251‑1 (MQ=255)
                                                                                                                                                                                                           aGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCATTTTACCTGCTTTTATGTCTTCAACCTTCCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAAGGCTGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTTCTAGTAGCAGTaaaa                                  <  2:214473/251‑1 (MQ=255)
                                                                                                                                                                                                            gTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGCGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCACTATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAAc                                 >  2:143870/1‑251 (MQ=255)
                                                                                                                                                                                                               cGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCt                              <  2:145004/251‑1 (MQ=255)
                                                                                                                                                                                                                      gcTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCcct                       <  2:237049/251‑1 (MQ=255)
                                                                                                                                                                                                                         ataCCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACATGCTACTAGTAGCAGTAAAACGCTCGCCcctgc                     <  1:53369/250‑1 (MQ=255)
                                                                                                                                                                                                                                           gCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAAc  <  1:318758/251‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                            
CACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATC‑GGGAATACTCATTTTTAATTGGTAACAGCAGTTTACCTGCTTTTATGTCTTCAACCTACCATATAGGTGGCGGTGG‑CACTGGCAATATAAAGCTGTTCTTCATTGTGTAATTCAACGCGGGCGACGGCGACTTTATTGCCTGCACGCAACAGGCTACTAGTAGCAGTAAAACGCTCGCCCCTGCCTGGGCGCAGATAATCAAC  >  NC_000913/4003897‑4004382

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