Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,931,833	A→G	6.1%	intergenic (+23/+32)	ygdG → / ← fucO	flap endonuclease/L‑1,2‑propanediol oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,931,833	0	A	G	6.1%	195.8 / 2.9	65	intergenic (+23/+32)	ygdG/fucO	flap endonuclease/L‑1,2‑propanediol oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (20/41);  new base G (2/2);  total (22/43)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGT  >  NC_000913/2931685‑2931982                                                                                                                                                     |                                                                                                                                                      gCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAc                                                                                                                                                      >  1:1218697/1‑150 (MQ=255) gCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAc                                                                                                                                                      >  1:1224129/1‑150 (MQ=255)  ccGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAcg                                                                                                                                                     <  2:116278/150‑1 (MQ=255)  ccGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAcg                                                                                                                                                     >  2:1535942/1‑150 (MQ=255)    gAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCCAAcgtc                                                                                                                                                   >  2:1283399/1‑150 (MQ=255)    gAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAcgtc                                                                                                                                                   <  2:1415872/150‑1 (MQ=255)           ggCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGACGTTTGCGGCGTATCCg                                                                                                                                            <  2:990194/150‑1 (MQ=255)           ggCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGACGTTTGCGGCGTATCCg                                                                                                                                            <  2:1327905/150‑1 (MQ=255)            gcgcAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCgg                                                                                                                                           <  1:961057/150‑1 (MQ=255)              gcAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCa                                                                                                                                         >  2:938394/1‑150 (MQ=255)              gcAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCa                                                                                                                                         <  1:829073/150‑1 (MQ=255)                    aaTTAGAAACGCATAAAGAGATGGCGTTTCTGTGGCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGCTGGTACGGTAAGGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATca                                                                                                                                   <  2:1220902/150‑1 (MQ=255)                         gAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCa                                                                                                                              >  2:865337/1‑150 (MQ=255)                                  aaaGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTa                                                                                                                     >  2:884795/1‑150 (MQ=255)                                  aaaGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTa                                                                                                                     <  1:865337/150‑1 (MQ=255)                                      agaTGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAg                                                                                                                 <  1:1399535/150‑1 (MQ=255)                                       gaTGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAcg                                                                                                                                                     <  2:833497/113‑1 (MQ=255)                                       gaTGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAcg                                                                                                                                                     >  1:833497/1‑113 (MQ=255)                                         tGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCg                                                                                                              >  1:106713/1‑150 (MQ=255)                                           gCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGt                                                                                                            >  1:577136/1‑150 (MQ=255)                                              tttCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTAGAGCTATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATg                                                                                                         <  1:1464317/150‑1 (MQ=255)                                                 cTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTa                                                                                                      <  1:542722/150‑1 (MQ=255)                                                   gtgCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTaaa                                                                                                    <  2:954572/150‑1 (MQ=255)                                                    tgCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAg                                                                                                   <  1:498375/150‑1 (MQ=255)                                                     gCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGc                                                                                                  <  2:739322/150‑1 (MQ=255)                                                       cgcgATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGGAAACGTCGTATGCGGCATTATCACATCAGCGCATTTACCAGGCGGGATGGGAAAGctg                                                                                                >  2:1484880/1‑149 (MQ=255)                                                       cgcgATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGctc                                                                                                >  1:875439/1‑150 (MQ=255)                                                        gcgATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGctct                                                                                               <  2:809035/150‑1 (MQ=255)                                                          gATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTAc                                                                                             <  1:1535942/150‑1 (MQ=255)                                                          gATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTAc                                                                                             <  2:384441/150‑1 (MQ=255)                                                          gATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTAc                                                                                             <  1:1190968/150‑1 (MQ=255)                                                          gATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTAc                                                                                             <  1:244018/150‑1 (MQ=255)                                                           atatTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACa                                                                                            <  1:938394/150‑1 (MQ=255)                                                                  ccGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCt                                                                                     <  2:1324339/150‑1 (MQ=255)                                                                    gCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCa                                                                                   <  2:677612/150‑1 (MQ=255)                                                                       tACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGc                                                                                >  2:43099/1‑150 (MQ=255)                                                                       tACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCCGATACGCCCCCAACGTGGTATCCGGCATTATCACATCCGCGCATTTACCAGGCGGGATGGTAAAGCTCTACAAAATCCTCAAGc                                                                                >  2:319692/1‑150 (MQ=255)                                                                               gATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTc                                                                        <  2:1218697/150‑1 (MQ=255)                                                                               gATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTc                                                                        <  2:1224129/150‑1 (MQ=255)                                                                                      atatCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCggg                                                                 <  1:620788/150‑1 (MQ=255)                                                                                         tCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTg                                                              <  2:658419/150‑1 (MQ=255)                                                                                         tCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTg                                                              >  1:128781/1‑150 (MQ=255)                                                                                           gACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGcc                                                            >  2:505126/1‑150 (MQ=255)                                                                                              ggCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCAcc                                                         <  1:721428/150‑1 (MQ=255)                                                                                                cAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCgg                                                       <  1:309610/150‑1 (MQ=255)                                                                                                cAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCgg                                                       <  2:468535/150‑1 (MQ=255)                                                                                                 aaTTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGt                                                      <  1:581996/150‑1 (MQ=255)                                                                                                   tttACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTAc                                                    <  2:380446/150‑1 (MQ=255)                                                                                                             ttGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACatcatc                                          <  1:1131897/150‑1 (MQ=255)                                                                                                              tGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATcc                                         <  1:220782/150‑1 (MQ=255)                                                                                                              tGCGGTTGGTACGGTAACGGCGAGCCCGATACGCCGCCAACGTCGTATCCGGCATTATCACATCAGCGCATTTAACAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACCTGATcc                                         >  1:105284/1‑150 (MQ=255)                                                                                                               gCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCa                                        <  2:1201429/150‑1 (MQ=255)                                                                                                                   ttGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTgc                                    <  1:1355041/150‑1 (MQ=255)                                                                                                                   ttGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTgc                                    <  2:875439/150‑1 (MQ=255)                                                                                                                     ggTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTgccg                                  <  1:692828/150‑1 (MQ=255)                                                                                                                       tACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTgccgcc                                >  2:1082246/1‑150 (MQ=255)                                                                                                                                  aatgCCGGATACGACGTTTGCGGCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGc                     >  1:1201975/4‑150 (MQ=255)                                                                                                                                      ccGGATACGACGTTTGCGGCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGa                 >  1:1344700/1‑150 (MQ=255)                                                                                                                                       cGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGaa                >  1:178479/1‑150 (MQ=255)                                                                                                                                         gATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATg              <  1:884795/150‑1 (MQ=255)                                                                                                                                             cgccgcAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCtt          >  2:1278608/1‑150 (MQ=255)                                                                                                                                             cgccgcAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCtt          <  1:1278608/150‑1 (MQ=255)                                                                                                                                                  cAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTg     <  1:1283399/150‑1 (MQ=255)                                                                                                                                                  cAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTg     <  2:128781/150‑1 (MQ=255)                                                                                                                                                  cAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGAGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTg     <  2:106713/150‑1 (MQ=255)                                                                                                                                                     aCGCCGTACCCGCCATTACCCCATCACCGCCCCTACCAGGCGGCCCGGTAAAGCTCTACAAAATCCTCAAGCCTTGCCTCACGCGGGTTGCCACCTGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGt  <  2:282606/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGT  >  NC_000913/2931685‑2931982

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