Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,161,142	0	C	A	23.7%	83.8 / 11.2	38	Q7K (CAA→AAA)	fabR	transcriptional repressor of fabA and fabB
Reads supporting (aligned to +/- strand):  ref base C (25/4);  new base A (0/9);  total (25/13)
Fisher's exact test for biased strand distribution p-value = 4.39e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAG  >  NC_000913/4160944‑4161380                                                                                                                                                                                                        |                                                                                                                                                                                                                                               ttATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAtt                                                                                                                                                                                                        >  1:644681/1‑239 (MQ=255)     cgaatatcctctCCGTTTAAGCGACATTTTTCCGTTGTACGTAAAAGAACCGGCAAAAGAATTGCAGTAAATATGTTTTATTGCTTTACCGTTCATTCACAATACTGGAGCAATCCATTATG‑TTCATTCTCGGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc                                                                                                                                                                                                    <  2:644681/231‑1 (MQ=255)                                          aaCGAAAAAGACCCGCCAAAAGATTTGCAGTAAATTGTTTTTATTGCTTTCCCTTTCATTCACAATGCTGGAGCAATCCAGTATG‑TTCATCCTCTGGTATAGTGCCAGCAGTCTTTTTGGAAAGGATTCAGACTTCGGGAGGGGGGTAAGAGCGAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGc                                                                                                                                                                                                    <  2:636268/201‑1 (MQ=255)                                                    aaCCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTg                                                                                                                                                                                          >  1:636268/1‑202 (MQ=255)                                                    aaCCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGt                                                                                                                                                     >  1:297079/1‑239 (MQ=255)                                                                                  ttATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCa                                                                                                                           >  1:540005/1‑235 (MQ=255)                                                                                          ttACCGTTCATTCACAATACTGGAGCAATCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAtt                                                                                                                                                                                               >  2:313820/1‑159 (MQ=255)                                                                                          ttACCGTTCATTCACAATACTGGAGCAATCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAAtt                                                                                                                                                                                               <  1:313820/159‑1 (MQ=255)                                                                                          ttACCGTTCATTCACAATACTGGAGCAATCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGGGCGGGAATCGGGCATTGCTCCCACCCCGTTTTTTc                                                                                                               >  2:20651/1‑239 (MQ=255)                                                                                                             acaggacAATCCCGTTTGTTTCTTTTTCTGTTAAAGTACCGGAAGTTCTTTTGAAAAGGATTCAGACATCGGGATGGGCAAAAGAGCACAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGc                                                                                                                                         <  2:270320/189‑1 (MQ=255)                                                                                                                 aGCAATCCATTATG‑TTCATTCTCTGGTATGGTCCCAGCAGTACTTTGGGAAAGGATTCAGACATCGTGAGGGGCGTAAGAGCAAAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGa                                                                                            <  1:438528/235‑1 (MQ=255)                                                                                                                     aTCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGaa                                                                                                                                                                                     >  2:71322/1‑142 (MQ=255)                                                                                                                     aTCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGaa                                                                                                                                                                                     <  1:71322/142‑1 (MQ=255)                                                                                                                     aTCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTTAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCCTTTTCGCGACGTATACGGACCgg                                                                                    >  2:369044/1‑239 (MQ=255)                                                                                                                     aTCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTgg                                                                                    >  1:338839/1‑239 (MQ=255)                                                                                                                      tCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTgat                                                                                                                               >  1:270320/1‑193 (MQ=255)                                                                                                                       acattATT‑TTCATTCTGTGATATTGCGCCAGTATTACTTTTGGCAATGATTCCGTCATCGTGGTGGGCGTAGGGGCCCAAAAAAAAGAAACACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGt                                                                                  <  2:540005/235‑1 (MQ=255)                                                                                                                           tATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTATCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTggg                                                                                   >  2:101214/1‑234 (MQ=255)                                                                                                                           tATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTggg                                                                                   >  1:244754/1‑234 (MQ=255)                                                                                                                            ttG‑TTCATCCTCGGGTATAGTCCAAGCAGTACTTTGGGCAAGGATTCAGACATCGTGATGGGCGAAAGAGCGCAAAAAAAAAAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGAc                                                                             <  1:20651/238‑1 (MQ=255)                                                                                                                                   aTTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggt                                                                                                                                                                                                                    >  2:657339/1‑98 (MQ=255)                                                                                                                                   aTTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCtggt                                                                                                                                                                                                                    <  1:657339/98‑1 (MQ=255)                                                                                                                                         tGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtgatga                                                                 >  1:662541/1‑239 (MQ=255)                                                                                                                                                 tGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGtt                                                         >  2:553498/1‑239 (MQ=255)                                                                                                                                                 tGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAAATGGGCTTGGCCATGGTTGATGGGACCGGtt                                                         >  2:641676/1‑239 (MQ=255)                                                                                                                                                     agcagTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAAt                                                     >  1:206504/1‑239 (MQ=255)                                                                                                                                                             tgtgaacAAGGACACAGACTGCGTGAGGGGAGAAAGCTCAAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGcc                                             <  2:206504/233‑1 (MQ=255)                                                                                                                                                             tattgcaaaGGATTCACACATCGTGATGGACGAAAGGCCAAAAAAAAAAGAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGcc                                             <  2:338839/232‑1 (MQ=255)                                                                                                                                                                tGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCa                                            >  1:129407/1‑237 (MQ=255)                                                                                                                                                                       gATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATgc                                    >  1:41041/1‑238 (MQ=255)                                                                                                                                                                                 tCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtgat                                                                   >  1:7597/1‑197 (MQ=255)                                                                                                                                                                                 tCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGt                         >  1:361237/1‑239 (MQ=255)                                                                                                                                                                                  cGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACctc                                                                                                                       >  2:171560/1‑144 (MQ=255)                                                                                                                                                                                  cGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACctc                                                                                                                       <  1:171560/144‑1 (MQ=255)                                                                                                                                                                                  cGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTc                        >  1:462491/1‑239 (MQ=255)                                                                                                                                                                                     gATGGGCGAAAGCGCCAAAAAAAAAGAAAACCCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGc                     <  1:369044/239‑1 (MQ=255)                                                                                                                                                                                          gCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATc                >  1:442435/1‑239 (MQ=255)                                                                                                                                                                                                   cgcAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAg                                                                                                                                                                                                               >  1:526591/1‑39 (MQ=255)                                                                                                                                                                                                     caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt                                                                                                                                                                                                         <  1:360310/43‑1 (MQ=255)                                                                                                                                                                                                     caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCAt                                                                                                                                                                                                         >  2:360310/1‑43 (MQ=255)                                                                                                                                                                                                     caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcg                                                                                                     <  1:582254/143‑1 (MQ=255)                                                                                                                                                                                                     caacaaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcg                                                                                                     >  2:582254/1‑143 (MQ=255)                                                                                                                                                                                                        caaAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAg  >  1:488366/1‑239 (MQ=255)                                                                                                                                                                                                        |                                                                                                                                                                                                                                               TTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATG‑TTCATTCTCTGGTATAGTGCCAGCAGTACTTTTGGCAAGGATTCAGACATCGTGATGGGCGTAAGAGCGCAACAAAAAGAAAAAACCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAG  >  NC_000913/4160944‑4161380

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