Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,060	C→T	20.9%	intergenic (+266/+376)	gltP → / ← yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,060	0	C	T	20.9%	44.2 / 15.0	43	intergenic (+266/+376)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (21/13);  new base T (5/4);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.84e-01

ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCA  >  NC_000913/4295833‑4296274                                                                                                                                                                                                                                    |                                                                                                                                                                                                                        aTGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGTTTTTGTAGGTCGGATAAGGCGTTTACGCCCCAACCGACATCAAc                                                                                                                                                                                                                         >  1:121969/1‑228 (MQ=12)        aTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGAt                                                                                                                                                                                                                 >  1:148376/1‑228 (MQ=14)             aCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGATAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAT‑GCCTGATGcta                                                                                                                                                                                                             >  1:205486/1‑225 (MQ=14)                   gCCGCTGCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGc                                                                                                                                                                                                      <  1:4298/228‑1 (MQ=17)                        tttttctcgggCCGACGGCAGACAGCGCAATAGCCGTATTTAGCGGGATTTTGTGGGGCGGATAATGCGTTTATGCCGCATCAGACATAAACGCCTGATGCGACGTGTAACGCGTCTTATCAGGGCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAT‑GCCTGATGCGACGCTTGCCGCGt                                                                                                                                                                                                 <  2:148376/219‑1 (MQ=11)                        gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAATGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCCGCCAGCGCAATAGCCTGCTTTACCGTGATTTTGTAGGTCGGATACGGCGTTTACGCCCCACCCGACA‑CAACTGCCTAATACGACGATTGcccctgc                                                                                                                                                                                                >  2:377914/1‑224 (MQ=11)                             aTCAGGCCGACGCGAGACAGCGCAATGGCCTGATTTAGCGTGATTTTGTGGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAT‑GCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                                                                                            <  1:116570/228‑1 (MQ=255)                                 ggCCTACGCCAGACCGCGCAATAGCCTGATTTAGCTTGATTTTGGAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                       <  1:223033/229‑1 (MQ=21)                                  gctcagcccAGAAGGCGCAATGGCCGGATTTTGCGGTATTTTGTGGGTGGGATAAGGCGTTTACGCCTGATGCGCAAGCAATGCCTGATGCGAAGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                       <  2:195034/221‑1 (MQ=11)                                  gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                                                                                       <  1:454325/228‑1 (MQ=34)                                      aCGCCAGACAGCGCAATAGCCTGATTTGGCGTGATTTTGAAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGGCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACg                                                                                                                                                                                  <  1:152563/229‑1 (MQ=21)                                      aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAT‑GCCTGATGCGACGCTTGCCGCGTCTTTTCAGGccattc                                                                                                                                                                                  >  1:6762/1‑225 (MQ=18)                                       cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTTCCGCCGCATCCGACATCAAT‑GCCTGATGCGACCCTTGCCGCGTCTTATTAGGCATatct                                                                                                                                                                                 >  2:690258/1‑226 (MQ=255)                                            gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGCTTTTGTAGGTCGGATCAGGCGTTTACGCCGCATCCGCCATCAAT‑GCCTTGTGCGACGCTTGCCGCGACTTATCAGGCCTatcttaccc                                                                                                                                                                            >  1:355016/1‑221 (MQ=11)                                                 cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGCATAAGTCGTTTACGCCGCATCCGACATCAAC‑GCCGGCTCCGACGCCTGCCGCGGCCATTCATGCCTACGCCaacccccca                                                                                                                                                                       >  2:41490/1‑221 (MQ=11)                                                 cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGTACTCACC‑GCCTACTGCGACGTTTGCCGAGGCTTTAAAGGCCTCCCCAACACAgcgc                                                                                                                                                                       >  1:532236/1‑229 (MQ=11)                                                       aacctattttAGCGTGATTTTGTGGTGCGTATACGGAGTTTATGCCGCGTACACCATCAACGCCTGATGCGACGCTTACTGCGTCTTATCAGGCCTACGCGAGACAGCGAAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                 <  2:228334/222‑1 (MQ=11)                                                                                  cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGTGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCgg                                                                                                                                      >  1:414090/1‑229 (MQ=21)                                                                                  cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATCCGACGCTTGCCGCGTCTTATCAGGCCTACCCCCACCAGCGCATAAACTTTCATTAGCGTGATTTTGTAGTTCgg                                                                                                                                      >  2:135650/1‑229 (MQ=11)                                                                                               tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTCTCTGGCCTACGCCAAACAGCGCAAAAGCTTGCTTTACCGCGCTTTTGGAGGACGGAAAACACGTTTac                                                                                                                         >  2:272117/1‑228 (MQ=11)                                                                                               tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCCACCAGCGCATACACCTGATTTATCCGGATTTTGTAGCTCGGTAAAggggatact                                                                                                                         >  2:573334/1‑222 (MQ=11)                                                                                               tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACCCCCGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGACAAGGccgtaaa                                                                                                                         >  2:748880/1‑223 (MQ=255)                                                                                                      cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCa                                                                                                                   >  1:65926/1‑228 (MQ=37)                                                                                                                cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAT‑GCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                                                                                            >  2:9338/1‑145 (MQ=25)                                                                                                                cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAT‑GCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                                                                                            <  1:9338/145‑1 (MQ=25)                                                                                                                                          tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGCTTTAGCGTGCTTTTGTATGTCGGATCAGACGTTTCTGCCGCATCCGCCTTCAAATCCTGTTGCGACGCTTCAccctcc                                                                               >  1:352085/1‑223 (MQ=11)                                                                                                                                            ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCTTGATTTAGCTTGATTTTTTATGTCGGATAAGTCGTTTAACCCACATCCGCCATCACAGCCGGATGCTACCCTTTACGCGTCttt                                                                            >  2:431064/1‑228 (MQ=11)                                                                                                                                                             gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                                                                                             <  2:423171/115‑1 (MQ=39)                                                                                                                                                             gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGCTTATGCCGCCTCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTCTCAGGCCTACGCCAgg                                                            >  1:493090/1‑227 (MQ=19)                                                                                                                                                                  cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATCAGGCGTTTATGCCGCATCCGCCATCCACGCCTGATGCGACGCTTAACGCGCCTTTTCCAGCCAACGCCAAAccacga                                                      >  2:194917/1‑224 (MQ=11)                                                                                                                                                                       tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                                                                                             <  2:306241/105‑1 (MQ=40)                                                                                                                                                                        ggcaggATTTGGCGTTATTTTGGAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCACC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                <  1:704490/224‑1 (MQ=255)                                                                                                                                                                            tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAT‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                      <  2:616966/91‑1 (MQ=255)                                                                                                                                                                            tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                                                                                      >  1:616967/1‑91 (MQ=40)                                                                                                                                                                            tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                 <  2:271862/112‑1 (MQ=39)                                                                                                                                                                                tttgtttGATTTTGTGGGTAGGGTAGGGGTTTTTCGCCGCAGCCGCAATCGAC‑GCCTGATGCTACGGTTGCCGCGTCTTATCAGGCCTAAGTCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                        <  2:202758/223‑1 (MQ=11)                                                                                                                                                                                                 gTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATACGGCCTACGCCAAACACCGCAAAAACCTGATTTACCGAGATTTTGAAgcgc                       >  2:336498/1‑226 (MQ=11)                                                                                                                                                                                                               ttACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCACGTCTTATCAGGCCTACGCCCGACATCGCAATAACCTGATTTAGCGTGATTTTTTAGGCCGGATCAgggcttt          >  2:533392/1‑223 (MQ=11)                                                                                                                                                                                                                tACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAt                                                                                                                                    <  1:442170/105‑1 (MQ=40)                                                                                                                                                                                                                tACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGTTTTAACTCTATTTTAGAGGACGGATAAGGCGTTTac        >  2:248102/1‑228 (MQ=11)                                                                                                                                                                                                                tACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                      <  2:699594/167‑1 (MQ=35)                                                                                                                                                                                                                tACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATAAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                      >  1:699595/1‑167 (MQ=25)                                                                                                                                                                                                                tACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAAATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCTTCCCCCACACAGCGCAATAGCCTGATTTCACGTGATTTTGTAGATCGGAAAAGGCGTTTac        >  2:466047/1‑228 (MQ=11)                                                                                                                                                                                                                      gcATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGTTAAGGCGTTTACGCCGCa  >  1:237662/1‑229 (MQ=21)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                        ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAC‑GCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCA  >  NC_000913/4295833‑4296274

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: ‑