Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,435,707	0	C	A	28.6%	107.0 / 12.9	56	intergenic (‑88/+53)	uspF/ompN	universal stress protein F/outer membrane porin N
Reads supporting (aligned to +/- strand):  ref base C (31/9);  new base A (0/16);  total (31/25)
Fisher's exact test for biased strand distribution p-value = 4.91e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCG  >  NC_000913/1435457‑1435947                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  aaTACGCTAAACCCAGAGAGGTATAGTAGGGCATTGAAGGTATTACCCTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCTATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑TCCAGTTG‑CCACTTGTT‑GTGTAATGACGGAGTTCACATTGTTAATTT‑AGA‑Tc                                                                                                                                                                                                                                                   <  2:689296/251‑1 (MQ=255)   tACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCa                                                                                                                                                                                                                                                  >  1:719306/1‑250 (MQ=255)    aCGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCaaa                                                                                                                                                                                                                                                <  1:325104/251‑1 (MQ=255)     cGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCaa                                                                                                                                                                                                                                                 >  1:185342/1‑249 (MQ=255)      gCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAagg                                                                                                                                                                                                                                              <  1:362718/251‑1 (MQ=255)       cTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAagg                                                                                                                                                                                                                                              >  1:221990/1‑250 (MQ=255)          aaCCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTTCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCTATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑GCAAAggagga                                                                                                                                                                                                                                          <  2:530049/251‑1 (MQ=255)           aCCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTGCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGgaa                                                                                                                                                                                                                                         <  2:15484/251‑1 (MQ=255)                   agGCAGAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGGGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ATACTTTTT‑GTGGAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGc                                                                                                                                                                                                                                 <  1:733554/251‑1 (MQ=255)                       aTAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGa                                                                                                                                                                                                                              >  1:571207/1‑250 (MQ=255)                            agGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGc                                                                                                                                                                                                                        >  1:196438/1‑251 (MQ=255)                             gggCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGcc                                                                                                                                                                                                                       >  2:792347/1‑251 (MQ=255)                                 aGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑Gcc                                                                                                                                                                                                                  <  2:194876/251‑1 (MQ=255)                                  gTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTGATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑Gcc                                                                                                                                                                                                                  >  1:107346/1‑250 (MQ=255)                                           ttACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCgg                                                                                                                                                                                                        >  2:498251/1‑251 (MQ=255)                                            tACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCggg                                                                                                                                                                                                       >  2:771786/1‑251 (MQ=255)                                                     gAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGAt                                                                                                                                                                                              >  2:72415/1‑251 (MQ=255)                                                      aaaTGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGAtt                                                                                                                                                                                             >  2:178646/1‑251 (MQ=255)                                                           aaCCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGa                                                                                                                                                                                         >  2:659691/1‑250 (MQ=255)                                                                  gCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCTTTGATTTAATTCGGAATCGAAAATATCGATAGG‑GAAAAGAATCGTTCTGTTCATAAACCCTCCTGTTTTAGTATCCGAATAAAGTGAAAC‑GCCAGATG‑ACATTTTTT‑GTGTAAGGACGGAGTTCACATTTTAAATTA‑AGA‑TAAAAGGAGGAAAAATAAGAAAAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                                                                                                                                   <  1:74735/249‑1 (MQ=255)                                                                  ccatcatCAATTTTTGCCTCTCCCTAAACG‑GGCATAACCGCGGGTGAATTTAATTCTAAAGCGAAAATATGAATAGG‑GACAAGAAGCGTTCTGTTAAAAAAACCTCCTGTTTTAGTACCCGAATAAAGTGTAAC‑GCCAGAGG‑ACATTTTTT‑GTGAAAGGAGGGAGTTCACATTTTTAATTT‑AGA‑TAAAAGGAGGAAAAAAAAGAAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                                                                                                                                   <  1:684961/248‑1 (MQ=255)                                                                          aaTCTTGGCCTCTCCCTCACGGGGGC‑AAATCACGCGTTGAGTTAATTCTAAATCGAAAATATGGA‑AGGAGAAAAGAATCTTTCTGTAAAAAAACCCTCTTTTTTTAGTACCCGAAAAAAGTATAAC‑GCCAAATG‑ACTCTTTTT‑GTGTAATGACGAAGTCCATTTTTTTAATTA‑AGA‑TAAAGGGAGAAAAAAAAAAAAAAAAAACCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGAcc                                                                                                                                                                         <  2:490378/251‑1 (MQ=255)                                                                               ttcccTCTTCCTAACCTTGGC‑AAATCACGTGTTGATTTAATTTTGAACCGAAAATATCGATAGG‑GAAAAAAACCTTTCTTTTAAAAAACCCGCTTTTTTTAGTACCCGAAAAAATTGTAAC‑GCAAGAGA‑ACATTTTTT‑GTGAAATGACGGAGTTCAC‑TTTTTAATTTTAAA‑TAAAGGGAGAAAGAAAAAGCAAAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaaa                                                                                                                                                                     <  1:784893/247‑1 (MQ=255)                                                                               ttGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GCCAAAAATCGTTTTGTTAAAAAACCCTCTTGTTTAAGAATCCCAATAAAGAGTACC‑GCCAGAGG‑ACACTTTTT‑GGGTAATGAGGGAGTTACCATTTTAATTTT‑AGA‑TCAAAGGAGGAAAAATAAGCAGAAAAAGCCC‑GCAATAAAAGCGGGCAGGAGGATTTAGAACTGATAAACAAGCCCTaa                                                                                                                                                                      >  1:212024/1‑249 (MQ=255)                                                                                 cccttttcccaaaCGTGGC‑AAATCACGCGTTGAGTTAATGCTGAACGGGAAATATCGATAGG‑AACAAGATTCGTTTTGTTCATAAACCCTCTTTTTTAAGTATCCGAAAAAAGTGTAAC‑GCCAGATG‑ACATTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTA‑AGA‑TCAAAGGAGAAAGAAGAAGAAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGc                                                                                                                                                                   <  1:638540/239‑1 (MQ=255)                                                                                     cttctaaaccGTGGC‑TAATCCCGCGTTGAGTTAATTCTAAACCGAAAATTTCGATAGG‑GAAAAAAATCGTTTTGTTCAAAAAACCTCCTTTTTTAGTATCCGCATAAAGTGAAAC‑GCAAGATG‑ACTTTTTTT‑GTGAAATGACGGAGTTCACATTTTTATTTA‑AGA‑TAAAGGGAGAAAAAATAAGAAAAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGAc                                                                                                                                                                <  1:178646/240‑1 (MQ=255)                                                                                                  gC‑TAATCACGCGTTGATTTAATTCTAAATAGAAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAACCCTCTTGTTTTAGTACCCGCATAAAGTGTAAC‑GCCAGATG‑ACATTTTTT‑GTGTAATGACGGAGTTCAAATTTTTATTTA‑AAA‑AAAAAGGAGAAAGAAAAAGAAAAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTaga                                                                                                                                                 <  1:359201/251‑1 (MQ=255)                                                                                                      aatccggcgTTGAGTTAATTCAGAAGCGAAAATATCAATAGG‑AAAAAGAATCTTTCTTTTAATAAACCCCCCTTTTTTAGTATCCGAAAAAAGTGTACC‑CCCA‑ATGGACTTTTTTT‑GTGTAATGACGGATTTCAATTTTTTAATTT‑AGA‑AAAAAGGAGAAAAAATAAG‑AAAAAAAGCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTagag                                                                                                                                                <  1:807241/243‑1 (MQ=255)                                                                                                      aaTCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCCTAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTag                                                                                                                                                  >  1:686712/1‑247 (MQ=255)                                                                                                                tGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTACTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCA‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACACTATCATCGGTAGAGATGCCATTGGc                                                                                                                                     >  2:649462/1‑250 (MQ=255)                                                                                                                     tAATTATGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGAAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGTCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGAAGAGATGCCATTGGCAGCGt                                                                                                                                >  1:471359/1‑250 (MQ=255)                                                                                                                        ttCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGa                                                                                                                             >  2:875104/1‑250 (MQ=255)                                                                                                                            gAATCGAAAATATCAAGAGG‑GAAAAGAATCGTTTTGTTCAAAAACCCTCTTGTTTAAGTTCCCGCAAAAAGTGAAAC‑GCCGGAGG‑ACTTTTTTT‑GTTAAATGCCGGGTTTCACTTTTTTAATTT‑AAA‑CAAAAGGAGAAAAAAAAAGAAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCt                                                                                                                         <  1:771786/250‑1 (MQ=255)                                                                                                                                      tatCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCAt                                                                                                              >  2:380867/1‑251 (MQ=255)                                                                                                                                             tagagACAAGAATCGTTCTGTTCATAAAACCTCATGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGATCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATAGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAAc                                                                                                         >  2:60333/5‑250 (MQ=255)                                                                                                                                               g‑gACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACGGATAAACCAGACCTAAAGCGACAATATCATCGATAGAGATGCCATTGGCAGCGTAGAAGATGTCATCTTCATCCAACAGGt                                                                                                     >  2:42443/1‑251 (MQ=255)                                                                                                                                                              tttttaataaaccttctttttttAGTAACCACAAAAATTGAAAC‑GCCAGTTA‑ACATTTGTT‑GTGGAAGGACGGGTTTCATTTTTTAAATTA‑AGA‑TAAAAGGAGAAAAAAAAAGAAAAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTc                                                                                        <  2:571207/231‑1 (MQ=255)                                                                                                                                                                        acccccctttttaagtacccgaaaaaaGTGTAAC‑GCCAGATG‑ACTTTTTTT‑GTTAAATGAGGGATTTCACATTTTTAATTT‑AAA‑AAAAGGGAGAAAAAAAAAGCAAAAAAGGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTgg                                                                             <  1:792347/227‑1 (MQ=255)                                                                                                                                                                         ccctcctGTTTTATTTCCCGCA‑AAAGTTGAACCGCCAAATG‑ACCTTTTTT‑GTGTAATGACGGAGTCCA‑ATTTTTTATTT‑AAAATAAAGGGAGAAAAAAAAAGAAAAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGa                                                                            <  1:807952/250‑1 (MQ=255)                                                                                                                                                                             cctGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGCTG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATg                                                                        >  2:224287/1‑251 (MQ=255)                                                                                                                                                                                  tttattacccgtaaaaaGTGTAAC‑GCCAGATT‑ATATTTTTA‑GTAAAAGGAGGGAGTTCATTTTTTAAATTA‑AAA‑GCAAAGGAGAAAGAATAAAAAAAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGtttt                                                                    <  2:90161/236‑1 (MQ=255)                                                                                                                                                                                      gTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTg                                                               >  2:424777/1‑251 (MQ=255)                                                                                                                                                                                        acccgaaTAAAGTGAAAC‑GCAAGATG‑A‑ACTTTTTTGTGAAAGGACGGATTTAACATTTATAATTA‑AAA‑AAAAAGAAGAAAAAAAAAAAAAAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGa                                                              <  2:196438/244‑1 (MQ=255)                                                                                                                                                                                          cagaataaatgtaaaC‑GCAAGTTG‑ACTCTTTTT‑GGATAATGTCGGATTTAACATTTTAAATTA‑AGA‑AAAAAGGAGGAAAAAAAAGAAAAAAACGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGaaa                                                            <  2:881510/237‑1 (MQ=255)                                                                                                                                                                                            gCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAAt                                                           >  1:277863/1‑249 (MQ=255)                                                                                                                                                                                              aTAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGt                                                        >  1:670174/1‑250 (MQ=255)                                                                                                                                                                                                            ccAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAat                                           >  1:51151/1‑250 (MQ=255)                                                                                                                                                                                                              aGATG‑ACACTTTTT‑GTGTAAAGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCTC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCCTCCAACAGGTTGATTTTATAGTCCACGTAGGTGGACATGTTTTTATTGAAATAGCAAGTCGCGCCAATATc                                        >  2:200341/1‑251 (MQ=255)                                                                                                                                                                                                               agaaaacttttttt‑GTGAAATGAGGGAGTTCAAATTTTTATTTT‑AGA‑TAAAAGAAGAAAAAAAAAGAAAAAAAACCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATc                                        <  2:670174/242‑1 (MQ=255)                                                                                                                                                                                                                 tG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGc                                     >  1:228397/1‑251 (MQ=255)                                                                                                                                                                                                                  g‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAAc                                                                                                         >  1:402033/1‑182 (MQ=255)                                                                                                                                                                                                                   gacttttttt‑GTTTAAGGAGGTAGTTCACTTTTTTAATTT‑AGA‑AAAAAGGAGAAAAAATAAGAAAAAAACGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGc                                     <  1:498251/245‑1 (MQ=255)                                                                                                                                                                                                                               tgtAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGGGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGAAGGTGGACATGTTTTTATTGAATTAGTAAGTCGCGCAAATATCGGCGTATTTAAcc                           >  2:552321/1‑249 (MQ=255)                                                                                                                                                                                                                                gtAATGACGGATTTCAAATTTTTAATAT‑AGA‑ACAAAGGAGAAAAAATAAGCAAAAAAACCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGa                        <  1:846198/251‑1 (MQ=255)                                                                                                                                                                                                                                   aggagggagtcaacttttttAATTT‑AGA‑AAAAAGGAGAAAAAAAAAAAAAAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGTTCt                     <  1:380867/236‑1 (MQ=255)                                                                                                                                                                                                                                      aCGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCttt                   >  1:461244/1‑250 (MQ=255)                                                                                                                                                                                                                                                      tAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCg  >  2:583473/1‑251 (MQ=255)                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  AATACGCTAAACCCAGAGAGGCATAGTAGGGCAGTGAAGGTATTACCGTCAGGAAATGAACCTCTGCATCATCAATCTTTGCCTCTTCCTCAACGTGGC‑TAATCACGCGTTGAGTTAATTCTGAATCGGAAATATCGATAGG‑GACAAGAATCGTTCTGTTCATAAAACCTCCTGTTTTAGTATCCGCATAAAGTGTAAC‑GCCAGATG‑ACACTTTTT‑GTGTAATGACGGAGTTCACATTTTTAATTT‑AGA‑TCAAAGGAGGAAGAATAAGCAGAAAAAGCCC‑GCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGACATGTTTTTATTGAAATAGTAAGTCGCGCCAATATCGGCGTATTTAACCAGATCTTTATCATCAACACCTGCCG  >  NC_000913/1435457‑1435947

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