Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,296,381	+GC	intergenic (+587/+55)	gltP → / ← yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,380	1	.	C	94.4%	97.0 / ‑0.3	37	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (8/26);  total (11/26)
Fisher's exact test for biased strand distribution p-value = 7.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01
*	NC_000913	4,296,380	2	.	G	91.7%	82.1 / 3.8	37	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (3/0);  new base G (8/25);  total (11/26)
Fisher's exact test for biased strand distribution p-value = 2.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑A‑‑C‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTT  >  NC_000913/4296193‑4296618                                                                                                                                                                                              ||                                                                                                                                                                                                                                                  gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                                                                                                                                   <  1:133814/251‑1 (MQ=2)     tttATCAGTCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAGCTTAACCGtt                                                                                                                                                                                <  2:211929/249‑1 (MQ=11)             gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTGCCGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTATCGTGATTTTGTAGGTCGGATACGTCGTTT‑A‑‑CGC‑CCAGTCCGACATCAATGAGTGATGCGCCGCTTGCCGTGTCTTATCAGGCCTCTCTTAACCGTTGGTTCAttt                                                                                                                                                                       >  2:245679/1‑251 (MQ=11)                         cAGCGCAATAGCTTGATTTAGCGTGATTTTGTAGGTCGGATCAGGCGTTTATGCCGCATCCGGCATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                                                                                                                           <  1:352256/251‑1 (MQ=14)                               aaTAGCCTGAGTTAGCGTTGTTGTGTAGGTTGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAGTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGATATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc                                                                                                                                                     <  1:173777/251‑1 (MQ=14)                                aTAGCCTGATTTAGCGAGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAAGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAATGCGTTTCAG‑C‑C‑GCA‑TCCGACCTCAATGCGTGATTAGTCGCTGGCCGCGTCTTTTCATGCCTCTCTTACCTGTTGGATCAtttttaacaccctgtacacccg                                                                                                                                                    >  2:299081/1‑234 (MQ=11)                                  aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                                                                                                                                   <  1:396120/250‑1 (MQ=255)                                   gCCTGATGTAGCGTGATTTTGCATGTCGGATAAGGCGTTTATGCCGCATCCGACATCACCGCCCGATGCGATGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTCAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTa                                                                                                                                                 <  2:156864/251‑1 (MQ=14)                                      tGATTTCGCGTGATTTTGTCGGACGAATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAGCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGAGACGCTTGCCGCCTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCCAACTCTTCACACCCGGTATc                                                                                                                                               >  2:421381/1‑250 (MQ=14)                                           tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGATGCGACATTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACGCTTCACACCCGGTATCAAAcc                                                                                                                                          <  2:418730/250‑1 (MQ=17)                                             gcatcattttgtgaggcGGATAAGGGGCTTACGCCCCATCCGACCTCAACTCCTGATGCGACGTTTAACGTGTCGTTTCAGGCCTACGCTAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACTC‑C‑GCA‑ACCGACATCAAGGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACGCTTCACACCCGGTATCAAACCCtt                                                                                                                                       <  2:77448/235‑1 (MQ=17)                                                                   aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCa                                                                                                                 <  2:115318/251‑1 (MQ=255)                                                                     gCGTTTATACCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGt                                                                                                               >  2:356509/1‑251 (MQ=255)                                                                       gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCCAACCCTTCCATACAGCTCAGGTTCAACCAGTGc                                                                                                             >  2:243072/1‑251 (MQ=255)                                                                             gCCGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCCACCAGTGCAGCGcc                                                                                                       >  1:429087/1‑251 (MQ=255)                                                                               cGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACTCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCTACGCTTGCCGCTTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCtt                                                                                                     <  1:342537/251‑1 (MQ=21)                                                                               cGCATCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCtt                                                                                                     <  2:39927/251‑1 (MQ=255)                                                                                  aTCCGACATCAACGCCTGATGCGACGTTTAACGCGGTTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTg                                                                                                  <  2:295304/251‑1 (MQ=255)                                                                                   tCCGACATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGAAATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCACACTCTTCACACCCGGTATCAAACCCTTCCATACAGATCAGGTTCGACCCGTGCAgcgcactttgc                                                                                                 >  2:380222/1‑244 (MQ=18)                                                                                       aCATCAACGCCTGATGCGACGTTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTGCCATACAGCTCAGGTTCAACCAGTGCGGCGCCTTTTGCTTGt                                                                                             <  1:110368/251‑1 (MQ=255)                                                                                               gCCTGATGCGACGTTTAACGCGTCTTATCAGTCCTACGCCAGACAGCGCAATAGTCTGATTTATCGTGATTTTGTAGGTCGGATAAGGCGTGT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGCCGCTTGCCGCGTCGTATCAGGCGTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCg                                                                                     <  1:466533/251‑1 (MQ=255)                                                                                                          cgtttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑TCA‑GCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttc                                                                          <  2:440998/248‑1 (MQ=255)                                                                                                               aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCa                                                                      <  2:229067/250‑1 (MQ=255)                                                                                                                   cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGtt                                                                 <  1:243072/251‑1 (MQ=255)                                                                                                                        tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAg                                                             <  1:160604/250‑1 (MQ=255)                                                                                                                         atccggCCTACGCCTGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGTTAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAAGGCCTGATGCGTCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTTGTTAATTTTTCAAACTCGTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTGTGCTTGTTCTTCTCGATAATACCTTCTTCACCGTTTAAGaa                                                           <  2:388237/247‑1 (MQ=25)                                                                                                                                    gtgaGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAATGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTATCCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGGGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCTTTTAAGAACATCATTCCCg                                                <  1:223583/248‑1 (MQ=37)                                                                                                                                        gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCa                                            <  1:328609/251‑1 (MQ=255)                                                                                                                                                    gCCTGATTTAGGATGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGGTGGTTAATTTTTCAAACGCTTCACACCCGGTATCAAACCCTGCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTTCTGGTTCTTCTCGATAAAACCTTCTTCACCGGTTAAGAACATCAGTCCTGCCCAGGACGCGGAATa                                <  2:200394/251‑1 (MQ=25)                                                                                                                                                        gATTTAGCGTGATTTTGTAGGTCGGAGAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACAGCAATGCCTTATGCGACGCTTGCCGCGTCTTATCAAGCCTATCGTAACCGTGGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGGTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTGCACCGGTTAAGAACATCATTCCCGCCCAGTACTCGGAATAAcc                             <  2:415857/250‑1 (MQ=255)                                                                                                                                                        gATTTAGCGGGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTACCCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCg                            <  2:134361/251‑1 (MQ=255)                                                                                                                                                                gTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCttcttc                                                                       >  1:373842/1‑200 (MQ=255)                                                                                                                                                                gTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCttcttc                                                                       <  2:373842/200‑1 (MQ=255)                                                                                                                                                                gTGATTTTGTAGGTCGGATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGa                    <  2:49807/251‑1 (MQ=255)                                                                                                                                                                                gATAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGt    >  1:223260/1‑251 (MQ=255)                                                                                                                                                                                  tAAGGCGTTT‑ACGC‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGttt  >  1:306221/1‑251 (MQ=255)                                                                                                                                                                                  tAAGGCGTTT‑A‑‑CGCAGCA‑TCCGACATCAATGCATGCTGAGACGATTGCGGAGTCTTATCAGGCATATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTATCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCACAGTACTCGGAATAACCGGTTCGGGAAATTGCGGAGCTGCGttt  >  2:56285/1‑251 (MQ=255)                                                                                                                                                                                              ||                                                                                                                                                                                                                                                  GCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑A‑‑C‑C‑GCA‑TCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTT  >  NC_000913/4296193‑4296618

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