Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,296,481	G→A	V25M (GTG→ATG)	dolP →	division and outer membrane stress‑associated lipid‑binding lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,296,481	0	G	A	93.4%	169.5 / 1.4	62	V25M (GTG→ATG)	dolP	division and outer membrane stress‑associated lipid‑binding lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (26/31);  total (29/33)
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.98e-01

GCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAGTGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTCGCGCCAAACAGATTGCTATGGGCGTAGACGGTGCCAACGA  >  NC_000913/3296239‑3296725                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     gCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCGCAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTAcc                                                                                                                                                                                                                                              <  1:30922/251‑1 (MQ=255)    aTGACGGCGGCGAACTTGCAGGTTGGTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCgct                                                                                                                                                                                                                                           <  2:438999/251‑1 (MQ=255)       aCGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCgctgct                                                                                                                                                                                                                                        >  1:477648/1‑251 (MQ=255)                            tGTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACt                                                                                                                                                                                                                   >  1:378879/1‑251 (MQ=255)                             gTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTg                                                                                                                                                                                                                  >  2:231784/1‑251 (MQ=255)                               ttgggCCACAGGATTTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATGGCCTGTGCGATCTGATCGATATCACGCTTTTCCCTTACCAGGATGATTAAGGAGAATACATTAAGGCTCTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTTCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAGCTGa                                                                                                                                                                                                                 <  1:75168/248‑1 (MQ=255)                                aGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGACGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGGAACTGAcc                                                                                                                                                                                                               >  2:471510/1‑251 (MQ=255)                                aGGGCCACAGGATGTGGAGATCCGCATTACTTCGCATCGTAGTGCTCGCCTTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAACCGCCACTGAcc                                                                                                                                                                                                               >  2:200504/1‑251 (MQ=255)                                      acaGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCa                                                                                                                                                                                                         <  1:421064/251‑1 (MQ=255)                                        aGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGAACGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAgt                                                                                                                                                                                                       <  2:240062/251‑1 (MQ=255)                                        aGGATGTAGAGATCAGCATTCATTCGCATCGTAGTGCTAGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGCTAACACGCTTTTACCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTAACGTCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGAGTTGCCGCTGCAGTAGTGGGTACCGCTGCTGTGGGTAACAAAGCCGCAACTGACCCACGCAgg                                                                                                                                                                                                       >  2:463108/1‑250 (MQ=255)                                         ggATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAgtg                                                                                                                                                                                                      <  1:333744/251‑1 (MQ=255)                                           aTGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGCTAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCCGTGTc                                                                                                                                                                                                    >  1:329866/1‑251 (MQ=255)                                           aTGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTc                                                                                                                                                                                                    >  1:464929/1‑251 (MQ=255)                                             gtgGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCgg                                                                                                                                                                                                  >  2:140469/1‑251 (MQ=255)                                             gtgGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCg                                                                                                                                                                                                   <  1:378516/250‑1 (MQ=255)                                               gtagaTCCGCATTCCTTCGCTGCGTAGTGCTCTCATTCAGGAAATGCATCTGCTGTCGGTAAATTGCCTGTGCGATTTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATGCATGAAGGCATTATCGCCAATCGCAGTCCTTATGTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTTCCGTAATGGGTACCGCGGCTGTGGGTACCAAAGCCGCATCTTACCCACGCAGTGTCgcca                                                                                                                                                                                                <  2:178264/249‑4 (MQ=255)                                                   aTCCGCATGCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCa                                                                                                                                                                                            <  2:481768/251‑1 (MQ=255)                                                               tCGCCTCGTAGTGCTCGAATTCAGGAAATGCAGATGCTGACGGTAAATTTCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCAGTTTCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGTACGagg                                                                                                                                                                                 <  2:485498/250‑3 (MQ=255)                                                               tCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATg                                                                                                                                                                                 <  1:138869/250‑1 (MQ=255)                                                                     cGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCt                                                                                                                                                                          <  1:96908/251‑1 (MQ=255)                                                                     cGTAGTGCTCGCATTAAGGCAATGCATATGCCGACGGCAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCGGCAACTGTCCCACGCAGTGTCGGCACCCAGG‑CGACGATGGTACCCTg                                                                                                                                                                         >  2:249093/1‑251 (MQ=255)                                                                      gTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTg                                                                                                                                                                         >  1:78278/1‑251 (MQ=255)                                                                               ccATTCAGGAAATGCATATGCTGACGGTAAATTGCCTTTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAGTGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTgcgc                                                                                                                                                                <  1:463108/250‑1 (MQ=255)                                                                                     aGGAAATGCATATGCTGACGGTAAATTGACTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTACTGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAAc                                                                                                                                                          >  2:151718/1‑251 (MQ=255)                                                                                       gAAATGCATATGCTGACGGTAAAATGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAAGCGCTGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAg                                                                                                                                                        <  2:33089/251‑1 (MQ=255)                                                                                        aaaTGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAgc                                                                                                                                                       >  1:245644/1‑251 (MQ=255)                                                                                          aTGCATATGCTGACGGTAAATTGCCTGTGCGACCTGATCGATAACACGCTTTTCCCTCCCCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAgcgc                                                                                                                                                     <  2:419447/251‑1 (MQ=255)                                                                                                     gACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGaa                                                                                                                                            >  1:96716/1‑249 (MQ=255)                                                                                                     gACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGaa                                                                                                                                            >  2:135350/1‑249 (MQ=255)                                                                                                        ggTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACg                                                                                                                                       >  1:337806/1‑251 (MQ=255)                                                                                                                    gtgCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTaag                                                                                                                            <  1:258341/250‑1 (MQ=255)                                                                                                                             gATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATTCATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTaagaaag                                                                                                                        <  1:402204/245‑1 (MQ=255)                                                                                                                             aaTCGATAACACGCTTTTCCCTCACCAGGAGGATTAAGGAGAAGACATGAAGGCTTTATCGCCAATCGCAGTCCTTATTTCCGCGCTGTTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAgcgc                                                                                                                  <  1:589/250‑1 (MQ=255)                                                                                                                                 gATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCGGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCAt                                                                                                              <  2:477648/251‑1 (MQ=255)                                                                                                                                 gATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTGCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCATGTGGACGATGGTACCCTGGAAGTTCTCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCAt                                                                                                              <  1:466255/251‑1 (MQ=255)                                                                                                                                        cGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTa                                                                                                       <  1:200504/251‑1 (MQ=255)                                                                                                                                            tttCCCTCCCCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAAc                                                                                                     <  2:306832/249‑1 (MQ=255)                                                                                                                                            tttCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACg                                                                                                    <  1:342420/250‑1 (MQ=255)                                                                                                                                              cccctcctcAGGATGAGTAAGGATAATTCATGAAGGTCTTATGGACGATCACATGCTTTATTTCCGCGCTGTTGTTGCAAGGTTGTGTTTCCGCTGCCGTAATGGGTTCCGCTGCTGTGGGTACCAAAGCCGCAGGTGACCCACGCAGTGTCGGCAGCCAGGTGGACGATGGTACCCTGGAAGTGCGCTTGAACAGCGCATTGTCGAAAGACGAACTGATTAAGAATGAAGCGCGCATTAATGTATCGGc                                                                                                  <  2:211131/242‑1 (MQ=255)                                                                                                                                                          atgatTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGCCGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTCTCAGGGCaa                                                                                      >  2:153861/1‑250 (MQ=255)                                                                                                                                                             atTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAAc                                                                                                                                                                                                                    <  2:160878/121‑1 (MQ=255)                                                                                                                                                             atTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAAc                                                                                                                                                                                                                    >  1:160878/1‑121 (MQ=255)                                                                                                                                                                  ggAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTgctgct                                                                             <  1:140469/251‑1 (MQ=255)                                                                                                                                                                          cATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAACGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGCCAAAGTGCTGCTGGTTGGGc                                                                     >  2:9075/1‑251 (MQ=255)                                                                                                                                                                           aTGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTCAGAAAGAAGCGCGCATTAATGTAACGCCCTATCAGGGCAAAGTGCTGCTGGTTGGGc                                                                     >  1:410981/1‑250 (MQ=255)                                                                                                                                                                                        cGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCt                                                       >  1:145270/1‑251 (MQ=255)                                                                                                                                                                                          gcAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGCACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCGGTCACCAAATGCTg                                                      >  2:25465/2‑250 (MQ=255)                                                                                                                                                                                          ccAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGGGCTGCTGGTTGGGCAGTCACCAAATGCTg                                                      >  2:35100/1‑250 (MQ=255)                                                                                                                                                                                           cAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAGTGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCAtt                                                                                                                                                  >  2:139084/1‑157 (MQ=255)                                                                                                                                                                                           cAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAGTGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCAtt                                                                                                                                                  <  1:139084/157‑1 (MQ=255)                                                                                                                                                                                               cGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGAACGATGGTACCCTGGAAGTGCGCGTGCACAGCGCATTGCCGAAAGACGGACAGATTCAGAAAGAAGCGCGCATTAATGTAACGGACTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAAAGCTGACCttt                                                >  2:383066/1‑251 (MQ=255)                                                                                                                                                                                                        tATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTcg                                         <  1:135350/249‑1 (MQ=255)                                                                                                                                                                                                         aTTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGCCGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTCgcgc                                      >  2:48436/1‑251 (MQ=255)                                                                                                                                                                                                         aTTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCa                                                                <  2:289967/225‑1 (MQ=255)                                                                                                                                                                                                         aTTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCa                                                                >  1:289967/1‑225 (MQ=255)                                                                                                                                                                                                           ttCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTCGCGcc                                     <  1:496698/250‑1 (MQ=255)                                                                                                                                                                                                                  cggctgTTGCAAGGTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTCGCGCCAAACAGAt                             <  2:337806/249‑1 (MQ=255)                                                                                                                                                                                                                              ggTTGTGTTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATAGTCGAAAGACGAACAGATTAAGCAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCCAATGCTGCACTTTCGGCTCGCGCCAAACAGATTGCTATGGGCGt                 >  1:64399/1‑251 (MQ=255)                                                                                                                                                                                                                                 tgggtTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTCGCGCCAAACAGATTGCTATGGGCGTAg               <  2:96716/247‑1 (MQ=255)                                                                                                                                                                                                                                    gtTGCCGCTGCCGTAATGGGTACCGCTGCTGTGGGTAACAAAGCCGCTACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGTAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAAAAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTCGCGCCAAACAGATTGCTATGGGCGTAGACgg           <  1:139517/251‑1 (MQ=255)                                                                                                                                                                                                                                               cGTAATGGGTAACGCTGCTGTGGGTACCAAAGCCGCCACTGACCCAAGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGACGTGCGCGTGAACAGCGCATTGCCGAAAGACGAACAGATTAAGAAAGAAGCGCGCCTTAATGTAACGGCCTCTCAGGGCAACGGGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTCGCGCCCAACAGATTGCTATGGACGTAGACAGTGCCACCGa  >  2:184032/1‑249 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     GCTATGACGGCGGCGAACTTGCAGGTTTGTTAGGGCCACAGGATGTGGAGATCCGCATTCCTTCGCATCGTAGTGCTCGCATTCAGGAAATGCATATGCTGACGGTAAATTGCCTGTGCGATCTGATCGATAACACGCTTTTCCCTCACCAGGATGATTAAGGAGAATACATGAAGGCATTATCGCCAATCGCAGTCCTTATTTCCGCGCTGCTGTTGCAAGGTTGTGTTGCCGCTGCCGTAGTGGGTACCGCTGCTGTGGGTACCAAAGCCGCAACTGACCCACGCAGTGTCGGCACCCAGGTGGACGATGGTACCCTGGAAGTGCGCGTGAACAGCGCATTGTCGAAAGACGAACAGATTAAGAAAGAAGCGCGCATTAATGTAACGGCCTATCAGGGCAAAGTGCTGCTGGTTGGGCAGTCACCAAATGCTGAACTTTCGGCTCGCGCCAAACAGATTGCTATGGGCGTAGACGGTGCCAACGA  >  NC_000913/3296239‑3296725

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