Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,310,770	A→G	H32H (CAT→CAC)	yciI ←	protein YciI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,310,770	0	A	G	89.2%	92.2 / 4.0	37	H32H (CAT→CAC)	yciI	protein YciI
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base G (18/15);  total (20/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.77e-01

AAAAACCGTGCCCTTAAATATTGAATCT‑CTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCC‑TGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCT  >  NC_000913/1310537‑1311013                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    aaaaaCCGTGCCCTTAAATATTGAATCT‑CTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCC‑TGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGTCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACg                                                                                                                                                                                                                                    <  1:30109/251‑1 (MQ=255)      ccGTGCCCTTAAATATTGAATCT‑CTATCCAGAACACTTTCTTAAATGGTTTCGCGGAACGGTGTTCATAGACTCC‑TGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTa                                                                                                                                                                                                                               <  2:173285/251‑1 (MQ=255)       cgttcccTTAAATATTGAATCT‑CTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTAC‑TGACGCTACGTACGGGTCAGCATCGTCCCAGGCCTGCGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTAAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGCAAACACACGCCCTTCATCGTACAGTAACAGTAACCGTGCTa                                                                                                                                                                                                                               >  2:142468/5‑250 (MQ=255)        gttctcTTAAATATTGA‑TCTTCTTTTCAGAACACTTTCTTAAATGGTTTCACTGAACCTTTGTATTAG‑CTCCCTGCCGCTACGTACGGGTCAGCATCTGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCTCAGCCCCTGGTGCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCAGCGTGCAGTAACTGTAATCGTTCTaaa                                                                                                                                                                                                                             <  2:446214/246‑1 (MQ=255)               tAAATATTGAATCT‑TTATTCAGAACACTTTCTTAAATGGTTTCCCTGAAACGTGTTCATAGCCTCC‑TGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGcc                                                                                                                                                                                                                       <  2:291953/250‑1 (MQ=255)                    atTGAATCT‑CTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCC‑TGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCggccg                                                                                                                                                                                                                  >  1:410092/1‑247 (MQ=255)                      tGAATCT‑CTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCC‑TGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACACGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACgcac                                                                                                                                                                                                               >  2:109726/1‑248 (MQ=255)                            cgctATTCAGAACACTTTCTTAAATGGTTTCATTGAAACGTGTTAATAGACTCC‑TGCCGCTACGTACAGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTACATCCCGCAGCCCCTGGATCTTTACTGTCTACCGCTGGCATTGGCCCAGCTGTCAACAAACGACCTTCATAGTGCAGTCACTGTAAGCTAGCTAAGTGTGCCGGACTAAAGGaa                                                                                                                                                                                                           >  2:186470/3‑250 (MQ=255)                                   cAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCC‑TGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATACCGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGATAAATGAGCCGGACGAACGGAACGGCGttt                                                                                                                                                                                                   >  2:260155/1‑251 (MQ=255)                                          cTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCC‑TGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATAGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGAAGCACGGAAGGGGGTTTTTCGAgg                                                                                                                                                                                            >  2:393324/1‑251 (MQ=255)                                           tttCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCC‑TGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGc                                                                                                                                                                                           >  1:257897/1‑251 (MQ=255)                                                              aaaacttttttaTAGACCCC‑TGCCGCAACGTCCGGGTCAGCATCGGCCCAGGCCTGAGCTGCGTCCAGCGAGTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCGGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGTCTATCAGCTTTGCCTTGAGc                                                                                                                                                                        <  2:461743/240‑1 (MQ=255)                                                                           aGACCCC‑TGCCCCTACGTACGGGTCAGCATCGGCCCAGGACTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATac                                                                                                                                                           <  2:116809/251‑1 (MQ=255)                                                                             ccccc‑TGCCGTTACTTACGGGGCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCATTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATacaa                                                                                                                                                         <  2:425199/248‑1 (MQ=255)                                                                                               cgGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCtt                                                                                                                                        >  1:296714/1‑251 (MQ=255)                                                                                                  gTCAGCACCGGCCCAGGCCTGAGCTGCTGCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGTGCACAACGCCTTGtt                                                                                                                                     <  2:98157/251‑1 (MQ=255)                                                                                                           ggCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGt                                                                                                                            <  1:428596/251‑1 (MQ=255)                                                                                                            gCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCATAAACGGGTGAGCCAGTAAATCCCGCAGCTCTTGTATCGTTACTGTTTACCGCTGGCATTGGCCCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTAGCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTGCGGTAAAGtt                                                                                                                           <  2:221221/251‑1 (MQ=255)                                                                                                                   gcgtaGCTGCTTCCAGCGATTCAAATTCAACAATAACGGTTGAGCCAGAAAAGCCCGCAGCACCTTGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCCTCATCGTGCAGGAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCCCAACTCCTTGTTCGTGAAAGTTGTAAAGt                                                                                                                    <  2:246120/247‑1 (MQ=255)                                                                                                                      gAGCTTCTTCCAGCGATTCAAATTCAGCAATAACGGGTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTGACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCGTCATCGTGCAGTAACTGGAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTGTGTCGTGAGCGTAAAGAACATACAACACGGGCACAACTCCTTGTTCTGGAAAGTGGTAAAGTACg                                                                                                                 <  2:214620/251‑1 (MQ=255)                                                                                                                               ccAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGa                                                                                                         >  1:369504/1‑250 (MQ=255)                                                                                                                                      ttCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTGCCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAACGTACGTTATTTGAAAGGGCg                                                                                                  >  2:232666/1‑250 (MQ=255)                                                                                                                                                gCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCa                                                                                        >  1:108648/1‑250 (MQ=255)                                                                                                                                                        ggTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGTAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGaa                                                                                  >  1:337328/1‑248 (MQ=255)                                                                                                                                                                      aTCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCAGTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTa                                                                 <  2:296714/251‑1 (MQ=255)                                                                                                                                                                              gCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATg                                                         <  1:260155/251‑1 (MQ=255)                                                                                                                                                                                   tGGATCGAAACTGTCTCCCGCTGGCATGGGACCAGCTGTCAACCAACGACCTTCATAGTGCAGTAACTGTACACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACCAATCCTTGCTCGGGAAAGTTGTAAAGTACGTTATTTGAAAAGGCGAATATCTGCAACGGAACGAAGATGTCTTTGTTAAGGCCATGAATa                                                     >  2:465459/1‑250 (MQ=255)                                                                                                                                                                                             cTGTCTACCGCTGGCATTGGACCAGATGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGCGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCAGAACGTAAGGGt                                           >  1:321527/1‑250 (MQ=255)                                                                                                                                                                                              tGTCTACCGCTGGCATTGGACCAGCTGTCAACACACGCCCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCAGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACCTACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGATGCTCTGCAACGGAAAGATGCTGTCTTTGTTAAGGCCATGCATAATATAAGGGTa                                          >  2:454938/1‑250 (MQ=255)                                                                                                                                                                                                         ggCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCaaaa                              <  1:232666/251‑1 (MQ=255)                                                                                                                                                                                                            aTTGGACCAGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGATTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGAACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTGAGGCCATTCTTAAAGAAAGGGTAATTACGCAAAAAATg                           >  1:400304/1‑251 (MQ=255)                                                                                                                                                                                                               ggACCCGCTGTCAACAAACGACCTTCATCGTGCAGTAACTGTAAACGGGCTAAATGTGCCGGCAGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATACCATACAGCACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCTTGCATACAGTAAGGGGAATTACGC‑AAAAAttaacc                       >  2:436537/1‑246 (MQ=255)                                                                                                                                                                                                                           aacaaACGACCTTCATCGTGCGGTAACTGTAAACGAGCTAAATAAGCCGGACGAACGGAAAGGAGTTTATCGAGGCTATCCGCTTTGTATTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTAGGGCAAGTAGTAACGTAAGACATTTGACAGGACGAAGATCGGCAAAGGAAAGCTGATGTCTTTGTTCATGCCATGCATAAAGTACGGGAAATTACACCAAAAATGAAATTTTAActaaac            >  2:397281/1‑247 (MQ=255)                                                                                                                                                                                                                           aacaaACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAAACTACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTTAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCActgatc            >  1:29280/1‑251 (MQ=255)                                                                                                                                                                                                                                 cGACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTTCCGGACGAACTGAAAGGCGTTTTTCGAGGCGATCAGCATTGTCTTGAGCGTAAATAACATACAACACGGGCTCAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCActgatcctgatc      <  2:419303/251‑1 (MQ=255)                                                                                                                                                                                                                                  gACCTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCg     <  1:402946/251‑1 (MQ=255)                                                                                                                                                                                                                                     cTTCATCGTGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGAACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCt  >  2:331297/1‑251 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                    AAAAACCGTGCCCTTAAATATTGAATCT‑CTATTCAGAACACTTTCTTAAATGGTTTCACTGAAACGTGTTCATAGACTCC‑TGCCGCTACGTACGGGTCAGCATCGGCCCAGGCCTGAGCTGCTTCCAGCGATTCAAATTCAGCAATAACGGTTGAGCCAGTAAATCCCGCAGCCCCTGGATCGTTACTGTCTACCGCTGGCATTGGACCAGCTGTCAACAAACGACCTTCATCATGCAGTAACTGTAAACGTGCTAAATGTGCCGGACGAACGGAAAGGCGTTTTTCGAGGCTATCAGCTTTGTCTTGAGCGTAAATAACATACAACACGGGCACAACTCCTTGTTCGGGAAAGTTGTAAAGTACGTTATTTGAAAGGGCGAAGATCTGCAACGGAAAGATGATGTCTTTGTTAAGGCCATGCATAAAGTAAGGGTAATTACGCCAAAAATGACATTTTCACTGATCCTGATCGTCT  >  NC_000913/1310537‑1311013

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