Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,402,938	G→A	A301T (GCT→ACT)	hflK →	regulator of FtsH protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,402,938	0	G	A	98.4%	189.3 / ‑7.0	63	A301T (GCT→ACT)	hflK	regulator of FtsH protease
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (30/32);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GGAAGAGACGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAG  >  NC_000913/4402697‑4403181                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     ggAAGAGACGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTg                                                                                                                                                                                                                                             <  1:472136/251‑1 (MQ=255)        ccGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGgc                                                                                                                                                                                                                                      <  1:328795/250‑1 (MQ=255)        aCGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGCTTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGCAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGgc                                                                                                                                                                                                                                      >  2:342651/1‑251 (MQ=255)        aCGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGgc                                                                                                                                                                                                                                      <  1:29890/251‑1 (MQ=255)                    tATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGCAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGCAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGCAACTCAGGGTGAAGTGGCGCGCTTTGCTa                                                                                                                                                                                                                           >  1:460271/1‑250 (MQ=255)                              gTATCACGCAGATGGCAGTCAACTTCCAGGCAGATCGTCCGCCGGAAGAC‑GTAAAAGCGGCGTTTGACGATGCGATTGCAGCGCGTGAAAAAGAACAGCAATACATTCGTGAAGCAGAAGCGGATACCAACGCCGTTCAGCCGCGAGCGAACGGTCAGGCGCAACGTATCCACGAAGAGGAGCGTGCGTACAAGGCCCAGACCCTCCTGGAAACTCAGGGTGAAGCGGCGCGCCTTGCTAAACTTCTGc                                                                                                                                                                                                                  >  2:246244/1‑249 (MQ=255)                               tATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCgg                                                                                                                                                                                                               >  2:330884/1‑251 (MQ=255)                               tATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCgg                                                                                                                                                                                                               >  2:167376/1‑251 (MQ=255)                                     actgctgGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTAAAAGGCCCAGACCATCCTGGAAACTCCGGGAGAAGTGGCGCGCTTTTCTAAACTTCTTCCGGAatat                                                                                                                                                                                                          >  2:159647/2‑250 (MQ=255)                                          tgGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGcc                                                                                                                                                                                                    >  1:241299/1‑251 (MQ=255)                                           ggACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCcg                                                                                                                                                                                                   <  1:244276/251‑1 (MQ=255)                                                       ccAGGCTGCTCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGGCGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCGCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAAtt                                                                                                                                                                                        <  2:267549/250‑1 (MQ=255)                                                         aGGCTGCTAGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAATCAGAAGCGTATCCCACCGACGTTCAGCCGCGTGCGAACGGTAAGGCGCAACGTATCCTCGAAGAGGCGCTGGCGTACAAGGCCCAGACCATACTGGCAACTCAGGGTGAAGAGGAGCGCTTTGCTAAACTTCTGCCGGCatattacgcagggcacggacttccc                                                                                                                                                                                     >  2:353517/1‑230 (MQ=255)                                                         aGGCGGCTCGTCCGCCGGAAGAATGTAAA‑GCGGCGCTTGACGAAGCTATTGTCGCGCGTGAAAACGACCAGCAAACCATTCGTGAAGCAGAAGCGTATACCAACGAAGTGCAGCCGCGTGCGAACGGGCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGGAGGTGCGCGCTTTGCTAAACTTCTGCCGGAATATAATGCCGCGCCGGAAATTACt                                                                                                                                                                                     <  2:436912/251‑1 (MQ=255)                                                            ctgctCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCCGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCAGCGCCGGAAATTACTcgc                                                                                                                                                                                  >  1:26278/1‑251 (MQ=255)                                                                      gccgGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGt                                                                                                                                                                        <  1:454556/251‑1 (MQ=255)                                                                         ggAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGtata                                                                                                                                                                     <  1:330884/251‑1 (MQ=255)                                                                            agaa‑gTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACGCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCg                                                                                                                                                                  <  1:292045/251‑1 (MQ=255)                                                                              aa‑gTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGGGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGCGCGTCTGTATATCg                                                                                                                                                                  >  1:281886/1‑249 (MQ=255)                                                                                        cggcgTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTTGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACTATGGaa                                                                                                                                                       <  2:439754/251‑1 (MQ=255)                                                                                            gTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑g                                                                                                                                                  <  1:296710/251‑1 (MQ=255)                                                                                               tGACGATGCGATAGCCGCGCGTGAAAACGAACAGCAAGACATGCGTGAAGCGGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACGCGCGAGCGTCTGTATATCGAGACGATTGAAAAA‑gtg                                                                                                                                                <  2:237208/250‑1 (MQ=255)                                                                                                   gATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAAAGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCCTCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTAATCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GGGTTggg                                                                                                                                           >  1:255059/1‑251 (MQ=255)                                                                                                      gCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTa                                                                                                                                         <  1:344491/250‑1 (MQ=255)                                                                                                         aTTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAat                                                                                                                                                                                                            >  2:131699/1‑181 (MQ=255)                                                                                                         aTTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAat                                                                                                                                                                                                            <  1:131699/181‑1 (MQ=255)                                                                                                               gcgcgTGAAAACGAACAGCAATACATTCGTGAAGCCGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCCTCCTGGACACTCAGGGGGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAACGCCGCGCCGGAAATTACTCGCGCGCGTCTGTATATCGAGACGATGGAAAAC‑GTGTTGGGTAACCCCCGCa                                                                                                                                >  2:157403/1‑250 (MQ=255)                                                                                                                     gAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAACCTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGAAGATGGCAAAA‑GGGTTGGGTAACACCCGCAAAGTGCt                                                                                                                         >  2:105817/1‑251 (MQ=255)                                                                                                                       aaaCGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTgg                                                                                                                       <  1:313340/251‑1 (MQ=255)                                                                                                                        aaCGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGt                                                                                                                      >  1:469236/1‑251 (MQ=255)                                                                                                                            aaCAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACATATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAAGT‑TTGGGTAACACCCGCAACGTGCGGGTTAAc                                                                                                                  >  2:235998/1‑251 (MQ=255)                                                                                                                                  aatacaggcgggAAGCAGAAGCGTATACCAACGAAGGTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCAAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATaaa                                                                                                            <  2:255059/240‑1 (MQ=255)                                                                                                                                        ttCGAGAAGCAGAAGCGTATACCCACGAAGTTCAGCCGCGTGCGAACGGTAAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGACGCGCCGGAACATACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAACGGGCGGGTTAACGATaagaggggg                                                                                                      >  2:281979/1‑244 (MQ=255)                                                                                                                                          cGTGAAGCAGAAGCGTATACCCACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACg                                                                                                                                                             >  2:66691/1‑195 (MQ=255)                                                                                                                                          cGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAGCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACg                                                                                                                                                             <  1:66691/195‑1 (MQ=255)                                                                                                                                          cGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGGGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCACAGACCATCCTGGATACGCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCTAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGGTAACGATAAAGGTGGCaa                                                                                                    <  2:287984/251‑1 (MQ=255)                                                                                                                                           gTGAAGCAGAAGTGTATACCAACGAAGGTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAAc                                                                                                   <  2:448761/251‑1 (MQ=255)                                                                                                                                            tGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAAtata                                                                                                                                                                                                         <  1:449631/149‑1 (MQ=255)                                                                                                                                            tGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAAtata                                                                                                                                                                                                         >  2:449631/1‑149 (MQ=255)                                                                                                                                            ggAAGCAGAATCGTACTCCAACGAAGTTCAGCCGCGTTCGAACGGTCAGGCGCAACGTATCCGCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCTCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTTTGTATATCGAGACGATGGAAAAA‑GGGTTGGGTATCACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAAcc                                                                                                  <  2:264465/250‑1 (MQ=255)                                                                                                                                             gAAGCAGAAGCGTATTCCAACGAAGTCCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAAcc                                                                                                  <  1:281979/250‑1 (MQ=255)                                                                                                                                             gAAGCAGAAGCGTATACCAAAGACGTTCTGTCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGTCCCAGACCCTCCTGGAAACTCAGGGTGAAGTGTCGCGCTTTGCTAAACTGCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGTGACGATGGAAAAA‑GTGTTGGGTAACACACGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCt                                                                                                 <  2:303907/251‑1 (MQ=255)                                                                                                                                                      gCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCt                                                                                        <  2:241299/251‑1 (MQ=255)                                                                                                                                                               aaCGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGa                                                                               <  2:49767/251‑1 (MQ=255)                                                                                                                                                                        cAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACg                                                                                                                 >  1:130154/1‑208 (MQ=255)                                                                                                                                                                        cAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACGTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACg                                                                                                                 <  2:130154/208‑1 (MQ=255)                                                                                                                                                                           ccGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCCGATGCTGac                                                                   >  2:122444/1‑250 (MQ=255)                                                                                                                                                                                 gccaaCGGTAAGGCGCAACGTATCCCCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACGCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAggtg                                                              <  2:433923/247‑1 (MQ=255)                                                                                                                                                                                     aCGGTCAGGCGCAACGTATCCGCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAAc                                                         <  1:122444/251‑1 (MQ=255)                                                                                                                                                                                     aCGGGCAGGCGCAACGGATCCTCGAAGAGGCGTGTGCGTACAAGGCCCAGCCCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAAc                                                         <  2:150375/251‑1 (MQ=255)                                                                                                                                                                                      cGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACg                                                        >  1:185687/1‑251 (MQ=255)                                                                                                                                                                                          cAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGcccc                                                    >  1:407522/1‑251 (MQ=255)                                                                                                                                                                                                 aaCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGGGGTAACGCCCCTgcggcg                                             >  2:158633/1‑251 (MQ=255)                                                                                                                                                                                                   cGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGa                                            >  1:124293/1‑250 (MQ=255)                                                                                                                                                                                                      aTCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTCAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGACAAA‑GAGTTGGGTAACACCCGCACAGTGCTGGTTAACGATAAAGGGGGCAACCTGCTGGTTCTGCCGTTAGACCCGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAag                                          >  2:372238/1‑249 (MQ=255)                                                                                                                                                                                                      aTCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAag                                          >  1:331177/1‑249 (MQ=255)                                                                                                                                                                                                         cTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCg                                      >  2:6590/1‑250 (MQ=255)                                                                                                                                                                                                                            aCAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATctg                  >  1:444576/1‑251 (MQ=255)                                                                                                                                                                                                                            aCAAGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATctg                  >  1:391159/1‑251 (MQ=255)                                                                                                                                                                                                                              aaGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATctgct                >  1:310974/1‑251 (MQ=255)                                                                                                                                                                                                                              aaGGCCCAGACCATCCTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACCTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATctgct                <  1:480824/251‑1 (MQ=255)                                                                                                                                                                                                                                    cAGACCATCCTGGAAACGCAGGGTGAAGTGGCGCGCTTGGCTAAACTTCGGCCGGAATATAAAGCCGCGCTGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGGTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGAGAACGGTGCCAGCAATCTGCTGCGTCt          <  2:460271/251‑1 (MQ=255)                                                                                                                                                                                                                                            ccTGGAAACTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAg  <  2:185687/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     GGAAGAGACGATTCGTCCGTATGACATGGGTATCACGCTGCTGGACGTCAACTTCCAGGCTGCTCGTCCGCCGGAAGAA‑GTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAA‑GTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAG  >  NC_000913/4402697‑4403181

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: ‑