Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,171,5570CG29.4% 72.2 / 53.0 52intergenic (+13/‑80)rrlB/rrfB23S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (14/22);  new base G (5/10);  total (19/33)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.95e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTT  >  NC_000913/4171401‑4171801
                                                                                                                                                            |                                                                                                                                                                                                                                                       
aCGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTa                                                                                                                                                         <  1:338326/251‑1 (MQ=31)
  gAAACTTTCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAAGTTGCCTGGCGGCAGTAgc                                                                                                                                                       <  2:278271/251‑1 (MQ=14)
   aaaCTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGCGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAgcg                                                                                                                                                      >  2:409098/1‑251 (MQ=14)
        tGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGAAGAAGCGCgg                                                                                                                                                   <  2:164916/249‑1 (MQ=14)
          ccccGAGATGAGTTCTCCCTGACCCTTTAAGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCggtggt                                                                                                                                               <  2:77287/251‑1 (MQ=31)
                       tctcCCTGACCCTTTAAGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGAccc                                                                                                                                   <  1:408625/250‑1 (MQ=31)
                              gACCCTTTAAGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGcc                                                                                                                             <  1:266577/250‑1 (MQ=14)
                              gACCCTTTAAGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGcc                                                                                                                             <  1:390557/249‑1 (MQ=31)
                               ccccTTTAAGGGTCCTGAAGGAACGTAGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATGTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGaa                                                                                                                          <  2:60681/250‑1 (MQ=14)
                                   tttAAGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTc                                                                                                                       >  1:704544/1‑250 (MQ=31)
                                     tAAGGGTCCTGAAGGAACATTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAg                                                                                                                     <  1:363318/250‑1 (MQ=14)
                                      aaGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGaa                                                                                                                   <  1:475371/251‑1 (MQ=31)
                                            ccTGAAGGAACGTTGAAGACGACGACGTTGNAAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGaa                                                                                                              <  2:360661/250‑1 (MQ=11)
                                                    aaCGTTGAAGACGACGACGTTGATAGGCCGGGTGTAGAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAg                                                                                                     <  2:91854/251‑1 (MQ=11)
                                                         tGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGcc                                                                                                 <  1:651687/250‑1 (MQ=31)
                                                             gacgacgacGTTGAGAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAGCCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTTGCGCCGATgg                                                                                            <  2:274914/251‑1 (MQ=14)
                                                                gacgacgTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAg                                                                                                                                                        >  2:62927/1‑188 (MQ=32)
                                                                gacgacgTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAg                                                                                                                                                        <  1:62927/188‑1 (MQ=32)
                                                                 acgacgTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTCCAACGCCGAAGCTGTGTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAgt                                                                                        <  2:167973/251‑1 (MQ=14)
                                                                    acgTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGtg‑tg                                                                                    <  1:184235/251‑1 (MQ=11)
                                                                    acgTTGATAGGCAGGGTGTGTAAGAGCAGCGATGCGTTGAGCTCACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGCAGCTGTTTTGGCGGAT‑GAGAGAAGAGTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTCGCGGCAGTAACGCGGGGGTCACACCTGACCCCATGCCGCACGCCGGAGGGAAAGGCCGCAGCGCCGGTggaaggggg                                                                                     >  2:721195/1‑244 (MQ=11)
                                                                       ttGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGCTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑Tgggg                                                                                 >  1:518486/1‑251 (MQ=14)
                                                                            aGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATTAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTGGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGtctc                                                                             <  2:291131/250‑1 (MQ=14)
                                                                            aGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGtctc                                                                             <  1:67947/250‑1 (MQ=31)
                                                                             ggCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCAGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTccc                                                                           <  1:399475/251‑1 (MQ=14)
                                                                                cGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCAt                                                                        <  1:118501/251‑1 (MQ=31)
                                                                                  ggTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGc                                                                      <  2:224545/251‑1 (MQ=11)
                                                                                      tgtAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACATTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGGAGCGACGATGGTAGGG‑GGGGGTCTCCCCATgaga                                                                    >  2:382228/1‑246 (MQ=11)
                                                                                       gtAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCgag                                                                   >  2:562666/1‑249 (MQ=11)
                                                                                             gcAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGa                                                           <  1:93464/251‑1 (MQ=11)
                                                                                             gcAGCGATGAGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTGTGGCGGAT‑GAGAGAAGATGTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGa                                                           <  2:396281/251‑1 (MQ=2)
                                                                                               agagATGAGTTGAGCTAAACGGTAATAATGAACCGCGAGGCTTAACCTTACAACGCCGACGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAAAGCAAAAGCGGTCTGAATAAAA‑AGAATTAGCCTTGCGGCAGTAGCGCGGTGGACCCACCTGACCCCATGCCGAACTCAGAAGGGAAAGGCCGTAGCGCCGATGGTGGGGGTGGGG‑CTCACCATGCGAGAGTAGGGaa                                                          >  2:594414/4‑250 (MQ=11)
                                                                                                 cGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTg                                                       <  1:586462/251‑1 (MQ=31)
                                                                                                     gCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGAACCCATGCCGAACTCAGAAGTGAAACGACGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGcc                                                     >  1:32652/1‑249 (MQ=11)
                                                                                                      cGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCGGACCCCATGCCGAGCTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAgg                                                  <  1:89065/251‑1 (MQ=11)
                                                                                                       gTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGc                                                 >  1:537836/1‑251 (MQ=11)
                                                                                                       gTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCACCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGc                                                 <  1:398118/251‑1 (MQ=11)
                                                                                                          gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATc                                              <  1:432187/251‑1 (MQ=11)
                                                                                                          gAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCAt                                               <  1:217438/250‑1 (MQ=31)
                                                                                                            gCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGGCCCAACTGCCCCCATGCCGAACTCAGAAGGGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGAAGGGAACGGCCCGGCATCa                                             >  2:164167/1‑250 (MQ=11)
                                                                                                               aaCCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAAt                                          >  1:436668/1‑250 (MQ=31)
                                                                                                                    ttACTAATGACCCGGGAGGCTTAACCTTACAACGCCGAAGCTGGTGTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACCGATTAAATCAAAACGCAGAAGCGGTCTGA‑TAAAACAGGATTTGCCTGGCGGCAGTAGAGGGGAGATCCCCCCTGACCCCATGCCGAACTCCGAAGGGAAACGCCGTTGAGCAGCTGGTAGTG‑TGGGGCCTCCCCATACGAGAACAGGGACAGGCAAGGCATCAAATAAAACg                                    >  2:475258/2‑251 (MQ=2)
                                                                                                                       cTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGaa                                  >  2:527432/1‑250 (MQ=35)
                                                                                                                         aaTGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGAG‑TGGGGTCCCCCCATGCGAGAGTAGGGAACTGCCCGGCATCAAATAAAACGAAAgg                               >  2:385914/1‑251 (MQ=25)
                                                                                                                         aaTGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGAGTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAgg                               <  2:498435/251‑1 (MQ=25)
                                                                                                                                    aGGCTTAACCTTACAACGCCGAAGATGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGaaa                    <  1:409098/251‑1 (MQ=21)
                                                                                                                                            ccTTACAACGCCGAAGCTGTCTTGGCGGAT‑GAGAGAAGATTTTCAGACTGATACAGCTTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGGCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGAGCCTATGGTAGTG‑TGGGGGCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCGCAGTCTGAAGACTGGGc            <  2:417613/251‑1 (MQ=11)
                                                                                                                                                 cAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCGTTc       >  1:361705/1‑251 (MQ=37)
                                                                                                                                                     gCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGt                                                                >  1:217267/1‑190 (MQ=11)
                                                                                                                                                     gCCGAAGGTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGt                                                                <  2:217267/190‑1 (MQ=11)
                                                                                                                                                     gCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGCCTGGGCCTTTCGttt   >  1:536940/1‑251 (MQ=37)
                                                                                                                                                       cGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAACCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGACCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAACGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCCTCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGtttt  >  1:32643/1‑250 (MQ=255)
                                                                                                                                                            |                                                                                                                                                                                                                                                       
ACGAAACTTGCCCCGAGATGAGTTCTCCCTGACCCTTTAAGGGTCCTGAAGGAACGTTGAAGACGACGACGTTGATAGGCCGGGTGTGTAAGCGCAGCGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGCTGTTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGA‑TAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTG‑TGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTT  >  NC_000913/4171401‑4171801

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