Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,585,015	0	T	A	23.9%	109.1 / 11.1	46	I86K (ATA→AAA)	yrhB	putative heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base T (29/6);  new base A (0/11);  total (29/17)
Fisher's exact test for biased strand distribution p-value = 9.28e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TATCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGT  >  NC_000913/3584768‑3585262                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                            tATCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATaa                                                                                                                                                                                                                                                           >  1:370078/1‑250 (MQ=255)  aTCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATaaa                                                                                                                                                                                                                                                          <  2:371800/250‑1 (MQ=255)    cACGACGAATTCGCGAAAGCGAACCATTCCCTTGATGATGCAGATCTACCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACA‑G‑CTGGAAGAATATCTACAGGATTACGAAATaaaaa                                                                                                                                                                                                                                                        >  2:296933/1‑250 (MQ=255)      cgacgCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAgg                                                                                                                                                                                                                                                     >  1:54124/1‑251 (MQ=255)         cgCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑t                                                                                                                                                                                                                                                  >  2:494134/1‑250 (MQ=255)          gCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TAc                                                                                                                                                                                                                                                >  2:240672/1‑251 (MQ=255)             ttCGCGAAAGCGAACCATTACCTTGATGATGCAGAGCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCtt                                                                                                                                                                                                                                             >  2:227289/1‑251 (MQ=255)                   aaaGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCtt                                                                                                                                                                                                                                       <  2:33893/251‑1 (MQ=255)                   aaaGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCtt                                                                                                                                                                                                                                       >  1:544521/1‑251 (MQ=255)                         aaCCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑g                                                                                                                                                                                                                                >  2:51614/1‑251 (MQ=255)                          aCCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑g                                                                                                                                                                                                                                >  1:369330/1‑250 (MQ=255)                                      atgatgCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGACGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTATCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAGAGC‑TACCGTCAGCTTGCCCC‑G‑ACAAAATAGCCct                                                                                                                                                                                                                  >  2:357605/1‑251 (MQ=255)                                             aGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTccc                                                                                                                                                                                                            >  1:401292/1‑250 (MQ=255)                                             aGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCa                                                                                                                                                                                                           <  1:35402/251‑1 (MQ=255)                                              gATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCAc                                                                                                                                                                                                          >  1:286039/1‑251 (MQ=255)                                                          gtcATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGc                                                                                                                                                                                              >  1:111448/1‑251 (MQ=255)                                                                aCTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAAcc                                                                                                                                                                                         <  2:370078/250‑1 (MQ=255)                                                                               tttAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGAGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGtt                                                                                                                                                                        >  1:79230/1‑251 (MQ=255)                                                                                                   cTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCgg                                                                                                                                                    >  2:98824/1‑251 (MQ=255)                                                                                                       ttCGAAGCAGGAGATTTTCTCGAAACTGGAGTGGGGGCCGAGCGTTTAGTGGGTACCGCACTTTTTATTATTGAAAAAGAAAGTGGGAAAATTCTTTCTC‑GGGAAA‑GGGAAAACCCG‑CTGGAAAATTATCTACAGAATTACAAAAAAAAAAGGGC‑TCCTTTCGG‑TTGCCCTTG‑AAAAAATAGCCCTTTTCCCAGGAAGTGGGCCGTTACCTTT‑AGTTCCCCGCCGTTAAATCATCCACCGAGCGCGGGTc                                                                                                                                                 <  1:202865/250‑1 (MQ=255)                                                                                                             gcaggagaTTTTCTCGAAACTGGAGAGGGGGCCGCGCGTTTGGCTGTAAACGCACTTTTTATTATTAATAAAGACAGTGGTAAAATTCATCCTCTGGGAAC‑GGCAAACCC‑G‑CGGGAAGAATTTCTACAGATTTACAAAAAAAAAAAGGC‑TCCCTTCGGCTTCCCCT‑G‑ACAAAATACCCCCCTCCCCACGAAGGGGCCCGCAACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcg                                                                                                                                          <  1:86569/247‑1 (MQ=255)                                                                                                                            gAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCg                                                                                                                           >  2:318568/1‑251 (MQ=255)                                                                                                                                cgggagaggggcccGCGCGCTTAGCTGGTCCCCCCCTTTTTATTATTGAAAAAAACGGTGGTAAAATTCATTCTC‑GGGAACCGGCAAACCC‑G‑CTGAAAAAATACCTACAGGATTACAAAAAAAAAAAGGC‑TCCCTCCGCCTTCCCCT‑A‑CCAAAATACCCCCCTCCCCCCGAAGGGGCCCGCTACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                                                                       <  1:291576/239‑1 (MQ=255)                                                                                                                                        ggggccggggtcTTAGCTGGTACCGCACCTTTTATTATGGAAAAAAACGGTGGAAAAATTATTTCTGTGGGAAC‑GACAAACCC‑G‑GTGAAAAAATATCTACAGAATGAC‑AAAAAAAAAAGGCCTCCTTTCGCCTTCCCCG‑G‑AAAAAATACCCCTCTCCCCACGAAGGGGCCCGCAACCCTT‑GTTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAAcc                                                                                                               <  1:404639/240‑1 (MQ=255)                                                                                                                                          ggCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCAt                                                                                                             >  1:205488/1‑251 (MQ=255)                                                                                                                                                  gcTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCttt                                                                                                     >  1:357795/1‑251 (MQ=255)                                                                                                                                                      aggtGGTACCGCACTTTTTATTTTTAAAAAAGACGGTGGTAAATTTCATTCGATGG‑AAC‑CGAAAAACCCG‑TTGAAAAAAAATCTCCGGGATTAAGAAAAAAAAAGGGC‑TTCTTTCGCCTGCCCCT‑G‑ACAAAATCGCCCTCCCCCCAGGAGGGGGCCCGCAAACCTT‑ATTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                                                                                  <  1:446682/247‑1 (MQ=255)                                                                                                                                                            aaaCGCACTTTTTATTATGAAAAAAAACGGTGTTGAATTTAATTCTC‑GGGAAC‑CGGAAAACC‑GGCGGAAAAAATTTCTACAGGTTTAAAAAAAAAAAAAGCC‑TCCCCCCGGCTTCCCCA‑A‑AAAAAACACCCCTCTCCCCACGAAGAGGGCCGCTACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGcc                                                                                           <  2:357795/250‑1 (MQ=255)                                                                                                                                                                     ttttATTATTAAAAAAGAGAGTGGTAAAATTCATTCTCGGGGACC‑GACAAAACC‑G‑CGGGAAGAATATCTACAGATTTAAAAAAAAAAAAAGCC‑TACCTTCGGCTTCCCCG‑A‑ACAAATTACCCCTCCCCCCCCGAGGAGGGCCCCTACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                                                                   <  2:54124/250‑1 (MQ=255)                                                                                                                                                                             ttGATAAAGACAGTGGTGAAAGTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACACTGCCGTTAAATCATCCACCGAGCGCGGGTCGGCTGCGCCGTACAACTCACCGTCCGGCCCAACCATACAGCTTTGTGTACTGCACATCTCCTCTTTCAgcc                                                                          >  2:258754/1‑250 (MQ=255)                                                                                                                                                                                     gACAGTGGTGAAATTCATTCTCGGGAAAC‑GGAAAACCC‑G‑CTGAAAGAATATCTACAGTATTACAAAAAAAAAAAGGC‑TCCCTTCGGCTTCCCCT‑G‑AAAAAATACCCCTCTTCCCACGAAGAGGCCCGCTACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTc                                                                  <  1:285624/251‑1 (MQ=255)                                                                                                                                                                                      aaaGTGGTAAAATTCTTTCTCGGGGAAC‑GACAAACCC‑G‑CTGAAAGATTATCTACAGGATTACAAAAAAAAAAAGGC‑TCCCGCGGGCTCCCCCT‑G‑AAAAAACCGCCTTCTCCCCACGAAGAGG‑CCGGTTACCTTTAGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt                                                                 <  1:430352/249‑1 (MQ=255)                                                                                                                                                                                         gtggtgAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCTCCTCTTTCAGCGCCACTTTCTGAc                                                              >  2:340750/1‑251 (MQ=255)                                                                                                                                                                                                   tcattcTCGGGGAGC‑GGCAAACCC‑G‑CTGGAAGAATATCTACAGGATTACAAAAAAAAAAAGGC‑TACCTTCGGCTTCCCCT‑G‑ACAAAATACCCCTCTCCCCACAAAGAGGCCCGCTACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt                                                     <  1:353823/250‑1 (MQ=255)                                                                                                                                                                                                               aaC‑GGAAAAACC‑G‑CGGAAAGAATATCTACAGGATTACAAAAAAAAAAAGGC‑TCCCTTCGGCTGGCCCT‑G‑ACAAAATAGCCCTCTCCCCACGAAGGGGCCCGCTACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAg                                        <  1:548709/251‑1 (MQ=255)                                                                                                                                                                                                                   ggCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGt                                     >  1:266571/1‑251 (MQ=255)                                                                                                                                                                                                                      aaaccc‑G‑CTGAAAAAATATCTACATGATTACAAAATAAAAAAGCC‑TCCCTTCGGCTTCCCCT‑G‑AAAAAATCCCCCCCTCCCCACGAAGGGGCCCGCCACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATc                                  <  1:213136/247‑1 (MQ=255)                                                                                                                                                                                                                                tGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCAGGGCTa                           >  2:203285/1‑250 (MQ=255)                                                                                                                                                                                                                                 ggAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTaa                          >  2:234114/1‑250 (MQ=255)                                                                                                                                                                                                                                     gaaTATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCt                     >  1:99220/1‑251 (MQ=255)                                                                                                                                                                                                                                       atatCTACAGGATTACGCAATAAAAAAGGC‑TACCTTCGGCTTGCCTT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCttt                   >  1:487077/1‑251 (MQ=255)                                                                                                                                                                                                                                             ccAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGac             >  2:324718/2‑251 (MQ=255)                                                                                                                                                                                                                                             acagattaaaaaaaTAAAAAAGGC‑TCCATTCG‑TTTTCCCT‑GGACAAAATACCCCTCTCCCCACAAGGAGGGCCTATACCCTT‑AGTCCCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTc                <  2:286039/237‑1 (MQ=255)                                                                                                                                                                                                                                              cAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGac             >  2:413453/1‑250 (MQ=255)                                                                                                                                                                                                                                                       gAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCg   >  1:155715/1‑251 (MQ=255)                                                                                                                                                                                                                                                        aaaTAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGt  >  1:164119/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                            TATCACGACGCATTCGCGAAAGCGAACCATTACCTTGATGATGCAGATCTCCCGGTCGTCATTACTCTACATGGACGCTTTAGCCAGGGCTGGTATTTCTGTTTCGAAGCACGAGAATTTCTCGAAACTGGAGATGAGGCCGCGCGCTTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAAC‑GGCAAAACC‑G‑CTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGC‑TACCTTCGGCTTGCCCT‑G‑ACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTT‑AGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTTCCAGCAGCTTGAGCGTATCCGGGCTAAACCCTTTTTCGACACGCAGCTCGT  >  NC_000913/3584768‑3585262

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