Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,320,516	0	G	A	28.3%	76.7 / 13.9	54	S290F (TCT→TTT)	trpD	anthranilate synthase subunit TrpD
Reads supporting (aligned to +/- strand):  ref base G (32/6);  new base A (2/13);  total (34/20)
Fisher's exact test for biased strand distribution p-value = 2.50e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTGCGGGTTTTCAGTTGCTGGCGAACCGGCATCGCGTGGCGGAATCCGGTGTGATACTTCGGCGCAAAGAGGAAACATACACCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCACGCACCACCGCTGAAAACAGCTGGTGGCTTTCTTGTTGGCTAAGCGT  >  NC_000913/1320271‑1320751                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            gggTGCGGGTTTTCAGTTGCTGGCGAACCGGCATCGCGTGGCGGAATCCGGTGTGATACTTCGGCGCAAAGAGGAAACATACACCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGGAGAAat                                                                                                                                                                                                                                         >  2:68935/1‑250 (MQ=255)                  tGCTGGCGAACCGGCATCGCGTGGCGGAATCCGGTGTGATACTTCGGCGCAAAGAGGAAACATACGCCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTtgctgc                                                                                                                                                                                                                        <  2:112454/250‑1 (MQ=255)                        cGAACCGGCATCGCGTGGCGGAATCCGGTGTGATACTTCGGCGCAAAGAGGAAACATACACCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCAcc                                                                                                                                                                                                                 <  1:344546/251‑1 (MQ=255)                         gAACCGGCATCGCGTGGCGGAATCCGGTGTGATACTTCGGCGCAAAGAGGAAACATACACCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTTCTGCCGTCAccc                                                                                                                                                                                                                >  2:207023/1‑250 (MQ=255)                           aCCGGCATCGCGTGGCGGAATCCGGTGTGATACTTCGGCGCAAAGAGGAAACATACCCCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATCCCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑AGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGcc                                                                                                                                                                                                              <  2:469727/251‑1 (MQ=255)                           aCCGGCATCGCGTGGCGGAATCCGGGGTGATACTTCGGCGCAAAGAGGAAACATACACCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGGTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTGAAAATATTGATACTGTTGCCGCCCTCACCGcc                                                                                                                                                                                                              >  1:132872/1‑251 (MQ=255)                                cATCGCGTGGCGGAATCCGGTGTGATACTTCGGCGCAAAGAGGAAACATACACCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTAc                                                                                                                                                                                                         >  1:79826/1‑251 (MQ=255)                                          cGGAATCCGGTGTGATACTTCGGCGCAAAGAGGAAACATACACCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTTATACTCTTGCTGCCGTCACCGCCAGTACCGACGATATc                                                                                                                                                                                               >  2:290782/1‑251 (MQ=255)                                             aagccGGTGTGATACTGCGGCGAAAAGAGGAAACATACACCAAACACATCCAGCGCCTGGCGGGTTTTAGGGGGGTTCATATAAAGATTAATCCCAAACGCC‑CCCAGCAGATCGGACAAACCAGATTTGCTGGAGCCGCTACGGTTGCCG‑TGTTCCGCCACTTTCAGCCCAAAGGCCGCGGCAACAAACGTTCTGGGGGTATAAATTTTGAACCTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGc                                                                                                                                                                                            <  1:330849/248‑1 (MQ=255)                                                                   cAAAGAGGAAACATACACCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCggg                                                                                                                                                                      >  1:354695/1‑251 (MQ=255)                                                                       gagGAAACATACACCTAACTCATCAAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATAATGTTGCTGCCGTCACCCCCAGGACCGCCGAAATCAGCAAACAGATAATCCCGGCGCGGGAACg                                                                                                                                                                  >  1:370922/1‑251 (MQ=255)                                                                         ggAAACATACACCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGACTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACACACGCACTGGCGGTAGAAATATTGATACCGTTGCTGCCGTCCCCCCCCGTACCGACGATATCAGCAAACAGATAATTCGGGCGCGGGAACggg                                                                                                                                                                >  2:4376/1‑250 (MQ=255)                                                                           aaaCATACACCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGAACCGACGATATCAGc                                                                                                                                                                                            >  1:252535/1‑221 (MQ=255)                                                                               aTACACCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCg                                                                                                                                                                                                                      <  1:316479/191‑1 (MQ=255)                                                                                       aaCTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGAGGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGACGCGGCGACAAACGCACTGGCGGTAAAAATATTTAAAATGTTGCAGGAGTCACCGGCAGGAACCAAGATAAAAGCAAACAGATAACCCGGGCGCGGGAACgggggggcgtgttctcg                                                                                                                                                  >  2:444778/1‑236 (MQ=255)                                                                                       aaCTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACTAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTCTTGCTTCCCTCACTGCCCGTACCGACGCTATCTGCAAACAGATAATCAGGGCGCCCGCAAGGCGCTGCGTTTTcc                                                                                                                                                    >  2:363457/1‑249 (MQ=255)                                                                                       aaCTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTcc                                                                                                                                                    >  1:77870/1‑249 (MQ=255)                                                                                         cTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGt                                                                                                                                                 >  1:283480/1‑250 (MQ=255)                                                                                           cATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCCGTACCGCCGATATCCGCAAACAGATAATCCGGGCGCGGGGACGGCGCTGCGTTTTCCAGTAgc                                                                                                                                              >  2:245636/1‑251 (MQ=255)                                                                                              ccagcgctgggggcgATTTTCCGGCTTTCATATCAAGAT‑AAAACCCAACCCCCGCCACCAGATCGGACGAACCAGTTTTCCTGGAGCCGTTAGGGTTGCCG‑TGTTTCCCCCTTTTCCGCCCCCGGCCCGCGGCGACAAACCCCCTGGCGGTAAAAATATTGATACTTTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAgcg                                                                                                                                             <  1:96931/237‑1 (MQ=255)                                                                                                aGCGCCTGGCGCAATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCggg                                                                                                                                         >  1:420413/1‑250 (MQ=255)                                                                                                       ggCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGcc                                                                                                                                  >  1:54221/1‑251 (MQ=255)                                                                                                            gATTTATCGGCGTTCATATCAAGCTTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTTCTGCCGTCCCCGCCCGTACCGACGCTATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCCGTAGCGCGGATGCTGCCCCggg                                                                                                                             >  2:80552/1‑250 (MQ=255)                                                                                                             ttttATCGGCGTTCATACAAAGATAAATCCCGACCGCC‑GCCAGCAGATCGGACAAACCAGATTTACGGGGGCCGCTACGGTTGCCG‑TGTTCCCCCCCTTTCACCCCACAGCCCGCGGCGACAAACGCACTGGCGTAAAAAATATTGATCCTGTTGCTGCCGTCTCCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCg                                                                                                                               <  1:515368/247‑1 (MQ=255)                                                                                                             aTTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTTATACTGTTGCTGCCCTCACCCCCCGTACCGACGGAATCCGCAAACAGGTAATCCGGGCGCGGGAACAGCGCTGCGTTTTCCCGGAGCGCGGTTGCTGCACCgcg                                                                                                                             >  2:515368/1‑248 (MQ=255)                                                                                                                 acgggCGTTCAAATCAAGATAAACACCGACCGCC‑CCCAGCAGATGGGGCAAACCAAATTTCCGGGAGACGCTACGGTTGCGT‑TGTTCCGCCATTTTCAGCCCACGGGCCGCGGCGACAAACGCATGGGCGGTAAAAATATTGATACTGTGGCGGCCTTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCCGCCCCGGCGAtc                                                                                                                         <  1:415532/247‑1 (MQ=255)                                                                                                                  ccGGCGTTCATAACAAGATTAATACCGACCCCC‑GCCAGAAGACGGGACGAACCATATTTTCGGGAGACGCT‑CGGGTCCCGGTTTTGCGCCTCTTTCAGCCCACGGCCCGCGGCAACAAGCGCACGGGCGGGAAAAATATTGAGCTTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGAtctc                                                                                                                       <  1:391882/250‑1 (MQ=255)                                                                                                                    ggCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATAATGTTGCTGCCGTCAACGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGGACCGCGCCGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGt                                                                                                                     >  2:537325/1‑251 (MQ=255)                                                                                                                          cATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGGACGGCGCTGCGTTTTCCCGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCgggg                                                                                                               >  2:499173/1‑250 (MQ=255)                                                                                                                                aaGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCCCCGCCAGTACCGACGATATTAGCAAAAAGATAATCCCGGCGGGGGAACGGCGCTGCGTTTTTCAGTAACGCGGGTTCTCCCCCCGCGGTCCCGTTAGGGGGCTCAc                                                                                                         >  2:321017/1‑251 (MQ=255)                                                                                                                                  gATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGGTACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGGTGCTGCCCCGGCGAACTCGTTCGGGGGCTCACcg                                                                                                       >  1:550789/1‑251 (MQ=255)                                                                                                                                          cccgaccccccccAGCATATCGGACGAACCAGTTTTACGGGAGACGCTGCGGTTGCCG‑TTTTCCCCCACTTTCACCCCACGGGCCGCGGCCACAACCGCCCGGGCGGTAGAAATATTTATACTGTTGCTGCCGCCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAAtt                                                                                                 <  1:207023/239‑1 (MQ=255)                                                                                                                                             aaCGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCCGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGGATTTTCAt                                                                                            >  2:510925/1‑251 (MQ=255)                                                                                                                                                gcc‑gccAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGGAGAAATATTGATAACGTTGCTGCCGTCACCGCCAGTACCGACGATATAAGCAAACAAATACACCGGGCGCgg                                                                                                                                                                       >  2:198895/1‑173 (MQ=255)                                                                                                                                                    cccAGCAGATCGGAAGAACAAGATTTCCGGGAGCCGCTACGGTTGCCG‑TTTTCCGCCACTTCCACCCCACAGGCCGCGGCGAAAAACTCACGGGCGGTAAAAATATTGATACGTTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCAc                                                                                      <  1:510925/250‑1 (MQ=255)                                                                                                                                                                              aCTGGAGCCGCTACGGATCCCG‑CGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAACTATTGCCAATGTTTATGCCGTCACCCCCAGTACCCCCCATATCAGCAAAAAGATAAACCAGGACCGGGAACGCCGCCGCGGTTTCCAGTAGCGCGGGTGCTGCCCCCGCGCTCGCCTTCGGGGGCTCCCCGCGCATTTTCATTCTCCCCCGcgccgccccccgctgttcacgct                                                            >  2:464109/1‑235 (MQ=255)                                                                                                                                                                                tGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGGTGCTGCCCCGGCGATCTCGCTCGGGGGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTc                                                          >  2:337162/1‑251 (MQ=255)                                                                                                                                                                                 ggAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCCGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCCGCGTTTTCCAGTAGCGCGGGTGCTGCCCCGGCGATCTCGGTCGGGTGCTCACCGCGAATTTTCATGCTCACCCACGCCGCCGCCAGTTGTTCCGGCTTc                                                          >  2:198463/1‑250 (MQ=255)                                                                                                                                                                                     ccgttgggttggcgg‑tgtTTCCCCACTTTCACCCCACAGGCCGGGGCGACAAACGCACGGGGGGTAAAAATATTGATTCTGTGGCGGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTc                                                     <  2:79826/237‑1 (MQ=255)                                                                                                                                                                                           gggTTGCCG‑TGTTCCGCCCCTTCCGGCCCACGGGCCGGGGAGAAAACCGCACTGGCGGTAAAAATATTGTTACGGGTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCAc                                                <  1:450652/249‑1 (MQ=255)                                                                                                                                                                                            ggTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCATCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCACGc                                              >  1:314574/1‑251 (MQ=255)                                                                                                                                                                                                           cccACTTCCCGCCCACGGGCCGCGGCGACAACCGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCACGCACCACCGCTGaaaa                                <  1:104119/250‑1 (MQ=255)                                                                                                                                                                                                              aCTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAAAGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCACGCCCCACCGCTGgaaaa                               >  2:37270/1‑248 (MQ=255)                                                                                                                                                                                                                  ccAGCCCACAGGCCGGGGGAACAAACGCCGGGGGGGTAAAATTATTGATCCTGTGGCTGCCGTCTCCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCACGCACCACCGCTGAAAACAGCtg                          <  1:465626/249‑1 (MQ=255)                                                                                                                                                                                                                      cccACAGGCCGCGGCGACAAACGCACGGGCGGAAAAAATATTGATACTGTGGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCACGCACCACCGCTGAAAACAGCTGGTGGc                     <  1:337162/251‑1 (MQ=255)                                                                                                                                                                                                                         ccgggCCGGGGCGAAAACCGCACGGGCGGAAAAAATTTTGATATTGTTGCGGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCACGCACCACCGCTGAAAACAGCTGGTGGCttt                  <  1:411852/248‑1 (MQ=255)                                                                                                                                                                                                                          cggcccGCGGCAAAAAACGCACGGGGGGTAAAAATATTGATATTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCACGCACCACCGCTGAAAACAGCTGGTGGCTTTc                 <  1:165336/247‑1 (MQ=255)                                                                                                                                                                                                                               cggggcaaaaaaCGCACGGGCGTTAAAAATATTGATTCTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCACGCACCACCGCTGAAAACAGCTGGTGGCTTTCttgtt            <  1:338738/242‑1 (MQ=255)                                                                                                                                                                                                                                  ggcgACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGGTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCCTGCCCACCAGCGCCCCCGCCCGTTGTTTCCGCTTCAGCTCGCCACCCACCACCCCCGAAAACAGCTGGGGGCTTTCTTGTTGGc         >  2:140799/1‑251 (MQ=255)                                                                                                                                                                                                                                     gACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCACGCACCACCGCTGAAAACAGCTGGTGGCTTTCTTGTTGGCTa       >  1:143534/1‑250 (MQ=255)                                                                                                                                                                                                                                     gACAAACGCACGGGCGGTAAAAATATTGATCCTTTCGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCACGCACCACCGCTGAAAACAGCTGGTGGCTTTCTTGTTGGCTaa      <  1:4376/251‑1 (MQ=255)                                                                                                                                                                                                                                      aCAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCAcc                                                                                     >  1:15634/1‑171 (MQ=255)                                                                                                                                                                                                                                      aCAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAACGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCACGCACCACCGCTGAAAACAGCTGGTGGCTTTCTTGTTGTCTaa      >  2:53777/1‑250 (MQ=255)                                                                                                                                                                                                                                         naCGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGCCGATATCAGCAAACAGATAATCCCGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCCCGGGTGCTGGCCCGGCGATCTCGTTCGGGCGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCCTCAGCTCGCCACGCCCCACCCCTGAAAACAGCTGGTGGCTTTCTTGTTTGTTAAGCGt  >  2:260996/2‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            GGGTGCGGGTTTTCAGTTGCTGGCGAACCGGCATCGCGTGGCGGAATCCGGTGTGATACTTCGGCGCAAAGAGGAAACATACACCTAACTCATCCAGCGCCTGGCGCGATTTATCGGCGTTCATATCAAGATTAATACCGAACGCC‑GCCAGCAGATCGGACGAACCAGATTTACTGGAGACGCTACGGTTGCCG‑TGTTTCGCCACTTTCAGCCCACAGGCCGCGGCGACAAACGCACTGGCGGTAGAAATATTGATACTGTTGCTGCCGTCACCGCCAGTACCGACGATATCAGCAAACAGATAATCCGGGCGCGGGAACGGCGCTGCGTTTTCCAGTAGCGCGGTTGCTGCCCCGGCGATCTCGTTCGGGTGCTCACCGCGAATTTTCATGCTCACCAGCGCCGCCGCCAGTTGTTCCGGCTTCAGCTCGCCACGCACCACCGCTGAAAACAGCTGGTGGCTTTCTTGTTGGCTAAGCGT  >  NC_000913/1320271‑1320751

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