Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,135,781	0	G	T	23.1%	87.4 / 11.4	52	R197R (CGC→CGA)	pitB	metal phosphate:H(+) symporter PitB
Reads supporting (aligned to +/- strand):  ref base G (13/27);  new base T (11/1);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 4.99e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGTAAACC  >  NC_000913/3135533‑3136029                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                             aGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCACCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGTCGACGAAGCCAGCAGGGGCAACCCCCACCAGTACCAGCATTACCAGTCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCTCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCg                                                                                                                                                                                                                                                            <  1:390523/251‑1 (MQ=255)      aGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCtt                                                                                                                                                                                                                                                       <  1:404118/251‑1 (MQ=255)           gCTGTTGCAGGTAGTGTGCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGGCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCATCCgg                                                                                                                                                                                                                                                   <  1:152346/250‑1 (MQ=255)           gCTGTTGCAGGTAGAGTTCAAAGTTGGTAACGGCATCGCGGGTGCGGGTACTTTCATAGCCGGACGGATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGCCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAACCAATCAGCGCAATACGCGTCCAGAATGGCGGATTTTG‑TATTGCTCTTTTTTTTTTTTGcttcttctgg                                                                                                                                                                                                                                                    >  2:103599/1‑242 (MQ=255)           gCTGTGGCATGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTTCGATCTTCCGGa                                                                                                                                                                                                                                                  <  2:36162/251‑1 (MQ=255)                                    ggTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTTAG‑T‑TTGC‑CCTTTTTCTTTTTTCGGTCTTCCGGAATGCGGGGGATA‑AGG‑TCCCGCtttt                                                                                                                                                                                                                       >  1:110341/1‑251 (MQ=255)                                       aaCGGCATCCCGGGGACGGGTAATTTCATAGCCGACCGAATTCATATTGACGACGAAGCCAGCAGGGTCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTc                                                                                                                                                                                                                     <  2:523792/250‑1 (MQ=255)                                          ggCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGATTTCG‑T‑TTTG‑CTTTTTTTTTTTTTTGCTCTTTCTGAGTGCGGTGAATT‑CGG‑GTACGCtttttttctc                                                                                                                                                                                                                 >  2:247637/1‑246 (MQ=255)                                               cgcgGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCt                                                                                                                                                                                                             <  1:588527/250‑1 (MQ=255)                                                      acgggtaATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCATCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGc                                                                                                                                                                                                     <  2:306965/251‑1 (MQ=255)                                                          gtaATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCg                                                                                                                                                                                                 <  1:131997/251‑1 (MQ=255)                                                                tCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTATG‑T‑TTTC‑CTTTTTTTTTTTTTCGAGCTTTCGGAATGCGGTGAATA‑AGG‑TCACGCTTTTTCGTCCCGCTCCCGAAGCGTCg‑ccgc                                                                                                                                                                                            >  1:191188/1‑246 (MQ=255)                                                                        ggACGCATTCATATTGACGACGAAGCCAGCAGGGGCAAGCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCgg                                                                                                                                                                                                                                                   <  2:636392/189‑1 (MQ=255)                                                                              aTTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTATG‑T‑TTTC‑CTTTTTTTTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TGACGCTTTTTTGTCCCCCTCCCGTAGCGTCG‑CAGCCGGAATATCAGGCCTc                                                                                                                                                                            >  1:382111/1‑251 (MQ=255)                                                                              aTTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTAAG‑G‑TTGC‑CTTTTTTTTTTTTTCGATCTTCCGGAATGCGGGGGATA‑CGG‑TCACGCTTTTTTGTCCCGCTCCAGTAGCGGCG‑CAGCAGGAATATCAGGCCTc                                                                                                                                                                            >  1:329592/1‑251 (MQ=255)                                                                                 aatatTGACGACGAAGCCAGACGGGGCAATCCCTCCCAGTCCTAGCATTACCAGGCCGATCCCGTTTTGTCCGTCGTTCGCGCCGGGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTGTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTGCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGTATATCAGGCCTCCCg                                                                                                                                                                         <  2:416305/250‑1 (MQ=255)                                                                                    atTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACA‑G‑TTGC‑CTTTTTTTTTTTTTTGTTCTTCCGGAATACGGGGGATA‑CGGGCCAC‑CTTTTTCGTCCCCCTCCAGGAGCGGCG‑CAGCAGGGATATCAgacatcccccc                                                                                                                                                                       >  1:233949/1‑240 (MQ=255)                                                                                                       ggggCAACCCCCCCCAGTACCAGCATTACCAGGCCGAGCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAgg                                                                                                                                                   <  2:329042/251‑1 (MQ=255)                                                                                                           cAATCCCCACCAGTACCAGCATTACCAGGCCGACCCCTGTTTGGCCGTCGTTCGCGCCGTGCGAAAACGCCACGGCCGCAGCGGAAACAATCAGCGCAATTCGCGGCCATAATGGCGGTTTCCG‑T‑TTGC‑TGTTTTTCTTTTTGCGAGCTTCCGGAATGCGGTGAATA‑CGG‑TCCCGCTTTTTCGTCCCGCTCCAGTAGCGGCG‑CAGCAGGAATATCAGGCCTCTCGCAATGACCAGGCCGACGATAGGGGaa                                                                                                                                               <  1:357036/251‑1 (MQ=255)                                                                                                                       gTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTTAG‑T‑TTGC‑CTTTTTTCTTTTTGCGCTCTTTCGGGATGCGGTGGATA‑CGG‑TTACGCTTTTTTGTCCCCCCACAGTCGCGTCG‑CAGCAGGAATATCTGGGCCCCCCCAATTACCCAGCCGACGATAGGGGGAA‑Caaaaaaag                                                                                                                                    >  2:452531/1‑243 (MQ=255)                                                                                                                        tACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACGGT‑TTGC‑‑TTTTTTTTTTTTGCGATCTTCCTGAATGAGGGGAATA‑CGG‑TCACGCTTTTTTGTCCCCCTCCCCGAGCGTCG‑CAGCAGGGAAATCAAGGCTCCCCCAAAGACCCGGCCGACGATAGGGGAAAACAA‑CAGCgagg                                                                                                                                 >  1:428004/1‑251 (MQ=255)                                                                                                                        tACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTTG‑CTTTTTTTTTTTTTTGGTCTTTCGGGATGCGGGGAATA‑CGG‑TCAATCTTTTTCGTCCCGCCCCCGTAGCGTCG‑CCGCAGGGATAATAGGCCGCCCGCCATGACCAGGCCGACCGCAGGGGAAA‑AAAACAGCgagg                                                                                                                                 >  1:279067/1‑251 (MQ=255)                                                                                                                                    cAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTTCG‑T‑TTGC‑CTTTTTTTTTTTTTCGATCTTCCGGGATGCGGTGGATA‑CGG‑TCACGCTTTTTCGTCCCGCCCCCGTAGCGTCG‑CAGCAGGGATATCAGGCCCCCCGCCATGACCCGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTgg                                                                                                                     >  1:118958/1‑251 (MQ=255)                                                                                                                                     aGGCCTATCCCTTTTTGGCCGCCTTCCGCGCCTTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGt                                                                                                                    <  2:118958/251‑1 (MQ=255)                                                                                                                                               cTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACA‑T‑TTTC‑CCTTTTTCTTTTTTAGATCTTCCGGAAAGCGGTGAAAAACGG‑‑CACGCTTTTTTGTCCCCCTCCAGTAGCGTCGGCGGCAG‑AATATCAGGACCCCCCCAATGACCACGCCGACGCCACGGGGAA‑CAAAAAACGAGGAGAAAATTTTTGTccactcacgc                                                                                                          >  1:316804/1‑242 (MQ=255)                                                                                                                                                 ttttGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCa                                                                                                         <  1:233536/250‑1 (MQ=255)                                                                                                                                                  tttGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑G‑TTTC‑CTTTTTTTTTTTTGCGGTCTTCCGGAATGCGGTGGATA‑CGG‑TCACGCTTTTTCGTCCCCCCCCAGGAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCCAAGACCCGGCCGGCGATTGGGGGAA‑CAAACAGCGAGGGGAAAAATTTTGTCGCTTCACGCAgg                                                                                                       >  1:620486/1‑251 (MQ=255)                                                                                                                                                   ttGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTTTTTTTTTGATCTTCCGGGATGCGGTGGATa                                                                                                                                                                                                                                      >  2:255574/1‑127 (MQ=255)                                                                                                                                                   ttGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATa                                                                                                                                                                                                                                      <  1:255574/127‑1 (MQ=255)                                                                                                                                                    tGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGt                                                                                                      <  1:554203/250‑1 (MQ=255)                                                                                                                                                                        cGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATg                                                                                 <  2:329592/251‑1 (MQ=255)                                                                                                                                                                          aaaaCGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGAt                                                                                  <  1:165479/248‑1 (MQ=255)                                                                                                                                                                             acgccacgccCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTAAG‑T‑TTGC‑CTTTTTTTTTTTTTCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CCGCAGGAAAATTAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGATTAAAGCATCCATCACCggtggaccg                                                                            >  1:44031/1‑243 (MQ=255)                                                                                                                                                                                 cacgccCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTg                                                                                                                      <  2:596083/205‑1 (MQ=255)                                                                                                                                                                                         cAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTTCGATCTTCCGGAATGCGGGGGATA‑AGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGCTAGGGGGAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTTACGCATCCATCACCGAAGAGCCCGGTAACACcgcg                                                                >  1:408853/1‑251 (MQ=255)                                                                                                                                                                                             ggAAACAATCAGCGCAATACGCGTCCAGAAAGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGCCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTgg                                                            <  2:335513/251‑1 (MQ=255)                                                                                                                                                                                                  cAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTaaa                                                       <  2:620486/251‑1 (MQ=255)                                                                                                                                                                                                      cAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGa                                                   <  1:247637/251‑1 (MQ=255)                                                                                                                                                                                                         cgcAATACGCGTCCAGAATGGCGGTTTTCG‑T‑TTTC‑CTTTTTTCTTTTTGCGCTCTTCCGGAATGCGGGGGATA‑CGG‑TCACGCTTTTTTGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGGCTCCCGCCATGAACAGGCCGACGATAGGGGGAA‑CAATCAGCGGGGAGGAAAATTTGGGCACTTCACGCAGGTTTAACGCATCCCTCTACGATGAGCCGGTTAACAGCGCGTTGGGTAAACCCATGc                                                >  1:462639/1‑251 (MQ=255)                                                                                                                                                                                                                 gcgTCCAGAAGGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATc                                        <  1:103599/251‑1 (MQ=255)                                                                                                                                                                                                                       aGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTGCGTCCCGCGCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCa                                  <  2:408853/251‑1 (MQ=255)                                                                                                                                                                                                                           tGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCg                                                                                    >  1:121069/1‑197 (MQ=255)                                                                                                                                                                                                                           tGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGGTTAACGCATCCATCACCg                                                                                    <  2:121069/197‑1 (MQ=255)                                                                                                                                                                                                                            ggCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCa                              >  1:142868/1‑250 (MQ=255)                                                                                                                                                                                                                             gCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAg                            <  1:452531/251‑1 (MQ=255)                                                                                                                                                                                                                               ggTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGgt                          >  1:515063/1‑251 (MQ=255)                                                                                                                                                                                                                                     cG‑T‑TTGC‑CTTTTTTCTTTTGGCGATCTACCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCg                     <  1:176938/250‑1 (MQ=255)                                                                                                                                                                                                                                        t‑ttGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAAc                  >  1:351703/1‑251 (MQ=255)                                                                                                                                                                                                                                             c‑cTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTgg               >  1:67330/1‑250 (MQ=255)                                                                                                                                                                                                                                                 tttttCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAACATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATTAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACCGGAGGc           >  1:555938/1‑251 (MQ=255)                                                                                                                                                                                                                                                      cTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGAAGAGCCGGTTAACAGCGCGTTGGTTAAACCGCTGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGt       >  2:589771/1‑250 (MQ=255)                                                                                                                                                                                                                                                           tGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGTAAAcc  <  1:120366/250‑1 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                             AGTTCAGGATGCTGTTGCAGGTAGTGTTCGAAGTTGGTAACGGCATCGCGGGTACGGGTAATTTCATAGCCGGACGCATTCATATTGACGACGAAGCCAGCAGGGGCAATCCCCACCAGTACCAGCATTACCAGGCCGATCCCTTTTTGTCCGTCGTTCGCGCCGTGCGAAAACGCCACGCCCGCAGCGGAAACAATCAGCGCAATACGCGTCCAGAATGGCGGTTTACG‑T‑TTGC‑CTTTTTTCTTTTTGCGATCTTCCGGAATGCGGTGAATA‑CGG‑TCACGCTTTTTCGTCCCGCTCCAGTAGCGTCG‑CAGCAGGAATATCAGGCCTCCCGCAATGACCAGGCCGACGATAGGGGAAA‑CAATCAGCGAGGAGAAAATTTTGGTCACTTCACGCAGGTTTAACGCATCCATCACCGATGAGCCGGTTAACAGCGCGTTGGTTAAACCGATGCCGATAATCGCACCAATCAAGGTGTGCGAACTGGAGGCCGGTAAACC  >  NC_000913/3135533‑3136029

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