Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,522,275	T→G	H607P (CAT→CCT)	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,522,275	0	T	G	98.6%	233.5 / ‑3.3	70	H607P (CAT→CCT)	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (39/29);  total (40/30)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-01

ACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTA  >  NC_000913/4522046‑4522506                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        aCCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGGAAATCCTGAGGGCTgg                                                                                                                                                                                                                                   >  1:260371/1‑236 (MQ=255)    catcatTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCt                                                                                                                                                                                                                                     >  1:326892/1‑231 (MQ=255)    catcatTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTCAGGGCt                                                                                                                                                                                                                                     >  1:658264/1‑231 (MQ=255)         ttATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCa                                                                                                                                                                                                                           <  1:247816/236‑1 (MQ=255)         ttATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCa                                                                                                                                                                                                                           <  1:475765/236‑1 (MQ=255)           aTCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCAtt                                                                                                                                                                                                                         <  1:245606/236‑1 (MQ=255)            tCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTc                                                                                                                                                                                                                                <  2:252418/228‑1 (MQ=255)            tCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTc                                                                                                                                                                                                                                >  1:252418/1‑228 (MQ=255)            tCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCAttt                                                                                                                                                                                                                        >  1:57790/1‑236 (MQ=255)            tCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGATTCATCTGGCAATTAAGGATCGGGAAGTAAATCCTGAGGGTTGGGTCTGGCAttt                                                                                                                                                                                                                        >  1:582331/1‑236 (MQ=255)                aCTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTACGGATCGGGAAGTAAATCCTGAAGGCTGGGTATGGCATTTact                                                                                                                                                                                                                     >  2:28220/1‑233 (MQ=255)                 cTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAatat                                                                                                                                                                                                                   >  2:233643/1‑236 (MQ=255)                                       agagTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGaa                                                                                                                                                                                             <  2:544202/236‑1 (MQ=255)                                         agTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAgg                                                                                                                                                                                           >  2:479379/1‑236 (MQ=255)                                                       ggTTTCTGGATGTTTTGAATTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGTCAATTTTGTTTACATCAAAAAAACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGc                                                                                                                                                                             <  2:582331/236‑1 (MQ=255)                                                         tgtCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATAtt                                                                                                                                                                                                                  >  1:33854/1‑197 (MQ=255)                                                         tgtCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATAtt                                                                                                                                                                                                                  <  2:33854/197‑1 (MQ=255)                                                         tgtCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGTAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTa                                                                                                                                                                           >  2:76253/1‑236 (MQ=255)                                                           tttGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCTGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTATGGCATTTAATATTGCAGTTGCCTCCTCCTGTGTAGGTAATTGTTCATCGCTAca                                                                                                                                                                         >  2:256982/3‑235 (MQ=255)                                                             tGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcc                                                                                                                                                                                                         >  2:724985/1‑202 (MQ=255)                                                             tGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGcc                                                                                                                                                                                                         <  1:724985/202‑1 (MQ=255)                                                             tGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACgg                                                                                                                                                                        >  1:133223/1‑235 (MQ=255)                                                             tGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGGTCATCGCTACgg                                                                                                                                                                        >  2:756281/1‑235 (MQ=255)                                                                  gTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCg                                                                                                                                                                  >  2:24683/1‑236 (MQ=255)                                                                     tttatttaTGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGaa                                                                                                                                                                                             <  2:662526/206‑1 (MQ=255)                                                                     tttatttaTGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGaa                                                                                                                                                                                             >  1:662526/1‑206 (MQ=255)                                                                     tttatttaTGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAgg                                                                                                                                                               >  2:67157/1‑236 (MQ=255)                                                                            atccTTACATTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATTCAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGa                                                                                                                                                        <  2:133223/233‑1 (MQ=255)                                                                                                        cgcAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGAcaatcaat                                                                                                                            >  2:683895/1‑236 (MQ=255)                                                                                                         gcAAAGCCGGGCCAATTTTGTTTACATCAAAAATACCACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGAcaatc                                                                                                                               <  2:260371/232‑1 (MQ=255)                                                                                                              gCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAAt                                                                                                                      <  1:233643/236‑1 (MQ=255)                                                                                                                  cgcCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTaagtgca                                                                                                                  >  2:555123/1‑231 (MQ=255)                                                                                                                   gcCAATTTTGTTTACATCAGAAATGCAACCGGTCCATGTCCCCCCCCTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAAtt                                                                                                                     <  2:57790/232‑1 (MQ=255)                                                                                                                            gTTTACCTCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATccc                                                                                                        <  2:326892/236‑1 (MQ=255)                                                                                                                                aCATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGt                                                                                                    <  1:24683/236‑1 (MQ=255)                                                                                                                                  atcatAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGc                                                                                                  <  1:67157/236‑1 (MQ=255)                                                                                                                                  atcatAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGc                                                                                                  >  2:517773/1‑236 (MQ=255)                                                                                                                                       aaaTACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGa                                                                                                                                                                                              <  1:756140/139‑1 (MQ=255)                                                                                                                                              aCCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGa                                                                                      >  1:408814/1‑236 (MQ=255)                                                                                                                                               ccGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATAtt                                                                                                                                                                                                                  <  1:189633/111‑1 (MQ=255)                                                                                                                                               ccGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATAtt                                                                                                                                                                                                                  >  2:189633/1‑111 (MQ=255)                                                                                                                                               ccGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTGGCCGTGAAGCTCGCGACAATCAAGTTTCATTTCAATTAAATCCCCGGTGCGTAATTTACCGa                                                                                      >  2:528415/1‑235 (MQ=255)                                                                                                                                                  gTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATgt                                                                                   >  1:447433/1‑234 (MQ=255)                                                                                                                                                       tGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCt                                                                              >  1:343704/1‑235 (MQ=255)                                                                                                                                                                 cTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATTGGTCCTCCGGCCaaagc                                                                   >  1:555161/1‑232 (MQ=255)                                                                                                                                                                 cTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAgcgc                                                                   >  1:99951/1‑236 (MQ=255)                                                                                                                                                                   ccgggCCTGAAGCATTTCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCtt                                                                 <  1:28220/233‑1 (MQ=255)                                                                                                                                                                    aCGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCtt                                                                 >  2:609956/1‑235 (MQ=255)                                                                                                                                                                        ccTGAAGCATTGCCCATTGCCGGGTATCAGCTGGCAATTCAGGATCGGGAAGGAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGc                                                            <  2:300180/236‑1 (MQ=255)                                                                                                                                                                                 ttGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATgg                                                   >  1:515178/1‑236 (MQ=255)                                                                                                                                                                                  tGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGc                                                  >  2:31904/1‑236 (MQ=255)                                                                                                                                                                                      cATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGGAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGa                                               >  1:688035/1‑235 (MQ=255)                                                                                                                                                                                         aGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTGATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCCCATGGCCGTTga                                            >  2:291911/1‑234 (MQ=255)                                                                                                                                                                                         aGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCa                                           <  2:408814/236‑1 (MQ=255)                                                                                                                                                                                         aGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCa                                           <  1:479379/236‑1 (MQ=255)                                                                                                                                                                                           ccGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCt                                                                              >  1:638242/1‑199 (MQ=255)                                                                                                                                                                                           ccGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCt                                                                              <  2:638242/199‑1 (MQ=255)                                                                                                                                                                                              ggTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCtt                                                                 >  2:702120/1‑209 (MQ=255)                                                                                                                                                                                              ggTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCtt                                                                 <  1:702120/209‑1 (MQ=255)                                                                                                                                                                                              ggTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACgcgc                                      >  1:56610/1‑236 (MQ=255)                                                                                                                                                                                               gTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCAcg                                         <  1:683895/232‑1 (MQ=255)                                                                                                                                                                                                   catcTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTa                                 >  2:132980/1‑236 (MQ=255)                                                                                                                                                                                                     tcTGGCATTTCAGGATCGGGAAGTAAATCCTGGGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGTTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGa                               <  1:724556/236‑1 (MQ=255)                                                                                                                                                                                                     tcTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGTGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCAcgc                                        >  2:297221/1‑227 (MQ=255)                                                                                                                                                                                                     tcTGGCAATTCAGGATCGGGAAGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCAcgc                                        <  1:297221/227‑1 (MQ=255)                                                                                                                                                                                                                         aaGTAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGGAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCAt           <  2:688035/236‑1 (MQ=255)                                                                                                                                                                                                                            tAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGTGCCAGTAGAAACGCCCGAGAAACGTGCATCgg        >  1:485110/1‑236 (MQ=255)                                                                                                                                                                                                                            tAAATCCTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATc          >  1:269861/1‑234 (MQ=255)                                                                                                                                                                                                                                  cTGAGGTCTGGGTCTGGCATTTAATATTGCAGTTGCATCCTCCGGTGAAGGTAATTGTTCATCGATACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTa  <  2:447433/236‑1 (MQ=255)                                                                                                                                                                                                                                  cTGAGGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTa  <  1:31904/236‑1 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                        ACCCATCATTATCCGACTGAGCAAATAACGAATGATGCAGAGTCGGATAATTCAGGTGTCTGGATGTTTTTATTTATGCTTTCAGTTTTTATTCATAAAATCGCGCAAAGCCGCGCCAATTTTGTTTACATCATAAATACAACCGGTCCATGTCCCACCACTCACGGCCTGAAGCATTGCCCATAGCCGGGTATCATCTGGCAATTCAGGATCGGGAAGTAAATCCTGATGGCTGGGTCTGGCATTTAATATTGCAGTTGCCTCCTCCTGTGAAGGTAATTGTTCATCGCTACGGGTTCCGAGGAAATTGACTTCGCCGTGAAGCTCGCGACAATCAATTTTAATTTCAATTAAATCCCCGGTGCGTAATTTACCGATGGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTA  >  NC_000913/4522046‑4522506

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