Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,023,189	C→T	R232Q (CGG→CAG)	yccS ←	putative transporter YccS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,023,189	0	C	T	97.4%	126.7 / NA	38	R232Q (CGG→CAG)	yccS	putative transporter YccS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/17);  minor base G (0/1);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCGGCTGCTGAT  >  NC_000913/1022920‑1023474                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                              cGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCTTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTCGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTca                                                                                                                                                                                                                                                                 <  2:224588/300‑1 (MQ=255)  gTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATNATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCNGGTcag                                                                                                                                                                                                                                                                >  2:580929/1‑300 (MQ=255)                           gcaTTAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAAt                                                                                                                                                                                                                                       >  2:126170/1‑300 (MQ=255)                           gcaTCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAAt                                                                                                                                                                                                                                       >  1:148826/1‑300 (MQ=255)                            caTCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACANCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATg                                                                                                                                                                                                                                      <  2:634277/300‑1 (MQ=255)                                   atGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCAt                                                                                                                                                                                                                               >  2:247175/1‑300 (MQ=255)                                     gCGTAAAAACGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGAAGCGTTCCCTGTTGACCACAATCGCCACGTCAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATcc                                                                                                                                                                                                                             >  2:344074/1‑299 (MQ=255)                                            aGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGAcc                                                                                                                                                                                                                      <  1:71150/300‑1 (MQ=255)                                                     aaaaTGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTTTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCgag                                                                                                                                                                                                             <  1:353429/300‑1 (MQ=255)                                                     aaaaTGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAANGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCgag                                                                                                                                                                                                             <  2:168756/300‑1 (MQ=255)                                                                 aTGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCAAGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCAATCAGCAGACCGTTGGCGAGAGCCAAATCGTa                                                                                                                                                                                                 >  2:1544/1‑300 (MQ=255)                                                                  tGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTAc                                                                                                                                                                                                <  2:148204/300‑1 (MQ=255)                                                                       aTAAGGCTGAAGCAACAAAATGCAGCTTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACTTCGCTGTGGCGAAAATGTTCACGCANTGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAATTAANAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATTGTACAGCgg                                                                                                                                                                                           <  2:571662/300‑1 (MQ=255)                                                                          aGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGc                                                                                                                                                                                        >  1:326867/1‑300 (MQ=255)                                                                            gCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCt                                                                                                                                                                                      <  1:129796/300‑1 (MQ=255)                                                                                 cGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCtt                                                                                                                                                                                 <  1:254252/300‑1 (MQ=255)                                                                                     acaaAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTAAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTCGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGGGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTCGGCGAGAGCCAAATCGCACAGCGGTGCCTGGCTTTGAt                                                                                                                                                                             >  2:128016/1‑300 (MQ=255)                                                                                                                         tcTGCATCGACATCATCCGCTGGAAACGGAACAGCACGTCACTTTGGCGAGAACGTTCACGCAATTTGGGGGACTGAATATGATGAA‑AGCGTGCACGCTCGTGAATAGCCGGGGCGATAAATTAATAATGCAGCGTGCTACGCGTCCCCTGTCGACCACGATCGCCACGGAAGCGGTTCAGAAGCGATAGTTTCGTCTGAGTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAGTATCAGGATCAAACaggctctatttgagcc                                                                                                                                         <  2:232152/299‑16 (MQ=255)                                                                                                                             cATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAg                                                                                                                                     <  1:327708/300‑1 (MQ=255)                                                                                                                                                    ggAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGc                                                                                                              <  1:379748/300‑1 (MQ=255)                                                                                                                                                     gAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCa                                                                                                             >  1:305527/1‑300 (MQ=255)                                                                                                                                                        cAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAAcg                                                                                                          >  1:54071/1‑300 (MQ=255)                                                                                                                                                            ctgcggggtgtgtataagagacagcGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGcc                                                                                                      >  1:15881/25‑300 (MQ=255)                                                                                                                                                                cGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATTA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGt                                                                                                  <  1:609906/300‑1 (MQ=255)                                                                                                                                                                cGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGt                                                                                                  <  1:547450/300‑1 (MQ=255)                                                                                                                                                                                        aTGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGa                                                                          >  1:419443/1‑300 (MQ=255)                                                                                                                                                                                          gTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAAc                                                                        >  1:234311/1‑300 (MQ=255)                                                                                                                                                                                                                gagCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTc                                                   >  1:420599/1‑300 (MQ=255)                                                                                                                                                                                                                     ggCACGCTCGTGAATATCCTGTGCGATAAAGTAATAATGCAGCGTGCGACGCGTTCCCGGTTGACCACGATCGCCACGGAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGACAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGAc                                              <  2:348718/300‑1 (MQ=255)                                                                                                                                                                                                                           cTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGAACAATAAGTGTAAGGACGTTGTa                                        >  2:253458/1‑300 (MQ=255)                                                                                                                                                                                                                           cTCGTGAATATCATGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGTGAGGGTTTCGTCTGATCCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATNTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTCGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTATCACGACGtcata                                        >  2:171239/1‑296 (MQ=255)                                                                                                                                                                                                                               tGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAg                                    >  1:327189/1‑300 (MQ=255)                                                                                                                                                                                                                               tGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGACGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAg                                    >  2:607384/1‑300 (MQ=255)                                                                                                                                                                                                                               tGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGACATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCCACGCGCCAGGTTGTCCTGNAGCGGGCGGACCGGGAACAGCAGATGAACCATAAGTGTCAGGACGTTGTACCAg                                    >  2:525327/1‑300 (MQ=255)                                                                                                                                                                                                                                 aaTATCCTGTTCGACAAAGTAATAATACAGCGTGCGACGCTTTCCCTGTTGATCACAATCGCCACGTAAGCGGGGCAGTAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGGTGGCGAGAGCCAAATCATACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTGGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGAc                                  <  2:154375/300‑1 (MQ=255)                                                                                                                                                                                                                                    atCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCTGTTGTCCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCANACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCANACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGc                               <  2:571896/300‑1 (MQ=255)                                                                                                                                                                                                                                                   gTAATAATGCAGCGTGCGACGCGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATa                >  1:313130/1‑300 (MQ=255)                                                                                                                                                                                                                                                                 tGCGACGAGTTCCCTGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATTAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCATATACATCGGCTGCTGAt  <  2:53961/300‑1 (MQ=255)                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                              CGTTATCGCGCATCCGCTCCAGCGCAGCATCAATATGCGTAAAAGCGCGCTCAAAATGCGGATCATGTTGATAAGGCTGACGCAACAAAATACAGCGTGACAGTTGCTGGCACGCCTGGCCCTGCATCGACATCAGCCGCTGAAAACGGAACAGCACGTCGCTGTGGCGAAAATGTTCACGCAATGTTTGATACTGAATATGA‑GAAGAGCTGGCACGCTCGTGAATATCCTGTGCGACAAAGTAATAATGCAGCGTGCGACGCGTTCCCCGTTGACCACGATCGCCACGTAAGCGGGTCAGCAGCGAGAGTTTCGTCTGATTCAATGTCGCCATCAGCAGACCGTTGGCGAGAGCCAAATCGTACAGCGGTGCCTGGCTTTGATCTTCAATATCAGGATCAAACATGCGCGACTTGAGCTCAAGATAACGCGCCAGTTGTTCATAGCAACGCGCCAGGTTGTCCTGCAGCGGGCGGACCGGGAACAGCAGATGACCAATAAGTGTCAGGACGTTGTACCAGACGGCACCGGCCAGCAGATACATCGGCTGCTGAT  >  NC_000913/1022920‑1023474

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