Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,541,357	C→T	E412K (GAA→AAA)	narZ ←	nitrate reductase Z subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,541,357	0	C	T	97.2%	122.8 / ‑6.6	36	E412K (GAA→AAA)	narZ	nitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (17/18);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGATGCGCGGATCATGCGCCCGGGAACGTAGCTACCATCGTCGCGAGGCTCCAGCATTACCAGCATCGGCATGTCGCTGTAG  >  NC_000913/1541114‑1541592                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            cACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACcg                                                                                                                                                                                                                                      <  1:339011/251‑1 (MQ=255)          cATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCCATg                                                                                                                                                                                                                             >  1:88779/1‑251 (MQ=255)             cTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTTTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAACTGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATggct                                                                                                                                                                                                                          >  1:101062/1‑249 (MQ=255)                 cAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCATTGgattcc                                                                                                                                                                                                                       >  2:157251/1‑250 (MQ=255)                   gCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGgattcc                                                                                                                                                                                                                       >  1:390761/1‑248 (MQ=255)                    cccGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGgattcc                                                                                                                                                                                                                       >  1:155791/1‑247 (MQ=255)                     ccGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCTCCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCACGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGgattccagatt                                                                                                                                                                                                                  <  1:145246/251‑1 (MQ=255)                        cgATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGTCAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGgattccagattcc                                                                                                                                                                                                                >  2:321103/1‑250 (MQ=255)                        cgATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGgattccagattccc                                                                                                                                                                                                               >  2:176666/1‑250 (MQ=255)                               ggCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCt                                                                                                                                                                                                        <  2:106613/251‑1 (MQ=255)                               ggCCCTAATTCCCCAGTCCCAATTCATAAACGCTGACCACCGCAAAGGTGTTGCCATCGGTTAACGTCAGATTTTTAGCCGGCAATTGGCGCACCAGCTCAGAATTATGTTTTACGTTGCAAAAATTCTGACTTTCTATACCACCAAAGTAAGCGAAGGCCTCGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACTATTCGGTGATCGTACAGGCGGCAATGGATTCCAGATTCCATTTGcc                                                                                                                                                                                                         <  2:343332/250‑1 (MQ=255)                                                        aTAAACGCTGACCACCGGACAGGTGTTGCCATCGACTACCGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGCCCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAAAGGATTACAGATTCCATTTTCCTTTTTCTCCCCAGCGGAAACCAATCg                                                                                                                                                                               >  2:244382/1‑251 (MQ=255)                                                            aCGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAAcc                                                                                                                                                                           <  2:90586/251‑1 (MQ=255)                                                                     accGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCAc                                                                                                                                                                  <  2:347889/251‑1 (MQ=255)                                                                             ggtattagcATCGTCTAACGTCAGGTTTTTAACTGGCAATTGGCGCTCCAGCACCTGATTGTTTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCTCGCCTGCAACAGCGTCATGTGGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTCCCGGCGGCAATGGATTCCAGATTCCATTTGCGTTTTTCTCCCCAGCGTAAACCAATCGAACCGTTCGGCACTACCAAtt                                                                                                                                                          <  1:257599/243‑1 (MQ=255)                                                                                                gTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGc                                                                                                                                       >  1:343220/1‑251 (MQ=255)                                                                                                     gTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTGTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTGCGGCACTACCAATTCACCTGCGGTATTAACTGCTACgg                                                                                                                                  <  2:101062/251‑1 (MQ=255)                                                                                                       ttttAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTTCTACTGtt                                                                                                                                >  2:109103/1‑251 (MQ=255)                                                                                                        tttAACGGGCAATTGGCGCACCAGCGCCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCTTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGtt                                                                                                                                <  1:147689/250‑1 (MQ=255)                                                                                                         ttAACTGGCAATTGGCGCTCCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGttt                                                                                                                               <  2:203037/250‑1 (MQ=255)                                                                                                                   aTTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCgg                                                                                                                     >  1:192183/1‑250 (MQ=255)                                                                                                                                       ttGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGTGGAAGGCCACGCCTGCAACAGCGTCATGTTGGCCGAGCATGGTTAACGACAATTCGGTTTTCGTACCGGCGGCCATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTTCCAATTCACCTCCGGTATTAACTGCTACGGTTTTCCCCTGCGGATTGTTTCTTTCGCCCAGtcc                                                                                                >  2:19185/1‑251 (MQ=255)                                                                                                                                          tgtTTTACACTGCGAAAATGCGGATTTTCAATGCCGCCCAACCAGCGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGGGCAGGGTTAACGACCATTCGGATTTCGTACCGGCGGCAAAGGATTCCAGATTCCAGTTTCCTTTTTCTCCCCAGGGGAAACCAATCGAACCGCTCGGCCCTACCAATTCACCCGCGGTATTAACTTCTACGGTTTTTCCCCGCGGATTGTTGCTTTTGCCCAGtcctac                                                                                             >  2:67500/1‑248 (MQ=255)                                                                                                                                                       gAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGAt                                                                                <  1:70198/251‑1 (MQ=255)                                                                                                                                                            agCGGATTTACAATGCGTCCAAAGTAGTGGACGCCCACGCCTGCAACAGCCTCATGTTGACCGAGTAGGGTTAACGTCAATTCGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAGCCAATCGAACCGTTCGGCACTACCAACTCACCTGCGGGATTAACTGCGCCGGTTTTCCACTGCGGATTGTTGCTTTCCACCAGTCCATCCACCACGTCAGATGcgcg                                                                           <  2:60147/250‑1 (MQ=255)                                                                                                                                                                     tCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGATGCGCGGATCATgc                                                                   <  1:41925/250‑1 (MQ=255)                                                                                                                                                                                         gAAGGCCACGACTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGATGCGCGGATCATGCGCCCGGGAACGTAGCTACCAt                                              <  1:325022/251‑1 (MQ=255)                                                                                                                                                                                          aaGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGATGCGCGGATCATGCGCCCGGGAACGTAGCTACCAtc                                             <  1:319472/251‑1 (MQ=255)                                                                                                                                                                                          aaGGACACGCCGGCAACAGCGTCATGTTGACCGAGCAGGGGTAACGACAAGGCGGTTTTCGTACCGGCGGCAATGTATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTg                                                                                                                                        <  1:380408/160‑1 (MQ=255)                                                                                                                                                                                                    cTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGATGCGCGGATCATGCGCCCGGGAACGTAGCTACCATCGTCGCGAGGc                                   >  2:91714/1‑251 (MQ=255)                                                                                                                                                                                                       cAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGATGCGCGGATCATGCGCCCGGGAACGTAGCTACCATCGTCGCGAGGCTcc                                >  1:392743/1‑251 (MQ=255)                                                                                                                                                                                                       cAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAAAGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGATGCGCGGATCATGCGCCCGGGAACGTAGCTACCATCGTCGCGAGGCTcc                                >  1:394393/1‑251 (MQ=255)                                                                                                                                                                                                                       aCCGAGCAGGGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGATGCGCGGATCATGCGCCCGGGAACGTAGCTACCATCGTCGCGAGGCTCCAGCATTACCAGCATCg                >  2:326463/1‑251 (MQ=255)                                                                                                                                                                                                                         cGAGCAGAGTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCCAGTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGATGCGCGGATCATGCGCCCGGGAACGGAGCTCCAATCGTCGCGAGGCTCCAGCATTACCAGCATCGGc              <  1:326464/251‑1 (MQ=255)                                                                                                                                                                                                                               gggTTAACGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGATGCGCGGATCATGCGCCCGGGAACGTAGCTACCATCGTCGCGAGGCTCCAGCATTACCAGCGTCGGCATGTCg        <  2:200150/251‑1 (MQ=255)                                                                                                                                                                                                                                     aCGACAATTCGGTTTTCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGATGCGCGGATCATGCGCCCGGGAACGTAGCTACCATCGTCGCGAGGCTCCAGCATTACCAGCATCGGCATGTCGCTGTAg  <  2:392743/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            CACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAATGCCGCCAAAGTAGGGGAAGGCCACGCCTGCAACAGCGTCATGTTGACCGAGCAGGGTTAACGACAATTCGGTTTCCGTACCGGCGGCAATGGATTCCAGATTCCATTTGCCTTTTTCTCCCCAGCGGAAACCAATCGAACCGTTCGGCACTACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGATGCGCGGATCATGCGCCCGGGAACGTAGCTACCATCGTCGCGAGGCTCCAGCATTACCAGCATCGGCATGTCGCTGTAG  >  NC_000913/1541114‑1541592

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