Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	628,868	T→G	27.8%	G106G (GGT→GGG)	entA →	2,3‑dihydro‑2,3‑dihydroxybenzoate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	628,868	0	T	G	27.8%	115.8 / 4.4	64	G106G (GGT→GGG)	entA	2,3‑dihydro‑2,3‑dihydroxybenzoate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (14/27);  new base G (19/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 2.98e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTG  >  NC_000913/628642‑629105                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     aaGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATTGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGTGGGTAACGTCGGCGGGGC‑GGTTTAccct                                                                                                                                                                                                                                         >  2:8431/1‑237 (MQ=255)      aCAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGGCCGATCAGCTCAGTAAAGAGGGCTGGCAGCAGACTTTTGCGGGTAACGTCGGCGGGGC‑GGTTTAACTGT‑Tcc                                                                                                                                                                                                                                    >  1:1151462/1‑238 (MQ=255)             ttGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGTGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAAcc                                                                                                                                                                                                                           >  2:33815/1‑239 (MQ=255)                   aaGCGTTCCCTCAGGAGCAATATCCCTTTGCGACCCGAGTGGTGGGTGTTGCCGACGCTGAGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGGAACGGAGGGACTGGGCGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAAc                                                                                                                                                                                                                     <  1:841032/239‑1 (MQ=255)                    aGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGGGCGACCGATCAGCTCAGTAAAGAGGGCTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGGGC‑GCGTAACCTGT‑TCCAG‑CAAACCATTAAc                                                                                                                                                                                                                     >  1:212128/1‑238 (MQ=255)                    aGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCCGCTCAGTAAAGAGGACTGGCAGCAGACTTTTTCGGGTAACGTCTGCGGGGC‑GTTTTACCTCT‑TCCCG‑CAAACAATGaaac                                                                                                                                                                                                                    >  1:803061/1‑237 (MQ=255)                    aGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGGGG‑GGTTAACCTGG‑TCCAG‑CAAAACATGGAcc                                                                                                                                                                                                                    >  1:538480/1‑239 (MQ=255)                     gCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTTCGGTTAACGTCGGCGGTGC‑GTTTAAACTGT‑TCCAG‑CCAAACATGAACCa                                                                                                                                                                                                                   >  2:989680/1‑239 (MQ=255)                     gCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCa                                                                                                                                                                                                                   >  2:556964/1‑239 (MQ=255)                               aGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTTCGGTTAACGTCGGCGGGGG‑GTTTAACCCGT‑TCCCG‑CAAAACATGAACCA‑GTTTTGCCg                                                                                                                                                                                                         >  2:277169/1‑238 (MQ=255)                               aGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCCCAGGAAAGAGGACAGGCAGCAGACGTTTGCGGTTAACGTCGGCGGGGCGGTTTAACC‑GG‑TCCAA‑CAAACCATGacaaaggttctccg                                                                                                                                                                                                          >  1:784597/1‑225 (MQ=255)                                  aGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCCCAGTAAAAAAGGCTGGCAGCAGACTTTTGCGGGTAACGGCGGCGGGGC‑GGTTGACCCTT‑TTCCG‑CAAAcccaggccccagttcccccctc                                                                                                                                                                                                       >  1:1124042/1‑218 (MQ=255)                                    cAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACCCATGGGGCCGGCCGAGCAGGTCAGTGAAGAGGACTGGCGACAGACCTTTGCGGTTAAAGTCGGGGGGGC‑GTTGAACCTGT‑TCCaa‑ccaacccctaacccacttgcccacccacc                                                                                                                                                                                                    >  1:464280/1‑216 (MQ=255)                                       tatCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCACTGGGGCGACCGATCAGCTTCGTAAAGAGGACTGGCAGCCGACTTTTGCGGTTACAGGCGGGGGGGG‑GGTTTAGCTGTATCCAG‑CAAACAACACAACA‑TTTTCG‑CGTCCg‑cgggg                                                                                                                                                                                               >  1:833563/1‑239 (MQ=255)                                         tCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAg‑cggggcg                                                                                                                                                                                             <  2:1125617/239‑1 (MQ=255)                                         tCCCTGTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAg‑cggggcg                                                                                                                                                                                             <  1:255465/239‑1 (MQ=255)                                              ttGCGACCGAGGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAACCGGAGCGGCGGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGGGTGATCGATCCGCTCAGTAAAGAGGACTGGCAGCGGACCTTTGCGGTTAACGTCGCCCGGGC‑GGTTAACCCGT‑TCCCGCCAACCCTTGAACC‑‑GGTTCGCCGCCCC‑CCgggccgggcg                                                                                                                                                                                        >  1:546850/1‑239 (MQ=255)                                               tGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGGGG‑GTTTAAACTGT‑TTCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGgggcggggcgg                                                                                                                                                                                       >  2:779810/1‑238 (MQ=255)                                                   aCCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGGATTTTACGCATGGGCGCGGCCCCTCAGCTCAGTCAAGAGGACCGGCCGCCGGCTTTTGCGGTTTACGTCGGCGGGGC‑GTTTTACCCTT‑TCCCG‑CCAACCACGAACCAGGTTTCCCCG‑CCG‑ccgggcgggggcggtgt                                                                                                                                                                                   >  1:64806/1‑232 (MQ=255)                                                    ccGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGt                                                                                                                                                                                   <  2:607877/238‑1 (MQ=255)                                                       aaGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAGCGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGGAAAAAGGACTGGCAGCAGACTTTTGCGGTTAACGTCCGCGGGGG‑GGTTTACCCGT‑TCCAG‑CAAACCATGAACCA‑GTTTTGCCCTCCG‑cgggggggggggggtgtcgcc                                                                                                                                                                               >  1:1110695/1‑228 (MQ=255)                                                                    tGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGCATTTTACGCCTGGGCGCGACCGATCAGCCCAGTAAAGAGGACTGGCAGCCCACTTTTTCGGTTAACGTCCGCGGTGC‑GTTTTACCTCT‑TCCCC‑CCAACCCTGAACCC‑GTTTCGCCGTCAA‑CGGGGCGGGGGGGCTGTCACTGTGGccttcgcc                                                                                                                                                                   >  1:1027235/1‑231 (MQ=255)                                                                    tGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCCCAGTAAAGAGGGCTGGCAGCAGACTTTTTCGGGTAACGTCGGCGGGGG‑GTTTAAACCGT‑TCCAG‑CAAAACAATGACCC‑GTTTTCCCGTCAC‑cggggcgggggggttttatactttggcctccga                                                                                                                                                                   >  1:351010/1‑215 (MQ=255)                                                                                      cgcgCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGt                                                                                                                                                <  2:1125612/238‑1 (MQ=255)                                                                                         gcAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCCGACTTTTTCGGTTCACGTCGGCGGGGG‑GGTTTACCTGG‑TCACC‑CAAAACCAGAACCA‑CATTTCCCGTCCG‑CGGGGCGGGGCGCATGTCTCTCTGGC‑GCTCGGcgccccccgcccgcccccgcatg                                                                                                                                            >  1:651070/1‑220 (MQ=255)                                                                                         gcAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTAtt                                                                                                                                             <  1:1189851/238‑1 (MQ=255)                                                                                         gcAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGGTAACGTCGGCGGGGC‑GTTTAACCTGT‑TCCAG‑CAAAACATGAACCA‑GGTTTGCCCTCAG‑CGGGGCGGGGGGGTTGTTCCTGTGGC‑GTCCGACCCCCCGCACCCGCcccggaatg                                                                                                                                            >  2:39395/1‑234 (MQ=255)                                                                                           aaGTGTGTCAGCGACTGTTAGCTGAAACGGGGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCATCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGc                                                                                                                                          <  2:546850/239‑1 (MQ=255)                                                                                              tgtgtCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GGTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGGTTGTTACTGTGGC‑GGCCGGCGCCGCGCACACCCCGCGGATTTGCAtt                                                                                                                                       >  2:793689/1‑238 (MQ=255)                                                                                              cgtgtCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGAATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GCTTAACCCCG‑TCCAG‑CAAACCATGGAACA‑GTTTCGCCCGCAG‑CGGGGGGGGGCGGGTTTGACTGTGGG‑GCCCGGCGCCCCGcaccaccgccgtgctgggc                                                                                                                                         >  1:605910/2‑221 (MQ=255)                                                                                               gtgtCAGCGACTGTTAGCTGAAGCGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATg                                                                                                                                       <  2:1110695/238‑1 (MQ=255)                                                                                               gtgtCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCAt                                                                                                                                        <  2:1151462/237‑1 (MQ=255)                                                                                               gtgtCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATg                                                                                                                                       <  2:746339/238‑1 (MQ=255)                                                                                                tgtCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCCGCAGACTTTTGCGGGTAACGTCCGCGGGGC‑GGTTTACCCTT‑TCCCG‑CAAACCATGACACA‑GTTTTGCCCGCCG‑CGGGGGGGGGCGGCTTTTCCTTTGGG‑GCCCCGCGCCCCCCACACGCCGCGTACTgggcgtgg                                                                                                                                     >  1:743143/1‑233 (MQ=255)                                                                                                 gtCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAg                                                                                                                                     <  2:61539/238‑1 (MQ=255)                                                                                                                                   cgcgCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGc                                                                                                  <  2:1232121/238‑1 (MQ=255)                                                                                                                                           tCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGtt                                                                                          <  1:556964/238‑1 (MQ=255)                                                                                                                                                 cggcggGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGt                                                                                    <  1:823758/238‑1 (MQ=255)                                                                                                                                                    cggGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGG‑GGTTAACCTGT‑TCCAG‑CAAAACATGAACCA‑GTTTCGCCGTCAGCCGGGCCGGGGCGATT‑TCACCTTGGCGGTCCGA‑GCCCCCCACCCCCCCCCGACTGGCCCGGGTGCCTTTGGCGGATCGAAAGCGGCGGCTGAAAGGctcgggctgggggtccgg                                                                                 >  1:1031899/1‑222 (MQ=255)                                                                                                                                                    cggGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCgg                                                                                 <  2:429152/238‑1 (MQ=255)                                                                                                                                                          ttttACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGaa                                                                          <  2:128138/239‑1 (MQ=255)                                                                                                                                                             tACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTg                                                                       <  1:33815/239‑1 (MQ=255)                                                                                                                                                                  aTGGGCGCGACCGATCAGCTCAGTAAAGCGGACTGGCAGCAGACTTTTGCGGCTAACGTCTGCGGGGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGGACCA‑GGTTCGCCGTCAG‑CGGGGCGGGGCGCTTGTCACTGTGGC‑GGCCGGAGCCGCGCACCCGCCGCCGATTGGGATGAGTGCTTATGGGGCCTGGGAAGGGGC‑GCCGACAAAGCTGGCTGTGCGCGGCGCGGCGGAACTgggcgg                                                                  >  1:885983/1‑235 (MQ=255)                                                                                                                                                                  aTGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCggg                                                                  >  2:582694/1‑239 (MQ=255)                                                                                                                                                                    gggCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCCCCCCCGCGCCCCAGCCCCGCATTGGCCTGAGTGCCTTTGGCGCCTCGCAACGCGC‑CCTCAAAAGCCTGGCGTTGCGCGTTGGGCCGGGACTGGCCGGtg                                                                >  1:499016/1‑238 (MQ=255)                                                                                                                                                                            ccGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTg                                                        <  2:99320/239‑1 (MQ=255)                                                                                                                                                                            ccGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTg                                                        <  2:212128/239‑1 (MQ=255)                                                                                                                                                                            ccGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTg                                                        <  2:1144763/239‑1 (MQ=255)                                                                                                                                                                              gATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTgcg                                                      <  2:784597/239‑1 (MQ=255)                                                                                                                                                                                         tAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCCG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACCCGCCGCGGATTGGCATGAGTGCTTTTGGCGCCTCGGAAGCGGC‑GCGGGAAAGCCCGGCGTTGGGCGTCGGGCTGGGACTGGCgggggggggggggcgggctgaggtgg                                           >  1:908213/1‑216 (MQ=255)                                                                                                                                                                                              agGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGc                                                                               >  2:501306/1‑198 (MQ=255)                                                                                                                                                                                              agGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGc                                                                               <  1:501306/198‑1 (MQ=255)                                                                                                                                                                                              agGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGGTGCGGCGGGGGCCGTCAAGTGGgttt                                       >  2:461220/1‑237 (MQ=255)                                                                                                                                                                                              agGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGCGC‑GTTTTACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGGGATTGTGACTGCGGC‑GTCCGCCGCCGCGCACACGCCGCGGATTTGCATGAGGGCTTATGGCGCAACGAGAGCGGG‑GCCGCAAAGCCCGGCGCTGGGCGCCGGGCTGGGACCTGCGGGGCgcgggcgggggcgggggtcgggtgt                                       >  1:1037250/1‑217 (MQ=255)                                                                                                                                                                                                 aCTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGGTTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGGATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAACCTGGGGTTGAGCGTCGGGCTGGAACTGGCGGGGAGCGGGGGgcggcgggattggggtttttc                                    >  1:478266/1‑220 (MQ=255)                                                                                                                                                                                                   tGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGCTTGTCCCCGTGGC‑CTCCGGCGCCGCGCACCCGCCGCGTATTGGCGTGGGTGCTTATGGGGGATCGGAAGCGGC‑GCTCACAAGCCTGGCGCTGCGCGTCGGGCTGGGACTGGCGGGGTGCGGGGcgcgggggggagtggggtttctcgc                                 >  1:340777/1‑218 (MQ=255)                                                                                                                                                                                                          agtcTTTTGCGGTTAGCGGCGGCGGTGC‑GTTTAGCCTGT‑GCCAG‑CAAACCATGACCCA‑GTTGCGCCGTCAG‑CGGGGCGGGGTGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTcc                             <  1:8431/233‑1 (MQ=255)                                                                                                                                                                                                             cTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCcaccga                        <  2:538480/238‑1 (MQ=255)                                                                                                                                                                                                              ttttGCGGTTAACGTCGGCGGGGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCcaccgaca                      <  1:779810/239‑1 (MQ=255)                                                                                                                                                                                                              ttgtGCGGTTAACGGCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCcaccgaca                      <  1:39395/236‑1 (MQ=255)                                                                                                                                                                                                              tgttGCGGTTAACGTCGGCGGGGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCcaccgaca                      <  2:885983/237‑1 (MQ=255)                                                                                                                                                                                                                           gTCGGCGGAGC‑GTATAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAacg         <  2:245435/239‑1 (MQ=255)                                                                                                                                                                                                                             cggcggTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAacgca       <  2:64806/239‑1 (MQ=255)                                                                                                                                                                                                                                cggTGC‑GTTTAACCTGT‑TCCAG‑AAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAacgcacg     <  2:383808/238‑1 (MQ=255)                                                                                                                                                                                                                                   tGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCtg  <  2:651070/238‑1 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     AAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGC‑GTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAG‑CGGGGCGGGGCGATTGTCACTGTGGC‑GTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGC‑GCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTG  >  NC_000913/628642‑629105

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