Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,446,482	C→T	L791L (CTG→TTG)	tamB →	translocation and assembly module subunit TamB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,446,482	0	C	T	93.5%	125.5 / ‑0.0	47	L791L (CTG→TTG)	tamB	translocation and assembly module subunit TamB
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (23/20);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

GCGGAACTGTGCGTGCCGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGCTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGTGAGTGCCAAC  >  NC_000913/4446242‑4446718                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             gCGGAACTGGTCGTGCCGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTGCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTGCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGc                                                                                                                                                                                                                                    <  2:32854/251‑1 (MQ=255)        tgtgCGTGCCGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGAGGCGTTGCCCGTGGCGTTTc                                                                                                                                                                                                                               >  2:6580/1‑249 (MQ=255)                 cGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGCATCTTCACGGGTAAAGCGGACGTTGTCCGGGACACAACGCAAGAGGGGCTGACGCAGGGCCGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTAAACGACGCGGCGTCGCCCGTGGCGTTTCAGACACTgcat                                                                                                                                                                                                                    >  1:15322/1‑248 (MQ=255)                     aaCTATCGATGCGGGGGCCGAAGGGCGTGGGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAGAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGa                                                                                                                                                                                                                <  1:419538/251‑1 (MQ=255)                       cTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGAcg                                                                                                                                                                                                              >  1:177798/1‑251 (MQ=255)                           cgaggCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGAcggcg                                                                                                                                                                                                           <  2:476513/246‑1 (MQ=255)                             aTGCGGGGGCCGAAGGGCGTGCGGTGGTGAAACTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGaa                                                                                                                                                                                                        <  1:133648/251‑1 (MQ=255)                                     gCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGACAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGCCGGAATTGCGTa                                                                                                                                                                                                 >  1:473041/1‑250 (MQ=255)                                         aaGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTTCATaacaac                                                                                                                                                                                            >  2:138291/1‑251 (MQ=255)                                          aGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGACTTGCGTAACAAcc                                                                                                                                                                                           >  1:454268/1‑251 (MQ=255)                                                                   cGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCCACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGCATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATccc                                                                                                                                                                  >  2:95079/1‑250 (MQ=255)                                                                               gctctGCTGAAACCGTTTATGACTGCAACCCCTCAGGCCAGCGGTAGCTTCACGGTTAAAGCGTACGTTGCCTGGGACACCACGAAAGAGGTGCTGCCGCAGGGAAGTATCACCCTTTCGGTGCGT‑AACGTGCAGGTAACGCAAACAGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGTGGAATTGAGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATaa                                                                                                                                                      <  2:378153/247‑1 (MQ=255)                                                                               gCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGCCACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCACACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGAAAATGAATCTGACGGCGGCATTGCGTAACAAACGTGCCGAATTGGGCTGGCCCATCCGCCTGACCAATa                                                                                                                                                       >  1:235243/1‑250 (MQ=255)                                                                                ccATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGAGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGCTGCCCGTGGCGTTTCAGCCACTGAATCTGACGTCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAAc                                                                                                                                                     >  2:114908/1‑251 (MQ=255)                                                                                         aacctgttATGCCAGAAACGACTCAGGCCAGCGGTATCTTCACGAGTAAAGCGGGGGTGGCCTGGGACACCAAGAAAGAGGGGCGGCCGCAGGTCAGTATCACCCGTTCGGGG‑GTGAACGTGCAGGTAACGCAAACCGGCAACGATGCGGCGGTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAAGTGCGTAACAACCGTGCCGAATGGGGTTGGACCATCCGCCTGCCCAATAACGGCCAGttt                                                                                                                                            <  2:166310/245‑1 (MQ=255)                                                                                                    ccAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGAcaggg                                                                                                                                 >  1:74289/1‑250 (MQ=255)                                                                                                     cAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGAcaggtg                                                                                                                                <  1:144054/251‑1 (MQ=255)                                                                                                        aaaCCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGAcaggtgcag                                                                                                                             <  1:3223/251‑1 (MQ=255)                                                                                                        aaaCCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTGGCCCGTGGCGTTTCAGACACTGAATCTGCCGGCGGACTTGCGTAACAACCGTGCCGAATTGGGCGGGACCCTCCGCCTGACCAATAACGGCCAGTTTGATGGAaaggtgc                                                                                                                               >  2:40959/1‑249 (MQ=255)                                                                                                           ccACTCAGGCCAGCGGTATCGTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGAcaggtgcaggtg                                                                                                                          <  2:323538/251‑1 (MQ=255)                                                                                                                cAGGCCAGCGGTATCTGCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCGGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGa                                                                                                                     <  1:205252/251‑1 (MQ=255)                                                                                                                 aGGCCTGCGGTATCTTCACGGGTAAAGCGGTCGTTGCCTGGGAAACCACGAAAGAGGGGCGGCCGCAGGGCAGTATCACCCTGTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGGCAACGATGCGGCGTTGCCCGTGGCGTTTCATACACTGAAGCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGAt                                                                                                                    <  2:104660/251‑1 (MQ=255)                                                                                                                         ggTATCTTCACGGGTAAAGCGGACGTTGCATGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACAGATCCGCAAgg                                                                                                            >  1:103590/1‑251 (MQ=255)                                                                                                                                          aGCGGACGTTGCCTGGGACACCACGAACGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGAAAACCGGCAACGATGCGGCGTTGCCCGTGGCGGTTCAGACACTGAATCTGACGGCGGACTTGCGTACCAACCCTGCCGAATTGTGCTTGAACATCCGCCTGACCAATACCGGCCCGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCgtggtg                                                                                           >  2:157027/1‑251 (MQ=255)                                                                                                                                          aGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGCATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTtggtg                                                                                           >  1:255563/1‑251 (MQ=255)                                                                                                                                                cGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGACGCCGTAATCTTGGTGGAAACg                                                                                     >  2:371643/1‑251 (MQ=255)                                                                                                                                                 gttgcgggggACCCCACGAAAGAGGGGCTGCCGCAGGGCAGTATCTCCCTTTCGGGGCTT‑AACGTGCGGGTAACGCAAACCGCAAACGATGGGGCGTTGCCCGTGGCGTTTCAGACCCTGAATCGGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAAGAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGt                                                                                    <  2:454268/244‑1 (MQ=255)                                                                                                                                                    gCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATg                                                                                                                                        <  1:46796/196‑1 (MQ=255)                                                                                                                                                           acacCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGAGCAGGTCACGCAAACCGTCAACGATGCGGCGTTGCCCGCGGCGCTTCAGACACTGAATCTGACGGCGGCATTGCGTAACAACCGTGCCGCATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGATTAATGGACAGGTGCAGGTGACCGATCAGCAAGGACGCCGTAATCTTGGTGGCAACGACAattgccgt                                                                          >  1:216700/1‑245 (MQ=255)                                                                                                                                                                                  cAGGGCAGTATCACCCTTTCGGGGCGT‑ACCGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATa                                                     <  1:188089/249‑1 (MQ=255)                                                                                                                                                                                  cAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGGTCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGCTaa                                                    >  2:188089/1‑250 (MQ=255)                                                                                                                                                                                    gggCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGTAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAAcc                                                 <  1:413446/251‑1 (MQ=255)                                                                                                                                                                                    gggCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGGGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGCTTGATGGACAGGTGCAGGTGACCGGTCCGCAAGGCCCCCGTAATCTTGGTGGCAACGCCAATAACCGTAACTTCAACCATTCGATGATaagcc                                                 >  2:117186/1‑248 (MQ=255)                                                                                                                                                                                      gCAGTATCACCCTTCCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAAcccc                                               <  1:6580/251‑1 (MQ=255)                                                                                                                                                                                          tATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCt                                           >  1:395911/1‑251 (MQ=255)                                                                                                                                                                                              ccccTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGGGTTTCAGACACGGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATGGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTAc                                       <  1:314248/250‑1 (MQ=255)                                                                                                                                                                                                       gggCGT‑AACGTGCAGGTAACGCAAACCGCCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGa                              <  1:65085/251‑1 (MQ=255)                                                                                                                                                                                                        ggCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGaa                             <  1:138291/251‑1 (MQ=255)                                                                                                                                                                                                         gCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGaa                             >  1:46932/1‑250 (MQ=255)                                                                                                                                                                                                          cGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGaaa                            >  1:187176/1‑250 (MQ=255)                                                                                                                                                                                                                 gTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGCCGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACGATCCGCCt                                                                                                                                                               >  2:188810/1‑113 (MQ=255)                                                                                                                                                                                                                 gTGCAGGTAACGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCt                                                                                                                                                               <  1:188810/113‑1 (MQ=255)                                                                                                                                                                                                                   gCAGGTAACGCAAACCGTCAACGATGCGGCGCTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAAATGGGCTGGACCATCCGCCTGACCAATAACGGCCAGCTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCAGAAATCTTGGTGGCAACGACAATATCCGAAACTTAAACCTTGCGATGATAAACACCAGCTTTACCAGTGGGGAACAagactagg                   >  2:10284/1‑245 (MQ=255)                                                                                                                                                                                                                      ggTAAAGCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTCATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGAgggg                 >  1:176206/1‑250 (MQ=255)                                                                                                                                                                                                                         aaCGCAAACCGTCAACGATGCGGCGCTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATg              <  1:248918/250‑1 (MQ=255)                                                                                                                                                                                                                            gCAAACCGTCAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGGGCCGGGGACCGTGCCGCAAGGCCGCCGTAATCTTGGGGGCAAAGTCAATATCCGTACCTTAAACCGTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGTgg          >  2:145912/1‑250 (MQ=255)                                                                                                                                                                                                                                     cAACGATGCGGCGTTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGTGAGTGCCAAc  <  1:109238/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             GCGGAACTGTGCGTGCCGCAAACTATCGATGCGGGGGCCGAAGGGCGTGCGGTGGTGAATCTCAACCGCTTCGACCTCGCCATGCTGAAACCGTTTATGCCAGAAACCACTCAGGCCAGCGGTATCTTCACGGGTAAAGCGGACGTTGCCTGGGACACCACGAAAGAGGGGCTGCCGCAGGGCAGTATCACCCTTTCGGGGCGT‑AACGTGCAGGTAACGCAAACCGTCAACGATGCGGCGCTGCCCGTGGCGTTTCAGACACTGAATCTGACGGCGGAATTGCGTAACAACCGTGCCGAATTGGGCTGGACCATCCGCCTGACCAATAACGGCCAGTTTGATGGACAGGTGCAGGTGACCGATCCGCAAGGCCGCCGTAATCTTGGTGGCAACGTCAATATCCGTAACTTCAACCTTGCGATGATAAACCCCATCTTTACCCGTGGGGAAAAAGCAGCGGGGATGGTGAGTGCCAAC  >  NC_000913/4446242‑4446718

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