Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,869,637	C→A	L87L (CTC→CTA)	cbrA →	colicin M resistance protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,869,637	0	C	A	94.4%	90.8 / ‑6.4	36	L87L (CTC→CTA)	cbrA	colicin M resistance protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (19/15);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

GTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCAT  >  NC_000913/3869392‑3869883                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       gTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGGGCTGGATAACAAGCACCAGTGTGGTACTGAAGGTTTAAGCAAACCCTGCGGCGGTCTGCTGGCACCGGGCGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGCTGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCACCGCATCGCTAAcgcg                                                                                                                                                                                                                                                   >  1:173950/1‑251 (MQ=255)  tGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGCTGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGGCGACGTCGCCGCATCGCTAACGCGt                                                                                                                                                                                                                                                  >  1:190651/1‑251 (MQ=255)   ggCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGGCGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGCTGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGCAGTCGCCGCATCGCTAACGCGt                                                                                                                                                                                                                                                  >  2:41028/1‑250 (MQ=255)                  cGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGc                                                                                                                                                                                                                                   <  1:143430/250‑1 (MQ=255)                   ggCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCta                                                                                                                                                                                                                                 <  1:35951/251‑1 (MQ=255)                        ggCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATc                                                                                                                                                                                                                             >  1:342522/1‑250 (MQ=255)                          cAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAat                                                                                                                                                                                                                          <  1:371833/251‑1 (MQ=255)                                   ggTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGc                                                                                                                                                                                                                 <  2:364563/250‑1 (MQ=255)                                           cGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATGCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCttt                                                                                                                                                                                                         <  2:258612/251‑1 (MQ=255)                                             aaaGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGAACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGAAGC‑TTCCTGTCGAAGTGATCGCTAATCCGCAGATTTTCAGCGGCAAAACCGTCGGCGGCGCCGCATCGCTAAAGCGAAAATACCAGCGAAGCAATAGCAATAATAATCGCAATGCTTTCg                                                                                                                                                                                                       >  2:109914/1‑251 (MQ=255)                                                                 aGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTACAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCCGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGGTGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATTTTAATCGCCATGCTTTAGACTTGAGGATGAAATCGCTg                                                                                                                                                                                   >  1:320837/1‑251 (MQ=255)                                                                    ggATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCGGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATc                                                                                                                                                                                <  2:109820/250‑1 (MQ=255)                                                                          cgcTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTGTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATAAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGcc                                                                                                                                                                          <  1:142654/251‑1 (MQ=255)                                                                           gcTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGcc                                                                                                                                                                          >  1:130435/1‑250 (MQ=255)                                                                                                    aCTGAAGGTTTCAGCAAACCCTGCTGCGGTCTGCTGGCACCGGCCGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCt                                                                                                                                                >  1:109425/1‑251 (MQ=255)                                                                                                          ggTTTCAGCAAACGCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCTATCTATTATTCGCGATGGGCTGAACC‑TTCCTGTCGATGAGCTCGCTAATCCGCAGATTTTCAGCGTCCAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGCATTTTGGATGAAATCGCTGATCCCCGCCAGCGCTGAGGTTTATCCCGACAGCCTGTGCCg                                                                                                                                          >  2:179457/1‑251 (MQ=255)                                                                                                                                   tgctgGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAGCCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGCTaa                                                                                                                  >  2:250740/1‑250 (MQ=255)                                                                                                                                      gggcccggggcgcgGAGCGATTTTTTATTCGCGGTGGGCTG‑CGCCTTCCGGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAACCCGTCGACGTCGCCGCATCTCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGCCTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATg                                                                                                               <  2:4002/240‑1 (MQ=255)                                                                                                                                            cGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATTTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGATTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTc                                                                                                        <  1:92745/251‑1 (MQ=255)                                                                                                                                                   ggaGCGATCTTGTATTCGGGAGGGGCTGACGC‑TTCCGGTCGATGGGTTCGCTAATCCGAAGATTTTCAGCGCCAAAACCGTCGGCGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATAGCAATATGAAGCGCCATGTTTTCGACTTGTGGATGAAATCGCGGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGAGGATAAATGGCATGTCATTTTTc                                                                                                 <  2:186735/249‑1 (MQ=255)                                                                                                                                                    cAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGTTTTTCAGCGTCAAAACCGTCGAAGTCGACGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATACTTTCGACTTGTGGATGAAATCGCTGCTACCCGCCAGAGTTGAGGTTTATCACGACAGCCTGTGCCGCAAGCTCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCg                                                                                                >  2:32479/1‑251 (MQ=255)                                                                                                                                                    cAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCACAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGCTCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTct                                                                                                >  2:384756/1‑250 (MQ=255)                                                                                                                                                       cGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCACTATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCCCGACa                                                                                                                                                    >  2:361836/1‑196 (MQ=255)                                                                                                                                                       cGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACa                                                                                                                                                    <  1:361836/196‑1 (MQ=255)                                                                                                                                                         aTCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGc                                                                                             >  2:51713/1‑249 (MQ=255)                                                                                                                                                                       aTGGGCTGACGC‑TTCCTGTCTATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTGCGACTTTTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTTCCGAAAAATCTg                                                                                                                                 <  2:185977/199‑1 (MQ=255)                                                                                                                                                                       aTGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCACAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTTCCGAAAAATCTg                                                                                                                                 >  1:185977/1‑199 (MQ=255)                                                                                                                                                                             tGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCAAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGCCGGCTGGGAGCAACatgtt                                                                       >  1:403622/1‑248 (MQ=255)                                                                                                                                                                             ctgccgcatCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCCAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATCTCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCTTTGAGGTTTATCACGACAGCCTGTGCCGACAAATCTGGCGTGCGGATGATAAATGGCATGTCATTTTTCGGGCAGACGGCTGGGAGCAACctat                                                                        >  2:387338/9‑251 (MQ=255)                                                                                                                                                                                 gC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccg                                                                   >  1:51366/1‑251 (MQ=255)                                                                                                                                                                                     cccggcggattggATCGCGATTCCGCAGTTTTCCAGCGACAAAACGTTCGACGGCCCCGCCCCGTTAACGCGTAACTACCAGCGAAGCTATATAAATATTAACCGCCATGCTTTCGATTTGGGGATGAACTCCCTGACCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTAccgcc                                                                 <  2:342522/238‑1 (MQ=255)                                                                                                                                                                                             aTGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAGCCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTg                                                        <  1:250740/251‑1 (MQ=255)                                                                                                                                                                                                        aaTCCGCAAATTGTCCGCGTCAAAACCGTCGCCGTCGCCGCATCGCTACCGCGTAACTACCAGCGAATCTATATCAATATTAACCGCCATGCTTGCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGcc                                               <  2:130435/249‑1 (MQ=255)                                                                                                                                                                                                            cGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGgc                                         <  1:384756/251‑1 (MQ=255)                                                                                                                                                                                                              cAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGgcgc                                       >  1:319984/1‑251 (MQ=255)                                                                                                                                                                                                                                                   cGCTAACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCAATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGAGGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCAGTTGGCGCAAACTCGAAGGTGAGGCGACATCTCTACCCGGgtcat  >  2:351474/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       GTGGCGATTATTGGCCTCGGCCCGGCAGGATCGGCGTTGGCACGAAAGTTAGCAGGCAAAATGCAGGTGATCGCGCTGGATAAAAAGCACCAGTGTGGTACTGAAGGTTTCAGCAAACCCTGCGGCGGTCTGCTGGCACCGGACGCGCAGCGATCTTTTATTCGCGATGGGCTGACGC‑TTCCTGTCGATGTGATCGCTAATCCGCAGATTTTCAGCGTCAAAACCGTCGACGTCGCCGCATCGCTCACGCGTAACTACCAGCGAAGCTATATCAATATTAATCGCCATGCTTTCGACTTGTGGATGAAATCGCTGATCCCCGCCAGCGTTGAGGTTTATCACGACAGCCTGTGCCGAAAAATCTGGCGTGAGGATGATAAATGGCATGTCATTTTTCGTGCAGACGGCTGGGAGCAACATATTACCGCCCGCTATCTGGTCGGTGCCGATGGCGCAAACTCGATGGTGCGGCGACATCTCTACCCGGATCAT  >  NC_000913/3869392‑3869883

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