Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	447,538	A→G	F56F (TTT→TTC)	cyoE ←	heme O synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	447,538	0	A	G	97.3%	121.3 / ‑4.1	37	F56F (TTT→TTC)	cyoE	heme O synthase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (23/13);  total (24/13)
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.29e-01

ATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAACC  >  NC_000913/447293‑447774                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              aTCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAacac                                                                                                                                                                                                                                          <  1:236503/251‑1 (MQ=255)           cGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGAc                                                                                                                                                                                                                                <  1:252716/251‑1 (MQ=255)            gTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACg                                                                                                                                                                                                                               <  1:107774/251‑1 (MQ=255)                   gagTGGCGTTTCATGTAAAAGCTATAAACCCCGCCATAAACCCAAAAGCCCCTCCCCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AAc                                                                                                                                                                                                                       <  2:343892/251‑1 (MQ=255)                            ttCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGacac                                                                                                                                                                                                              >  1:365068/1‑251 (MQ=255)                                             aaCGCCGACATAACCCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGACAGCAATACCCAGGAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTGCGTCCTTCCCATCTTTCTGACGATAATCCTGGCGATGTAGTTGTTGAACTCACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTata                                                                                                                                                                                             <  1:33684/251‑3 (MQ=255)                                                   gACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAAc                                                                                                                                                                                        >  1:297396/1‑250 (MQ=255)                                                         acccACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGATACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGAt                                                                                                                                                                                 <  1:36670/249‑1 (MQ=255)                                                            caaaaaGACCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACAACAGCGAGACAGACGGAGAGATCAGGCCATTCACCAGCACCCGATGCTTCGTCCTTACCATCTTTCTGTCGATATCCCCGTCACTGTAGTTGTTGAATACACAAACCGACGC‑ACTAACAAGTGACATCACAACCCACGTGTAGATAAACAAGGaataat                                                                                                                                                                              >  2:396417/4‑251 (MQ=255)                                                             ccaAAACCCATCCCCCCCAGCCAGCAGTCCAGCGGATACGCGCCAAACCACAGCAGCATAAAGCCAGCAATGCCCAGCAACGTGGCGTACCCCAGCGATACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTTCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTATTTGGTGAACACACAACCAGAAGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATc                                                                                                                                                                             <  2:387672/250‑1 (MQ=255)                                                                aaGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGACCACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCCGCGTGTAGATAAACAGGGGATAATCAAt                                                                                                                                                                          >  2:223318/1‑251 (MQ=255)                                                                 aGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAAtg                                                                                                                                                                         >  1:126938/1‑251 (MQ=255)                                                                   cccATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAAtgct                                                                                                                                                                       <  1:223318/251‑1 (MQ=255)                                                                              cagccagcAGGCAAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGAACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGCACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGGTAAACCGGGGAGAATCAAtgctgcacatcgagg                                                                                                                                                            >  2:376645/1‑242 (MQ=255)                                                                                           aGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATcc                                                                                                                                               <  1:66889/251‑1 (MQ=255)                                                                                                          aaCCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCg                                                                                                                                    >  1:407984/1‑247 (MQ=255)                                                                                                              acaGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATAAGGCCTTTCACCAGCAACCGAATCTTAGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCAC‑AACAAGTGACACCCCAACCAGCGTGTAGATAACCAGGGGATAATCAATGATGCCCATCGAGGCCAGCAGGAATCCACAAATCACCGAGATCAgg                                                                                                                            >  1:298810/1‑251 (MQ=255)                                                                                                                       caaaGCCAGCAATACCCAGCAACGTGGCATATCCCAGCGAGACAGCAGGAGAGATAAGGCCTTTCACCAGCACACGATGCTTCGTCCCTTCCATCTTTCCGTCGATATCCCTGTCGATGTAGTTGTAGAACACACAACCCGACGTCAC‑AACCAGTGACACCCCAGCGAGCGTTTAGATAAACAGGGGATAATCAATGCTGACCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAg                                                                                                                   <  2:135758/250‑1 (MQ=255)                                                                                                                        aaaGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGCCGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGGGATCAGGTTGCCaacg                                                                                                                   >  1:68491/1‑248 (MQ=255)                                                                                                                                 aaTACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTAAAACCGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTAGATATCCCTGTCGATGTAGTTGTTGAACACACAACGCGACGCCAC‑AACCAGTGACACACCAACCAGCGCGTAGATAAACAAGGGCTAATCAATGCAGCCCTTTGAGGCCAGAAGGAATACACCAATCACCGAGATCAGGATGCAAAGGATGATGCCTg                                                                                                         >  1:199536/1‑251 (MQ=255)                                                                                                                                    aCCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCCCCGAAGCCAGGTTGCCAAAGAAGATGccgagtt                                                                                                      >  2:432374/1‑246 (MQ=255)                                                                                                                                           aaCGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTAc                                                                                               >  1:61812/1‑251 (MQ=255)                                                                                                                                              gaggCGTACACCAGCGAGACAGCAGGAGAGCTCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACCCACAACCCGACGCCAC‑AACCAGTGACACCACAACCAGCATGTAGATAAACACGGGATAATCAATGCAGCCCTTTGCAGCCAGCAGGAATCCCCCAATCACCGAGATCCGGTTGCCAAAGAAGATGACTGATTTAGTTAcatg                                                                                            >  2:46501/3‑248 (MQ=255)                                                                                                                                              gTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTg                                                                                            >  1:389199/1‑251 (MQ=255)                                                                                                                                                     acacCAGCGAGCCAGCAGGAGATATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGTGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTa                                                                                      <  1:408860/250‑1 (MQ=255)                                                                                                                                                            agagacatgagGAGAGATCAGGCCTTTCACCAGCTCCCGATTCTTCGTCCTTTCCATCGTTCGGTCGATATCCCTGACGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTaa                                                                              <  2:147337/242‑1 (MQ=255)                                                                                                                                                              agaCAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAACGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAc                                                                            >  1:119983/1‑251 (MQ=255)                                                                                                                                                                    cagGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGCTATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAAcatcat                                                                       >  1:65108/1‑250 (MQ=255)                                                                                                                                                                      ggAGAGATCAGGCCTTTCACCAGCTCCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAc                                                                    <  1:413813/251‑1 (MQ=255)                                                                                                                                                                         gagaTCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAAccgc                                                                 >  1:352551/1‑251 (MQ=255)                                                                                                                                                                           gaTCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATATAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGCCACCCCAACCAGCGTGTAAATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGCGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATAATACccgccc                                                               >  2:281456/1‑250 (MQ=255)                                                                                                                                                                               aGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGctct                                                           >  1:396998/1‑251 (MQ=255)                                                                                                                                                                                          cAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGcatcat                                                >  1:92274/1‑251 (MQ=255)                                                                                                                                                                                                    ttcttcGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCAtgttgt                                        >  1:243915/1‑249 (MQ=255)                                                                                                                                                                                                                          tGTCGATATCCCTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGCTGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATcc                   >  2:133387/1‑248 (MQ=255)                                                                                                                                                                                                                                     cTGTCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTAc      >  1:281634/1‑250 (MQ=255)                                                                                                                                                                                                                                        tCGATGTAGTTGTTGAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAAcc  >  1:408034/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              ATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATAAAGCCAGCAATACCCAGCAACGTGGCGTACACCAGCGAGACAGCAGGAGAGATCAGGCCTTTCACCAGCACCCGATTCTTCGTCCTTTCCATCTTTCTGTCGATATCCCTGTCGATGTAGTTGTTAAACACACAACCCGACGCCAC‑AACCAGTGACACCCCAACCAGCGTGTAGATAAACAGGGGATAATCAATGCTGCCCTTTGAGGCCAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTTAAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAACC  >  NC_000913/447293‑447774

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