Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	523,387	0	A	G	41.7%	38.4 / 11.4	36	M43V (ATG→GTG)	rhsD	protein RhsD
Reads supporting (aligned to +/- strand):  ref base A (10/11);  new base G (15/0);  total (25/11)
Fisher's exact test for biased strand distribution p-value = 6.94e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAGAACAAATAAGTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTG  >  NC_000913/523144‑523630                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           gagaACAACTAAGTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAACGGGAACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGGCCCCTTGTACAGGGTTTGGCAGGGGTGAGAAATGGGGAGCCCACCGGGGTGGCGTGGTTC‑GG‑GTGTCCGGGCGGGGATG‑ACtt                                                                                                                                                                                                                                                     >  1:115188/1‑251 (MQ=255)                         aaTATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑gctgct                                                                                                                                                                                                                         <  1:187750/251‑1 (MQ=255)                                     cATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTACCAATGTCCCGGGGTTGGCAGGTGTAAGAATTGGCGCGCCCCCCTGCGTGGCGTG‑GTCAGGT‑TGTCCGGGAGGG‑GTG‑ACTTC‑GGGCAACCCCGGAAAACCC‑G‑TGCTGGGGGCG‑agggg                                                                                                                                                                                                            >  1:385773/1‑247 (MQ=255)                                            tGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGGCGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGGG‑CTC‑GGGGTGGCCGGGCGGG‑GTG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCCGCTGGGGGCG‑AAGGGGG‑TGccc                                                                                                                                                                                                    >  1:28895/1‑251 (MQ=255)                                                         cATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGGAGGGGTAAGAATTTGCGCGCCCACCGGCGGGGCGTG‑CTC‑GGGGGGCCCGGGCGGG‑GTG‑ACCTT‑GGGCAACCC‑GGCAAAACC‑GCTGCTGGGGGCG‑GAgggggtggccggggggacggggc                                                                                                                                                                                        >  1:203263/1‑230 (MQ=255)                                                                     ttttGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGTGTAAGAATTGGCGCGCCCCCCGGCCTGGCGTG‑CTC‑GGTTTGTCCGGGCGGG‑GTG‑CCTTC‑GGGCCACCC‑CGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCCGCGAGACGGGCCTTGCGCTgcccc                                                                                                                                                                           >  1:128665/1‑250 (MQ=255)                                                                               aaCACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCCTTGTCCAGGGTTCGGCAGGGGGAAGAAATGGGGCGCCCACCCGCGTGGGGTG‑CTC‑GGGGTGGTCGGGGGGG‑GTG‑ACGTCTGGGAAACCG‑GGGAAA‑CC‑GCTGCTGGGGGGG‑GAGATGC‑TGCCCCGCGGGGCGGGCCATGCGCCGCCCCGCCCgccgcc                                                                                                                                                                 >  2:383749/1‑251 (MQ=255)                                                                                    aGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑GTG‑ACTTC‑GGGCAACCC‑GGTAAAACC‑CCTGCTGGGGGCG‑GAGGTGC‑TGCCCGGCGGGACGGGCCTTTCGCTgcccgggcccccgcccgtcc                                                                                                                                                            >  1:158157/1‑237 (MQ=255)                                                                                              aCAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGGGGGTCCGGGCGGG‑GTG‑GCTTC‑GGGCAACCC‑GGGAAATCCCGC‑GCTGGGGGCGGAAGG‑GC‑TGGCCCGCGCGACGGGCCTTGCGCTTCCCCGCCCGCTTCCGTTCATTCTCTcccc                                                                                                                                                  >  1:6595/1‑250 (MQ=255)                                                                                                aaTAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGGCCAGGGGTCGGCCGGGGGAAGAATTGGCGCGCCCACCGGCGTGGGGGG‑CTC‑GGGGTGTTCGGGGGGG‑GAG‑ACCTC‑GGGCAAACC‑GGGAAAACC‑CCTGCTTGGGGGG‑GAGGGGC‑GGCACGGCGGGGAGGGGCTTGCGCTGCCCGGGCCGCTgacgctctttttttcccgcc                                                                                                                                                >  1:283763/1‑232 (MQ=255)                                                                                                           gggCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGGGGCGGG‑CTC‑GGTGTGTCCGGGCGGG‑GTG‑ACTTC‑GGGGAACCC‑GGGAAATCC‑GCTGGTGGGGGGG‑AAGGGGG‑TGCCCGGCGGGACGGGCCCTGCGCTGGCCGGCCCGCCGCCGTTCCTTCTCCCCCCCACCCAcagccgc                                                                                                                                     >  1:93501/1‑248 (MQ=255)                                                                                                                 tttATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTTCCATTGTCCAGGGGTCGGCAGGGGGAAGAATTGGCGCGCCCACCGGCGGGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑GTGGACT‑C‑GGGCAACCC‑GGGAAATCC‑GCTTCTTGGGGCG‑GAGGTGG‑TGCCCGGCCGGACGGACCTTGGGCGGCCCGGCCCGCTGCCCTTCATTCTCTCCCGCACCTCCCGCAGCTCCCgg                                                                                                                               >  1:360868/1‑251 (MQ=255)                                                                                                                   tataggggCGGAAACCCAGGGGCGCGTCAGGGGAATAGGACCCAGTATGGCGG‑CCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGAc                                                                                                                             <  2:158547/243‑1 (MQ=38)                                                                                                                         aGCGGAAAACCAGCGGCGCGTCAGGGAGCTATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGGAGGGGGAAGGATTTGCGCGCCCAACGGAGTGGAGTG‑CTC‑GGGGTGGCCGGGGGGG‑GAG‑ACTTC‑GGGCAACCC‑GGTAAATAC‑GCTGCTGGGGGGG‑AAGGAGCGTGCC‑GAGGAGAAGGACCTTGCCCATCCCGGCCCCCCGCCGCTCGTTCTCTCCCGCCCCCACAGCAACTACCGGACGgaga                                                                                                                        >  1:368428/1‑250 (MQ=255)                                                                                                                          gCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGGTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGGGTGTCCGGGGGGG‑GTG‑ACTTC‑GGGGAACCC‑GGGAAATCC‑GCTGCTGGGGGGGGAAGG‑GC‑TGGCCGGCGAGACGGACCCTGCGCTGCCCCGCCCCCCGCCGCTCATTCTCTCCCCCACCCCCAACAGCCACCCGGCGAAGACg                                                                                                                      >  1:268972/1‑251 (MQ=255)                                                                                                                                      gggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTg                                                                                                          <  2:203263/250‑1 (MQ=255)                                                                                                                                      cggcgCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGGGTGTCCGGGCGGG‑GAG‑ACTTT‑GGGCAACCC‑GGTAAATCC‑CCTGCTGGGGGCG‑AAGGGGC‑TGCCCGGCGAGACGGGCCTTGCGCTGCCCCGCCCGCCGCCGTTCCTTCTCTCCCCCCCCTCAAGCAGCTACCCGACGAAGACGCCCGCACCgggg                                                                                                          >  1:167999/1‑249 (MQ=255)                                                                                                                                         cgcgTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGGGTAAGAATTGGCGCGCCCACCGGCGTGGCGGG‑GTC‑GGTGTGTCCGGGGGGG‑GAG‑AATTC‑GGGGAACCC‑GGGAAAACC‑GCTGCTGGGGGGG‑GAGGGGC‑TGCCCGGCGGGACGGGCCCTGGGCCGCCCGGGCCGCTGCCGGTTATTTTCTCCCGCCCCTCCAGCAGCTACCGGGAGAAGACGCCTGCACCggggggg                                                                                                       >  1:248546/1‑250 (MQ=255)                                                                                                                                                gggAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCCCCGGCGTGGCGTG‑CTC‑GGGGTGTCCGGGCGGT‑ATG‑ACTTC‑GGGCATCCC‑GGTAAATCC‑GCTGCGGGGGGCG‑AAGGCGC‑TGCCCGGCGAGACGGACCTTGCGCGGCCCTGCCCGCTGCCGGTCATTCTCTCCCGCACCTACACCAGCTACCGGATGAAGACGCTGGCACCGTTGGCCGTTTTCg                                                                                                <  2:247468/251‑1 (MQ=255)                                                                                                                                                 ggAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCgg                                                                                               <  2:215393/251‑1 (MQ=255)                                                                                                                                                                    gCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCt                                                                           <  1:493771/251‑1 (MQ=255)                                                                                                                                                                    cgagctccATTGTCCAGTGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCt                                                                           <  1:331124/245‑1 (MQ=255)                                                                                                                                                                          ccATTGTCCCGGGTTCGGCAGGTGTACGAATTGGCGCGCCCACCGGCGTGGAGTG‑CTA‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GATGCTGGGGGCG‑AAGGGGC‑TGCCCGGCGAGACGGACATTGCGCTGCCCGGCCCGCTGCCGTTCCTTATCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCt                                                                           >  2:331124/1‑245 (MQ=255)                                                                                                                                                                            aTTGTCCAGGGTTCGGCAGGAGTAAGAATTGGCGCGCCAACCGGCGTGGCGTG‑CAC‑GGTGTGTCCGGGCGGG‑GTG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGGCCGGCGAGACGGAACTTGCGCCGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCACCCGGGCGAAGACGCCTGCACCGGCGGGCGTTTTCGGCCCCGGCTGGGAAGCGGCCTCCTGatt                                                                    >  2:406864/1‑250 (MQ=255)                                                                                                                                                                              tGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCAACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAGAGCG‑CCTTCTgcaatcc                                                                 >  2:406031/1‑245 (MQ=255)                                                                                                                                                                                  cAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCTTCTGATATCCGCtt                                                             >  1:438075/1‑251 (MQ=255)                                                                                                                                                                                           gCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCTTCTGATATCCGCTTACAGCTACg                                                    >  1:97600/1‑251 (MQ=255)                                                                                                                                                                                                   aaGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCTTCTGATATCCGCTTACAGCTACGTGATGACg                                            >  1:271594/1‑251 (MQ=255)                                                                                                                                                                                                    aGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCCCCGGTGGGCGGTTTCGGCCCCGGCTGGAAAGCG‑CCTTCTGATATCCGCTTACAGCTACGTGATGACgg                                           >  1:10483/1‑251 (MQ=255)                                                                                                                                                                                                            cgcgcCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATAc                                   <  1:406031/251‑1 (MQ=255)                                                                                                                                                                                                                     cGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACa                          <  1:28184/251‑1 (MQ=255)                                                                                                                                                                                                                        cgtggcgtg‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAAcg                       <  2:93501/251‑1 (MQ=255)                                                                                                                                                                                                                              gtg‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACggcgggg                 >  2:34870/1‑250 (MQ=255)                                                                                                                                                                                                                                    c‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGCCGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCTTCTGATATCCGCTTACAGCTACGTGATGAAGGACTGATACTCAACGACAACGGCGGGCGGAgg            >  1:231748/1‑250 (MQ=255)                                                                                                                                                                                                                                    c‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAg             <  1:350773/250‑1 (MQ=255)                                                                                                                                                                                                                                               gggcgag‑ATG‑ACTTC‑GGGCAACCC‑GGTAACTCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGATGCCCGGCCCGCTGCCGTTCCTTCTCTCCCGCCCCTCCAGCAGCTACCGGACGCAGACGCCTGCACCGGTGGGCGCTTTCGGCCCCGGCTGGCAAGCG‑CCTTCTGATATCCGATTACAGCTACGTGATGACGGAATGCTACTCAACGACAAcggcggcgggggccgcccacctt  >  2:220691/1‑234 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                           GAGAACAAATAAGTTTATGTGAAAAATATATAAATACATTAGCTGGTCTTGTGTGTCATTTTATTTTTTTTTGTTGCTAACACAGGGATATGAACAATAACTAAAAGGGCACTTTATATGAGCGGAAAACCAGCGGCGCGTCAGGGAGATATGACTCAGTATGGCGGTCCCATTGTCCAGGGTTCGGCAGGTGTAAGAATTGGCGCGCCCACCGGCGTGGCGTG‑CTC‑GGTGTGTCCGGGCGGG‑ATG‑ACTTC‑GGGCAACCC‑GGTAAATCC‑GCTGCTGGGGGCG‑AAGGTGC‑TGCCCGGCGAGACGGACCTTGCGCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGCTACCGGACGAAGACGCCTGCACCGGTGGGCGTTTTCGGCCCCGGCTGGAAAGCG‑CCTTCTGATATCCGCTTACAGCTACGTGATGACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTG  >  NC_000913/523144‑523630

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