Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,040,6230GA63.2% 18.0 / 62.3 58noncoding (107/120 nt)rrfA5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (12/9);  new base A (18/18);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 7.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCC  >  NC_000913/4040378‑4040838
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                         
gATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGCCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCCTGCGAGAGTAGGAAACTg                                                                                                                                                                                                                      >  2:311461/1‑251 (MQ=14)
    cGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑Ag                                                                                                                                                                                                                 <  1:101046/251‑1 (MQ=14)
         aGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATc                                                                                                                                                                                                            <  2:428654/251‑1 (MQ=14)
         aGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTATTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTATCGCCGATGGTATTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATc                                                                                                                                                                                                            <  2:222706/251‑1 (MQ=9)
          gCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCa                                                                                                                                                                                                           <  1:83262/251‑1 (MQ=14)
                cGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATaaa                                                                                                                                                                                                     <  1:464153/251‑1 (MQ=17)
                cGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATaa                                                                                                                                                                                                      >  1:180680/1‑250 (MQ=17)
                 ggTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑a                                                                                                                                                                                                                  >  1:503964/1‑237 (MQ=14)
                 ggTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAAGTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑a                                                                                                                                                                                                                  <  2:503964/237‑1 (MQ=14)
                  gTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCACATaaacc                                                                                                                                                                                                   >  2:160322/1‑249 (MQ=17)
                                    gCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCGGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTATCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACTAATGGCTCAGTCGGAAg                                                                                                                                                                                 <  2:83675/251‑1 (MQ=255)
                                               cAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCttt                                                                                                                                                                      >  1:327571/1‑251 (MQ=38)
                                                aaCGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAATGCTCAGTCGGAAGACTTGGCCTTTc                                                                                                                                                                     <  2:429716/251‑1 (MQ=255)
                                                            gtTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑tgttg                                                                                                                                                       <  2:404947/251‑1 (MQ=14)
                                                             gtttGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGGGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCTTCAAATAAAACGAAAGGCTCATTCGGAAGACTGGGCCTTCCGTTTTATC‑tgttg                                                                                                                                                       <  2:360647/249‑1 (MQ=255)
                                                                tGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGAAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGttgtttg                                                                                                                                                   >  2:261755/1‑251 (MQ=14)
                                                                 ggCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGttgtttg                                                                                                                                                   >  1:186133/1‑251 (MQ=25)
                                                                 ggCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGGTTTATC‑TGttgtttg                                                                                                                                                   <  2:205292/251‑1 (MQ=25)
                                                                           gagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACg                                                                                                                                         <  1:444079/251‑1 (MQ=33)
                                                                               agaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGctct                                                                                                                                     >  1:137456/1‑251 (MQ=33)
                                                                                      cAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAg                                                                                                                               <  1:290396/250‑1 (MQ=18)
                                                                                        gCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGt                                                                                                                              >  1:312424/1‑249 (MQ=33)
                                                                                                cAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATc                                                                                                                    >  2:15923/1‑251 (MQ=18)
                                                                                                 aGATTAAGTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACTAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATca                                                                                                                   <  2:441834/251‑2 (MQ=18)
                                                                                                 aGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCTGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCGCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTGCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATcc                                                                                                                   <  2:548114/251‑1 (MQ=14)
                                                                                                  gaatAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGTCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCGCCTCATGCGAGAGTAGGGAAATGCC‑ATGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTGCGTTGTATC‑TGTTGTTTGTCGGTGAACGGTTTCCTGAGTAGGACAAATccg                                                                                                                  <  2:401446/248‑1 (MQ=11)
                                                                                                    ttAAATCAGAACGCAGAAGCGGTCTGATAAAAAAGAATTTGCCTAGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGACCTCAGAAGTGAAACGCCGTAGCACCGCTGGCAGTGTGGGGTCTCATCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGCAAGGCTCAGTCGGCATACTGGGCCTTTCGTTTTATC‑TATTGTTTGTCGGAGAACGCTCTCCTGAGACGGACAAATccgcc                                                                                                                >  2:331616/1‑251 (MQ=11)
                                                                                                      aaaTCAGAACGCAGCAGCGGTCTGATAAAACAGAATTTGACTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCgg                                                                                                              >  2:390032/1‑251 (MQ=18)
                                                                                                      aaaTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgccg                                                                                                               <  1:502933/250‑1 (MQ=18)
                                                                                                          cAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGc                                                                                                          >  2:137459/1‑251 (MQ=33)
                                                                                                                    aaGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc                                                                                                >  1:56458/1‑251 (MQ=33)
                                                                                                                        gggcTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGGTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGCTTTGCACGTTg                                                                                            >  2:426968/4‑251 (MQ=255)
                                                                                                                           cTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCg                                                                                          >  2:28007/1‑250 (MQ=18)
                                                                                                                               aaaaaCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCGGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCa                                                                                     <  2:3744/250‑1 (MQ=18)
                                                                                                                                aaaaCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCa                                                                                     >  1:1858/1‑250 (MQ=33)
                                                                                                                                     aGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGGCGTAGCGCCGATGGTTGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGcc                                                                               <  2:327571/251‑1 (MQ=18)
                                                                                                                                     aGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGcc                                                                               >  2:200518/1‑251 (MQ=18)
                                                                                                                                       aaTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCg                                                                             <  1:544308/251‑1 (MQ=18)
                                                                                                                                        aTTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCgg                                                                            >  1:211903/1‑251 (MQ=18)
                                                                                                                                           tGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTACGAAGCAACGGCCCGGAgg                                                                         >  2:44779/1‑251 (MQ=18)
                                                                                                                                             ccTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAgggg                                                                       >  1:285055/1‑250 (MQ=18)
                                                                                                                                                ggcggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGCTGGCAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGc                                                                    >  1:91803/1‑251 (MQ=18)
                                                                                                                                                ggcggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTg                                                                      <  2:45712/249‑1 (MQ=18)
                                                                                                                                                   ggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCggg                                                                 <  1:224977/251‑1 (MQ=33)
                                                                                                                                                     cAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCa                                                               <  1:15923/251‑1 (MQ=18)
                                                                                                                                                      aGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGTGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAg                                                              <  2:559579/251‑1 (MQ=18)
                                                                                                                                                       gTAGCGCGGAGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCCATGGTAGTGTGGAGTCTCCTCATGCGAGAGTAGGGACCTGCC‑AGGCATCAACTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAgg                                                             >  1:50638/1‑251 (MQ=255)
                                                                                                                                                         aGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGa                                                            <  1:390032/250‑1 (MQ=18)
                                                                                                                                                         aGCGCGGTGGCCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAc                                                           <  1:311461/251‑1 (MQ=18)
                                                                                                                                                           cgcgGTGGTCCCACCTGACCACATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCCGAAGACTGGGCCTTTCGTTTTCTC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgc                                                         >  1:102786/1‑251 (MQ=16)
                                                                                                                                                                 ggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccgcc                                                    >  1:411400/1‑250 (MQ=18)
                                                                                                                                                                                  tGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATc                                   >  2:360678/1‑250 (MQ=18)
                                                                                                                                                                                            aGAAGTGAAACGCCGTAGCGCCGCTGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGc                          >  1:167685/1‑249 (MQ=18)
                                                                                                                                                                                            aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGGGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGccg                        >  2:53195/1‑249 (MQ=18)
                                                                                                                                                                                                   aaccgccgtggcgccgtggGTAGGGGGGGGTCTCCTAAGGGGGGAGTAGG‑AACCGCCCAGGCACCAAATAAAACGAACGGCTCAGTGGGAAGAC‑GGGCCTTTTTTTTTTTCCTGTTGTTTGTCGGTAAACTCTCTCCTGAGTAGGAAAAACCCGCGGGGGGGGTATTTGAAGTTTGCAAAAAAACGGCCGGGAGGGTGGCGGGCAGGACTCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCa                 <  1:160322/248‑1 (MQ=2)
                                                                                                                                                                                                       gCCGTAGCGCCAATGGGAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCCGTCGGACGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGCGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCACATTAAGCAGAAGGCCATCCt             >  2:221237/1‑251 (MQ=255)
                                                                                                                                                                                                                ccGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGAAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTAACGTATgg    >  2:210621/1‑251 (MQ=18)
                                                                                                                                                                                                                  gatggtgggggggggCCTCCTCATGGGGGAGGAGGAAACGCCC‑AGGCATCAAAAAAAACGAAGGGCTCAGGGGGAAGACTGGCCTTTTCTTTTTATC‑TTTTGTTTGTGGGTGACCGCTCTCCGTAGTAGGAAAAACCCGCGGGGGGCGGATTTGAACGTTGCGAAGCAACGGCCCGGGGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGcc  <  1:198582/240‑1 (MQ=2)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                         
GATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCTCATGCGAGAGTAGGGAACTGCC‑AGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATC‑TGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCC  >  NC_000913/4040378‑4040838

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