Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,423,404	0	G	T	54.0%	16.5 / 45.1	37	R2L (CGA→CTA)	rzoR	putative prophage outer membrane lipoprotein RzoR
Reads supporting (aligned to +/- strand):  ref base G (10/7);  new base T (9/11);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGGTGCTCAGCAATACATAACCGAGCAGTGTTTAAAGTAAAATC  >  NC_000913/1423162‑1423648                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     gggCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTg                                                                                                                                                                                                                                               >  1:29437/1‑250 (MQ=255)           ttACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGct                                                                                                                                                                                                                                    >  1:58212/1‑251 (MQ=255)                   atgagcagaggcaccgcgattATCTCCGCTCGGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGtg                                                                                                                                                                                                                             <  2:59478/230‑1 (MQ=11)                        cGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGatgatg                                                                                                                                                                                                                       <  1:191467/251‑1 (MQ=255)                         gAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGt                                                                                                                                                                                                                      >  1:66571/1‑251 (MQ=255)                         agagtcaccgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAACGCCCAGCACGACAAAAATGCCAGCGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGAGGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGt                                                                                                                                                                                                                      >  1:41155/15‑251 (MQ=11)                            gtcaccgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGc                                                                                                                                                                                                                   >  2:360944/12‑251 (MQ=11)                             gcagcaTCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGcc                                                                                                                                                                                                                  >  2:301420/1‑251 (MQ=255)                                      ttatccctgctatGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCCTTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCAAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGCTGCGTCAGCGTGATGTTGCTGCGCTCGATGCACAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAACAAGATGCTCTGCGTGATGATGTTGCAGCTGgtcgt                                                                                                                                                                                                         >  2:331404/13‑251 (MQ=11)                                            ccGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcgtcg                                                                                                                                                                                                       <  1:63047/247‑1 (MQ=11)                                               ctctGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCTGCAATTACTGACATTCAGATGCGTCAGCGTGATGTGGCTGCGCTCGATGCAAAAGACACGAATGAGTTAGCTGATGCTAAAGc                                                                                                                                                                                                                                                 <  2:209568/202‑1 (MQ=11)                                               ctctGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGc                                                                                                                                                                                                                                                 >  1:209568/1‑202 (MQ=11)                                                tctGGATATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGc                                                                                                                                                                                               >  1:244226/1‑251 (MQ=11)                                                   ggTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGc                                                                                                                                                                                               <  2:44273/248‑1 (MQ=11)                                                     ttATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAAACATCACCATGATAAAGCCATTACCTACAAAGACCAGCGCGACAAAACTGCCAGCGAACAGACACTGGCGCACGCGGCAATTACTGACATGCAGAAGCCTCAGCGTGATGCTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGACAATGATGCTCTGCGTGATGATGTTACCGCTGGTCGTCGTCGGTTTCACATc                                                                                                                                                                                          >  1:363321/1‑251 (MQ=11)                                                           gCATCGTCGTCTGCCTGTCAGGGGGTGTTAATCATTACCTTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATTCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCTTTATGTTGCTGCGCGCGATGCAAAAGACACGTAGGAGTTAGCTGATGCGAAGGCGGAAAATGATGCTCTTCGTGATGATGGTGCCGCTGTTCGGCGTCGTTTGCACATCAAAGCa                                                                                                                                                                                    <  2:356071/251‑1 (MQ=255)                                                             atcatcgtcTGCCTGGCATGGGCTGTTAATCATTACCGTGATAACGCCATTTCCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGGTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCGCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCGCTGCGTGATGAGGTTGCAGATGGTCGTCGTCGGTTGCACATCAAAGCAGt                                                                                                                                                                                  <  2:324304/247‑1 (MQ=11)                                                                      tGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTgtcagtc                                                                                                                                                                         >  2:260662/1‑251 (MQ=11)                                                                       gCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTTGTCGTCGTCGGTTGCACATCAAAGCAGTCTgtcagt                                                                                                                                                                          <  2:29437/249‑1 (MQ=255)                                                                       gCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTgtcagt                                                                                                                                                                          <  2:314519/249‑1 (MQ=255)                                                                          tGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTg                                                                                                                                                                     >  2:170888/1‑251 (MQ=255)                                                                           gTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTg                                                                                                                                                                     <  1:170888/250‑1 (MQ=255)                                                                           gTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGACCTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGc                                                                                                                                                                    >  2:157309/1‑251 (MQ=11)                                                                                   ctgtTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGcc                                                                                                                                                            >  1:208725/2‑251 (MQ=255)                                                                                                               taccTACAAAGCCCAGCGCGCCAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGAGGCTTCAGCGTGAGGTGGCTGCGCTCAATGCAAAAGACACGAAGGAGTTAGCTGAGGCTAAAGCTGAAAATGATGTTCTGCGTGATGATGTTGCAGCGGGTCGGCGTCGGTTGCCCAGCAAAGCAGTCGGACAGTCAGTGCGGGAAGCCACCGCCGCCTCCGTCTTGGATAATGCa                                                                                                                                 <  2:41155/248‑1 (MQ=255)                                                                                                                      aaagcccAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTcccc                                                                                                                         <  1:66993/245‑1 (MQ=2)                                                                                                                            cAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGTCATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACt                                                                                                                   <  2:192962/251‑1 (MQ=14)                                                                                                                                gcgACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCGCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCa                                                                                                               <  1:157309/251‑1 (MQ=14)                                                                                                                                                   aCTGAAGTTGGCGAACGCGGCAATTACTGGCATGAAGATGCGTCAGCTTGATTTTGCTGCGCTCGATGCAAAATACATGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCGTCCGGCGTGGATAATTCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGatt                                                                                            <  1:46966/251‑1 (MQ=11)                                                                                                                                                      gAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGGGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGattat                                                                                          <  2:100702/250‑1 (MQ=14)                                                                                                                                                                          ttACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTg                                                                     >  2:166792/1‑251 (MQ=38)                                                                                                                                                                                               aGCGTGATGTTGCTGCGCTCGATGCAAAATACACCAAGGAGTTAGCTGATGCGCGAGCTGAAAATGATGCTCTGCGTGCTGATGTTGCCGCTGGACGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATACTGCAGCCTCCCCCCGACTGGCAGACACCCCTGAACGGGATTATTTCACCCTCCGGGCACGACTGGTAATGATGCAGGCCCAACtt                                                >  2:162545/1‑251 (MQ=14)                                                                                                                                                                                                                        aaaaTACACGAAGGAGTTAGCGGATGCGAAAGCGGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGGGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGGTGCTCAGCAATACATAACCg                       <  2:58212/251‑1 (MQ=255)                                                                                                                                                                                                                                  aaGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGGTGCTCAGCAATACATAACCGAGCAGTGttt             >  1:7482/1‑251 (MQ=255)                                                                                                                                                                                                                                         tAGCTGATGCGAAAGATGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCCAAGCAGTCTGTCCGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGCTAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGGTGCTCAGCAATACATAACCGAGCAGTGTTTAAAGTa       >  2:201255/1‑250 (MQ=255)                                                                                                                                                                                                                                                tGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGGTGCTCAGCAATACATAACCGAGCAGTGTTTAAAGTAAAATc  >  1:260035/1‑248 (MQ=255)                                                                                                                                                                                                                                                ggCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGGCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGGTGCTCAGCAATACATAACCGAGCAGTGTTTAAAGTAAAATc  <  2:260035/247‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     GGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGACTGGTAATGATGCAGGCCCAACTTGAAGGTGCTCAGCAATACATAACCGAGCAGTGTTTAAAGTAAAATC  >  NC_000913/1423162‑1423648

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