Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,403,010	C→T	R325C (CGT→TGT)	hflK →	regulator of FtsH protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,403,010	0	C	T	98.5%	198.4 / ‑5.1	66	coding (973/1260 nt)	hflK	regulator of FtsH protease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/36);  minor base . (0/1);  total (29/37)
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.77e-01

CGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGC  >  NC_000913/4402765‑4403252                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    cgccgGAAGAAGCAAAAGCGACGTTTGACGAGTCGATTGCCGCGCGTGAAAACGATCGGTAATACATTCGTGAAGCAGAAGCGTATACGAACGAAGTTCAGCCGCGTGCGAACGTTCTGGCGTAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCTCGCTTTGCTAAACTTCTGCCTGAAGATAAAGCCGCGCCGGAACGTACGCGCGAGTGT‑CTg                                                                                                                                                                                                                                               <  2:326088/251‑1 (MQ=255)   ccgGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAATGCCACGCCGGAAATTACTAGCGAGTGT‑CTGt                                                                                                                                                                                                                                              >  2:183380/1‑250 (MQ=255)   ccgGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGta                                                                                                                                                                                                                                             <  1:98962/251‑1 (MQ=255)    cgGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGtat                                                                                                                                                                                                                                            <  1:110841/251‑1 (MQ=255)    cgGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGtat                                                                                                                                                                                                                                            <  2:284585/251‑1 (MQ=255)        agaagTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCg                                                                                                                                                                                                                                        <  1:547292/251‑1 (MQ=255)          aagTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCgag                                                                                                                                                                                                                                      <  1:330643/251‑1 (MQ=255)            gTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGAc                                                                                                                                                                                                                                    <  1:291454/251‑1 (MQ=255)            gTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGAc                                                                                                                                                                                                                                    >  1:286115/1‑251 (MQ=255)                    ggcttgtGACGATGCGATGGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGaaa                                                                                                                                                                                                                            <  1:323784/245‑1 (MQ=255)                              gATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGAGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGATGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGAGCCGGAAAATAATCGCGAGTGT‑CTGTATATCGAGACGATAGAAAAAATGTgggg                                                                                                                                                                                                                  >  2:437463/1‑251 (MQ=255)                                tGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGGTgggg                                                                                                                                                                                                                 >  2:162962/1‑249 (MQ=255)                                    aatGCCGGGCGTGAAACCGGACAGCAAGAGATGCGTGAAGCAGAAGCGGATACCGACGAAGGTCAGCCGCGTGCGACCGGTCAGGCGAAACGTATCCGCGGAGAGGCGCGGGCGTACAAGGCCCAGTCCATCCGGGAAGCAAAGGGGGAATGGGCGAGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGATATTTCTCGCGAGTGT‑CTGTATATCTCTCCGATGTAAATAGTGTTGGGTAacac                                                                                                                                                                                                            <  2:23702/249‑1 (MQ=255)                                       gCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTGCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACCCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACAccc                                                                                                                                                                                                          <  1:115817/250‑1 (MQ=255)                                            gcgTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAg                                                                                                                                                                                                    <  1:219455/251‑1 (MQ=255)                                              gtcaaaaCGGACAGAAATACATTCGTAAAGCAGAAGCGTATACCAACGAAGTTCAGCAGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCGCAGGGTGCAGGGGGGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGTGTGT‑CTGTCTTACACGACGACGGACAAAGAGTTGGGTAAAAACCGCAAAGTg                                                                                                                                                                                                  >  2:93102/4‑251 (MQ=255)                                                taaaaCGCAGAGCAATACCGTCGTGAAGCATAAGCGTATACCCACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGACAGCTCAGGGTGAATTGGCGCGATTTGATAAACTTCTGCCGGAATATAAAGCCGCGCCGGAACTTACTCGCGAGTTT‑CTGTATAGCAAGACGATGGAAAAAGAGTTGGGTATCACACGcccaggtct                                                                                                                                                                                                >  2:77585/2‑243 (MQ=255)                                                  aaaCGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTgg                                                                                                                                                                                              <  1:403166/251‑1 (MQ=255)                                                     cGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGATCAGGGCGAAGGGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGCGTGT‑CTTTATATCGAGACGAAGAGAAAAGTGTTGGGTAAAACCCGCAAAGTGCTGGtt                                                                                                                                                                                            >  2:441892/1‑250 (MQ=255)                                                       aaccgCAATACATTCGTGGAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGACAAGACCATCCTGGAAGCTCAGGGTGAAGTTGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGTCGGAACTTACTCGCGCGTGT‑CTGTATATCGAGACGATGGAAAACGTGTTGGGGAACACCCGCAAATGGCTGGTTAAc                                                                                                                                                                                         >  2:88784/5‑251 (MQ=255)                                                         cAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGATCAGGAGCAACGTATCCTCGAAGAGGCGCGTGCGTAGAAGGCCCAGACCCTCCTGGAAGCTCAGGGAGAAGCGGCGCGCTTTGCTAATCTTATGCCGGACAATAAAGCCACGCCGACAATTACTCGCTGGTGT‑CTGGATATCGGAACGATGGAAATAGGGATGGGTaaagccggcaacaagctggttaacgc                                                                                                                                                                                       >  2:133397/1‑227 (MQ=255)                                                                 cATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGGGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTg                                                                                                                                                                               <  1:454786/251‑1 (MQ=255)                                                                    tCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGc                                                                                                                                                                             >  1:10573/1‑250 (MQ=255)                                                                    gcGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGGGCGAACGGTTAGGCGCAACGTATACTCGAAGAGGCGCGTGCGGACAAGGCCAAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTgg                                                                                                                                                                              <  2:343101/248‑1 (MQ=255)                                                                        gAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGTGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACCTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCt                                                                                                                                                                        <  1:153385/251‑1 (MQ=255)                                                                                      taCCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCt                                                                                                                                                                                                <  2:14950/213‑1 (MQ=255)                                                                                      taCCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCt                                                                                                                                                                                                >  1:14950/1‑213 (MQ=255)                                                                                          aaCGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGtt                                                                                                                                                         >  1:144808/1‑248 (MQ=255)                                                                                                       cgcgTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGTGACGATGGAACAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGAGGCAAcc                                                                                                                                                                         >  2:87895/1‑219 (MQ=255)                                                                                                       cgcgTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAAcc                                                                                                                                                                         <  1:87895/219‑1 (MQ=255)                                                                                                          gTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATAGCGAGAAGATGGAAAAAGTGTTGGGTAAAACCCGCAACGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAggg                                                                                                                                      >  2:422562/1‑250 (MQ=255)                                                                                                           ggacaacgttagggCGCAACGTCTCCTCGAAGGGGCGCGTGCGTACAAGGCCCAGTCCATCCTGGAAGCTCAGGGTTAAGTGGCGCGCTTTGCTAAGCTTCTGCCGGAATAGAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CGGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCGGAGGGTTCTGCCGTTAGACCAGATGCTGAAAggtg                                                                                                                                     <  2:218521/239‑1 (MQ=255)                                                                                                                     cAGGCGCAACGTATCCTCGACGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGAAAGCTCAGGGTGAAGAGGAGCGCTTTGCTAAACCTCTGCCGGAGTATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGCTACGATGGAAAAAGTGTTGGGTACCACCCGCAACGTGTTGGTTAACGATAAAGGGGTCACCCTGATGATTCTGCCGATATACCAgctgcttaagggagatagagacac                                                                                                                           >  2:407726/1‑233 (MQ=255)                                                                                                                      aGGCGCAATGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCAGCCGGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCt                                                                                                                          <  2:121455/251‑1 (MQ=255)                                                                                                                      aGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCt                                                                                                                          >  1:492729/1‑251 (MQ=255)                                                                                                                      aGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCt                                                                                                                          >  1:121455/1‑251 (MQ=255)                                                                                                                                   ccTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACGTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCg                                                                                                             <  2:203496/251‑1 (MQ=255)                                                                                                                                    cTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCg                                                                                                             >  1:195464/1‑250 (MQ=255)                                                                                                                                        aaGAGGCGCGTGCGTAAAAGGTCCAGACCATCCGGGAAGCTCAGGGTGAAGGGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCTTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAAc                                                                                                        <  1:337750/251‑1 (MQ=255)                                                                                                                                                       aCAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATct                                                                                          <  1:244701/250‑1 (MQ=255)                                                                                                                                                         aaGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATctgct                                                                                       <  1:262956/251‑1 (MQ=255)                                                                                                                                                          aGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATctgctg                                                                                      <  1:162962/251‑1 (MQ=255)                                                                                                                                                             cccAGTCCAGCCTTGAAGCTCCAGTTGAAGTGGGCCTCTTTGCTAAACTTCTGCCTGAATATAAAGCCGCACTGGAAATTACTCGCGCG‑GTGCTGTATATCGAGACTATGCATAAACCGTTGGGTAACACCCGCAATGTGCTGCTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTTGTAACTTCACTGCGGCGAAGAGCGATATCGGTGCCAGCAAGctgtggcga                                                                                   <  2:139813/251‑7 (MQ=255)                                                                                                                                                                   ccATCCTGTAAGCTCATGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCATCCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACTCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTgccg                                                                             <  1:515755/251‑1 (MQ=255)                                                                                                                                                                     aTCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTgccgcc                                                                           <  1:422562/251‑1 (MQ=255)                                                                                                                                                                          ggAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCct                                                                      <  1:535037/251‑1 (MQ=255)                                                                                                                                                                                tCAGGGAGAAGTGGCGCGCTTTGCTAAACTTCTGACGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACTCGAAAAGTGCTGGTTAACGATAAAGGGGGCAACTTGATGGTTATGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGACGAGCGATAACGGAGCCAGAACTCTGCAGAGACTGCCGCCAGcctcctcct                                                                >  2:387690/1‑251 (MQ=255)                                                                                                                                                                                 cAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATGTAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACCCCCGCAAAGTGCTGGTTAACGATCAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAACGGTGGTAACGCCCCTGAGGAGAAGAGCGATCACGGAGCCAGCAAACTGCTGCGTCTGCCGCCAGCCTCTtcctc                                                               >  2:569601/1‑251 (MQ=255)                                                                                                                                                                                      tGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCcaca                                                           >  1:14030/1‑250 (MQ=255)                                                                                                                                                                                          gtggggagttgtgctaacattcggCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGGTGGGTACCTCCCGCAAAGTGCTGGTTAACGAGAAAGGTGGTAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAg                                                      <  2:128082/228‑1 (MQ=255)                                                                                                                                                                                           tGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAAcccgt                                                     >  1:567112/1‑248 (MQ=255)                                                                                                                                                                                              cgcgcTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGa                                                  <  1:441892/251‑1 (MQ=255)                                                                                                                                                                                                    ttGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaa                                              <  1:437463/249‑1 (MQ=255)                                                                                                                                                                                                     tGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCTCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTtcc                                                                 <  1:367546/229‑1 (MQ=255)                                                                                                                                                                                                     tGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGagcaagc                                            >  1:351294/1‑250 (MQ=255)                                                                                                                                                                                                     tGCTAAACTTCTGCCGGAATATAAAGCAGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGAGGTAACGCCCCTGCGGCGAAGAGAGATAACGGTGCCAGCAATCTGCTGCGTCAGCCGCCATCCTCTtccagataggcagagcgtagtgtgg                                           >  2:367546/1‑229 (MQ=255)                                                                                                                                                                                                      gCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAgcaagcaa                                          <  1:125624/251‑1 (MQ=255)                                                                                                                                                                                                         aaaCTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAacac                                       <  1:183380/251‑1 (MQ=255)                                                                                                                                                                                                          aaCTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACAcg                                      >  1:73412/1‑251 (MQ=255)                                                                                                                                                                                                                 gcgggAATATAAAGCCGCGCCGGAAATTACCCGCGAGTGG‑CTGTATATCTAGACGAGGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTcca                               <  2:261759/248‑1 (MQ=255)                                                                                                                                                                                                                  caggAATATAAAGCCGCGTCGGAAATTACTCGCGTGTGT‑CTGTATATCGAGACGATGGCAAAAGTGTTGGTTAACACCCGCAAAGTGCTGGTTAACGATAAAGGAGGCAACCTGTTGGTTCTGCCGTTTGACCAGATGCTGAAAGGTGGCAACGACCCTGGGGCTATGAGAGAACACGTTGCCAGCAATTTACTGCGTCTGCAGCTAGACTCTGCCTCATCAACCACTGGATCAAGCAACAcgacctact                               >  2:155547/3‑243 (MQ=255)                                                                                                                                                                                                                   cGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTccacc                             >  1:420000/1‑251 (MQ=255)                                                                                                                                                                                                                       atatAAAGCCGCGACGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGAGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGGGGCAACATGATGGTTCTGCAGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGCGCGATAACGGTGCCAGCAATCTGCGGCGTCTGCCGCCAGCCACTTCCTCCACAACCAGTGGAGCACGCAACAAGTCGTCCACCAggc                         >  2:447504/1‑249 (MQ=255)                                                                                                                                                                                                                         ataAAGCCGCGCCGGAAATTACTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACTCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAg                       <  1:387690/251‑1 (MQ=255)                                                                                                                                                                                                                              gCCGCGCCTGAAATTACTCGCGAGTGT‑TTGTAAATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGTTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCGTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCCGTCAGGGCGa                  <  2:526383/251‑1 (MQ=255)                                                                                                                                                                                                                                              cTCGCGAGTGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGAGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGAGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTCCCTACACAACCAGCGGAGCCAGCAACACGTCGTCCACCAGTCAGGGCGATATTCTGGGCCAActc  >  2:281292/1‑249 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    CGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGT‑CTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGC  >  NC_000913/4402765‑4403252

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