Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,396,373	C→T	T103T (ACC→ACT)	amiB →	N‑acetylmuramoyl‑L‑alanine amidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,396,373	0	C	T	93.5%	116.3 / 1.6	46	T103T (ACC→ACT)	amiB	N‑acetylmuramoyl‑L‑alanine amidase B
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (25/18);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 5.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

ATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCG  >  NC_000913/4396144‑4396614                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  atatTCAGGTTTCTATCGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGGTTATGCGTTTAGCCATCAATGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGGTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTACAGATGCACAAACGCTACGGCTGGGGGGC‑GATCTTCCTGAAAACGGAAAAACCGAAGCg                                                                                                                                                                                                                              >  2:140234/1‑251 (MQ=255)      cAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGa                                                                                                                                                                                                                          >  1:396786/1‑250 (MQ=255)         gTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCg                                                                                                                                                                                                                      >  1:127509/1‑251 (MQ=255)               aaCGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCa                                                                                                                                                                                                                   <  2:396786/248‑1 (MQ=255)               aaCGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGaa                                                                                                                                                                                                                <  1:67551/251‑1 (MQ=255)               aaCGGTAAGTAACAGGCGCGGATAACGTTGAGTTGTATTGGCGATCATGATTATGCGTTTAGCCATCAAAGCAAACGCACCGGGGAGCTAGATATCAAACAAACGGGCGTGATTCAGGGTCTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAATCGGTAAACCCTGAGCGGTGAAGCGGCAg                                                                                                                                                                                                                  <  2:156470/249‑1 (MQ=255)               aaCGGGAATCAACAGGCGCGGCTACCGTTCAGTTTTATTGGCCATCCTGATTATGCGATTAGCCATCAACGCAAACGCACCGTGGCGCTCGCTATCATACAAACGGGCGTGATTCAAGGCCTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGTACGCCTAACGATGCACACACGCTCCGGCTGGGGGTCAGA‑CTTCCTGAAAACGGTAAAACCGAAgacgagaaacagaaga                                                                                                                                                                                                                 >  2:246591/1‑235 (MQ=255)                      aTCAACAGGAGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATAAACCAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAACGATGCACAAACGCTACGGCTGGGGGTG‑GATCTTACTGAAAACGGAAAAACCGAAGCGGGGAAGTGGCAGAATGgcagc                                                                                                                                                                                                         >  2:409234/1‑251 (MQ=255)                              acgtgttctctgccgagctgtataggCGATCCGGATTATGCGTTTAGCCACCAACGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCCGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTCAAGATGCATAAACGCTTCGGCTGGTTGAC‑GATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTacac                                                                                                                                                                                                 <  1:64275/227‑1 (MQ=255)                                  ggATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTgtc                                                                                                                                                                                             >  1:376652/1‑251 (MQ=255)                                             aGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCACAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAAAGGAAAAACCGAAGCGGGGAAGCGGGAGCATGGCCGCAAATTCACTTTCGTCTTTACGAt                                                                                                                                                                                  >  2:257257/1‑251 (MQ=255)                                                        cGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCCTTACGATTAACGCCGAtg                                                                                                                                                                       >  1:139357/1‑251 (MQ=255)                                                             cTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAACGGCAGAATGGCAGCAATTACCCTGTCGTCTTTTCGATTAACGcctagtgtccc                                                                                                                                                                  >  2:109480/1‑242 (MQ=255)                                                               gATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTgccgcc                                                                                                                                                                >  1:380471/1‑251 (MQ=255)                                                                      catttAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGcct                                                                                                                                                         <  1:89629/249‑1 (MQ=255)                                                                               aTCAAAGCAAACGCACCGTGGCGCTCGGTATCAAACAATCGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GAACTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGgcag                                                                                                                                                                                                          >  2:243092/1‑193 (MQ=255)                                                                               aTCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGgcag                                                                                                                                                                                                          <  1:243092/193‑1 (MQ=255)                                                                               aTCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGAACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTAAACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTccgccg                                                                                                                                                >  1:1184/1‑251 (MQ=255)                                                                                 cAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGAGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCgccgcc                                                                                                                                              >  1:232339/1‑251 (MQ=255)                                                                                       aaaCGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGt                                                                                                                                        >  1:261228/1‑251 (MQ=255)                                                                                                                 acaaacGGGCGTGATTCAGGGACTGCCGATGTTGTTCAGCGGCAATAATATGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGTTTAACGCCGATGTGCCGACACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCgg                                                                                                              >  1:358068/1‑251 (MQ=255)                                                                                                                  caaacGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGt                                                                                                             >  1:230835/1‑251 (MQ=255)                                                                                                                      cgggggggATTCAGTGCCTGCCGTTGTTGTTCAGCGGCAATATTCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGTTGGTGGTC‑GATCTTACTGAAAACGTTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTc                                                                                                         <  2:166101/244‑1 (MQ=255)                                                                                                                                 cggggtctgctgTTGTTGTTCAGCGGCAATAATCTGGTGAAGTCGTTCCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCGGGTGGTT‑GAGCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTATCGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGTAACGCCTGCGTTTGTCTTACCGCGCGt                                                                                              <  2:416785/245‑1 (MQ=255)                                                                                                                                 cagggccggCCGTTGTTGTTCAGCTGTAGTAATCTGGTGAATGAGATTCGCTCTGGCACGCCTAAATATGCAAAATCGCTACGGCTGGGTTTC‑GATTTTACTGAAACCGGTAAAGCCGAAGCGGTGAAGCGGCAGATTGCCTCTACGTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACTCGTTGACACGCCTGCGGTTGTCGCACcgcgcg                                                                                               <  2:314252/242‑1 (MQ=255)                                                                                                                                     gacggccgtttgtgttCAGCGGCAATAATCTGGTGAATGCGATCCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGGCAGc                                                                                          <  1:246591/239‑1 (MQ=255)                                                                                                                                     gACTGACGTTGTTGTTCAGAGGCAATAATCTGGTGCAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGGAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGCTGTGCCGCCACCGCCTCCTCCGCCGCCCGAGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTccgg                                                                                          >  2:265352/1‑248 (MQ=255)                                                                                                                                     gACTACCGTTGTTGTACAGTGTCAATAATCTGGTGAAGGCGATTCGTCCTGGAACGCCGAAAGATGCACAACCGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGTAGCGGTGAAGCGGCAGAATGGCAGCAATTCCACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCTCCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGc                                                                                          <  2:412917/251‑1 (MQ=255)                                                                                                                                                  gttCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGCAACGCCTGCGGTTGTCGCCCCGCGCGTCAGCGAACCGgcgagc                                                                              >  2:179697/1‑247 (MQ=255)                                                                                                                                                     cAGCGGCAAAATTCTGGGGAAGGCGATGCGCTCGGGAACGCTTAAAGATGCACAAACGCTACGGCTGGGGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGAGGCAGAATGGCAGCAATTACACTTTCGTCTTTACGATTAACGCCGATGTGCCGCCCCCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTGGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCTCAATc                                                                          <  2:261514/251‑1 (MQ=255)                                                                                                                                                             ataataTGGTGAAGGAGATGCGCCCTGGATCGCCTAAAGAGGCCCAAACTCTACTGCTGGTTGTC‑GATCTCACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCCATTACACTTTCTGCTTTACGATTAACGCCGATGTGCCCCCACCGCCGCCTCCGCCGCCCTGGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTaaa                                                                  <  2:182987/251‑1 (MQ=255)                                                                                                                                                               aatCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAAc                                                                <  1:109480/251‑1 (MQ=255)                                                                                                                                                                             gATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGACAACGGTAAAACCGAAGCGGGGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCAGATGGGCCGCCACCGCCACCTCCGCCGCCCGTGGTTGCAAAACGCGTTGAAACTCCTGCGGTTGTCGCCACGCGCCTCAGAGCACAGGCGCGCAATCCGCTTAAAACGGAAAGGAACCGcc                                                  >  2:228915/1‑250 (MQ=255)                                                                                                                                                                                 cGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTAc                                              <  1:134223/251‑1 (MQ=255)                                                                                                                                                                                   ctctGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGACACGCCTGCGGTTGTCGCACCGCGCGTCAGCGACCCGGCGCGCCATCCGTTTTAAACGGAAAGGACACGCACTACgg                                            >  2:85516/1‑251 (MQ=255)                                                                                                                                                                                       ggAACGCCTAAAGATGCACAAACCCTACGGCGGGGGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGgt                                          <  1:41196/249‑1 (MQ=255)                                                                                                                                                                                          aCGCCTAAAGATGCACAAACGCTACGGCTGGTGGTA‑GATCTTACTGAACACGGTACAACCGCAGCGGTGAAGCGGAAGAATGGCAGCAATTACACTGTCGTCTTTACGATTACCCCCGATGTGCAGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGACACGCCTGCGGGTGTCGCCCCGCGCGTCAGGGAACCGGGGGGCACTCCGTTTAAAACGGAAAGTAAcctaacaccggggggtat                                     >  2:344481/1‑235 (MQ=255)                                                                                                                                                                                           cGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACAGCCTCCTCCGCCGCCCGGAGTTGCGAAACGCGTTGAAACGCCGGCGGGTGTCGCCCCGCGCGGCAGCGAACCGGCGGCCAATCCGGTTAAAACGGAACGTAACCGAACTACGGGTGTTATc                                    >  2:93478/1‑251 (MQ=255)                                                                                                                                                                                               tAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATcagc                                 >  1:279932/1‑250 (MQ=255)                                                                                                                                                                                                      gCACAAACGCTAAGGCGGGGTTTC‑GATTTTACTGAAAACGGTAAACCCGAAGCGGTTAAGCGGCAGTATGGCAGCAATGACTCTGTCGTCGTTCCGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCTAAACGCGTTGAAACGCCTGCTGTTGGCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGGAACCGCACTACGGGTGTTATCAGCAGTAATa                          <  2:358068/250‑1 (MQ=255)                                                                                                                                                                                                         aaaaCGCTACGGCTGGTGGTC‑GATCTTTCTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGt                      <  2:279932/250‑1 (MQ=255)                                                                                                                                                                                                              gCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTCCACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTgg                                                                                                                                         <  2:107337/131‑1 (MQ=255)                                                                                                                                                                                                              gCTACGGCTGGTGGTC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTgg                                                                                                                                         >  1:107337/1‑131 (MQ=255)                                                                                                                                                                                                                        gtggtC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACcg                                                                                                   >  1:214121/1‑159 (MQ=255)                                                                                                                                                                                                                        gtggtC‑GATCTTACTGAAAACGGTAAAACCGAAGCGGGGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTCCGATTAACGCCGATGTGTCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACcg                                                                                                   <  2:214121/159‑1 (MQ=255)                                                                                                                                                                                                                             c‑GATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGAGTTATCAGCAGTAATACGGTAACGCGTCCGGCAgcgcgcg  >  2:126608/1‑251 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  ATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTC‑GATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCG  >  NC_000913/4396144‑4396614

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