Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,609,458	T→C	R324R (CGT→CGC)	hyfG →	hydrogenase 4, subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,609,458	0	T	C	98.8%	306.7 / NA	85	R324R (CGT→CGC)	hyfG	hydrogenase 4, subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (49/35);  total (49/36)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACA  >  NC_000913/2609311‑2609604                                                                                                                                                    |                                                                                                                                                   gATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCaa                                                                                                                                                  <  2:491723/150‑1 (MQ=255)  aTTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCaaaa                                                                                                                                                <  2:730461/151‑1 (MQ=255)    tCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAAc                                                                                                                                               <  2:946037/150‑1 (MQ=255)      tgctgGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCt                                                                                                                                             >  1:616691/1‑150 (MQ=255)         tgGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACgg                                                                                                                                         >  1:791077/1‑151 (MQ=255)          ggAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGt                                                                                                                                        <  1:321288/151‑1 (MQ=255)           gAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTc                                                                                                                                       >  2:740039/1‑151 (MQ=255)            aaGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCTGGTCGCGCAAAACCTACGGTCt                                                                                                                                      <  1:604730/151‑1 (MQ=255)             aGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTg                                                                                                                                     <  1:1252218/151‑1 (MQ=255)               tCGAACGGCTACGCAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGa                                                                                                                                    <  2:8464/150‑1 (MQ=255)               tCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGa                                                                                                                                    >  2:608629/1‑150 (MQ=255)               tCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGa                                                                                                                                    >  1:816036/1‑150 (MQ=255)                cGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAAt                                                                                                                                  >  2:877885/1‑151 (MQ=255)                cGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAAt                                                                                                                                  <  1:162222/151‑1 (MQ=255)                  aaCGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCt                                                                                                                                <  1:77038/151‑1 (MQ=255)                  aaCGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCt                                                                                                                                <  1:666948/151‑1 (MQ=255)                    cGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTg                                                                                                                               >  1:756041/1‑150 (MQ=255)                     ggCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGAt                                                                                                                             >  2:288884/1‑151 (MQ=255)                        tACGCAGTCATTTGCTTAACCTGGGCCTCGCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAAGTTTTCCGCGCGCGGGAAAAGTCGATGACGATGGAGGAATTGCTGATCGGGTCGGGTAAAACCTACGGTATGAAGCTGAtt                                                                                                                            <  1:240979/149‑1 (MQ=255)                         acacAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATtg                                                                                                                           <  1:1019031/149‑1 (MQ=255)                          cacaGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATtggtg                                                                                                                        <  2:46686/151‑1 (MQ=255)                           acaGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGCATCTGATtggtgg                                                                                                                       >  1:1202480/1‑151 (MQ=255)                            aaGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATtggtgg                                                                                                                       >  2:1180990/2‑150 (MQ=255)                                cATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTc                                                                                                                  <  2:245374/151‑1 (MQ=255)                                cATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTc                                                                                                                  >  2:12853/1‑151 (MQ=255)                                cATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTc                                                                                                                  <  1:857690/151‑1 (MQ=255)                                 aTTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGTGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcg                                                                                                                 <  2:708880/151‑1 (MQ=255)                                 aTTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcg                                                                                                                 >  1:432965/1‑151 (MQ=255)                                   ttGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcgcc                                                                                                               >  2:190293/1‑151 (MQ=255)                                      cTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTcgccgc                                                                                                             <  1:240770/150‑1 (MQ=255)                                        tAACCTGGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCGGATTGGTGGTGTTCGCcgcg                                                                                                            <  2:685735/149‑1 (MQ=255)                                         aaCCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGat                                                                                                          >  1:795433/1‑150 (MQ=255)                                         aaCCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGat                                                                                                          >  2:651386/1‑150 (MQ=255)                                           ccTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATAtt                                                                                                       >  2:89225/1‑151 (MQ=255)                                               ggCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATtctc                                                                                                    >  1:65453/1‑150 (MQ=255)                                                  ctctCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCaa                                                                                                  >  1:1083129/1‑149 (MQ=255)                                                  atctcCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCaa                                                                                                  >  1:1259116/2‑149 (MQ=255)                                                      ccTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGc                                                                                             >  2:469979/1‑150 (MQ=255)                                                          ccATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCTATATTCTCAAAGAGCAACGt                                                                                        <  1:993109/151‑1 (MQ=255)                                                          ccATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGt                                                                                        <  1:1280184/151‑1 (MQ=255)                                                           cATTTCGATTGCTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCCGGAAAAGACGATGACGATGGCGGAATTGCTCATTGGGTCGCGCAAAACGTACGGTCTGAATCTGATTGGTGGCGTTCGCCGCGATCTTCCCAACGATCacc                                                                                          >  1:106075/1‑146 (MQ=255)                                                             tttCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTg                                                                                     >  1:493514/1‑151 (MQ=255)                                                              ttCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTg                                                                                     <  1:699512/150‑1 (MQ=255)                                                                 gTTGGTTTTGATACCGGCTTTATGCAATTTTTACGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCaa                                                                                  >  1:411738/1‑150 (MQ=255)                                                                      ttttGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTg                                                                            >  1:1195289/1‑151 (MQ=255)                                                                       tttGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTg                                                                            >  2:1253553/1‑150 (MQ=255)                                                                         tGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGaa                                                                          >  2:788113/1‑150 (MQ=255)                                                                                     tATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGcgcg                                                               >  1:1178116/1‑149 (MQ=255)                                                                                     tATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGcgcg                                                               >  2:889662/1‑149 (MQ=255)                                                                                      aTGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGcgcg                                                               >  2:771749/1‑148 (MQ=255)                                                                                       tGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGcgcg                                                               >  2:1267011/1‑147 (MQ=255)                                                                                        gCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATg                                                          >  2:1108900/1‑151 (MQ=255)                                                                                         cAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATgc                                                         >  2:13515/1‑151 (MQ=255)                                                                                          aaTTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATgcg                                                        >  1:222433/1‑151 (MQ=255)                                                                                          aaTTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATgcg                                                        <  1:1196126/151‑1 (MQ=255)                                                                                              tttCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATgcgcgc                                                     <  2:224871/150‑1 (MQ=255)                                                                                              tttCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGcc                                                    >  2:386276/1‑151 (MQ=255)                                                                                              tttCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGcc                                                    >  1:679360/1‑151 (MQ=255)                                                                                               ttCCGCGTGCGGGAAAAGTCGGTGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGcc                                                    <  2:590927/150‑1 (MQ=255)                                                                                                tCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCg                                                   <  2:284208/150‑1 (MQ=255)                                                                                                tCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCg                                                   <  2:872961/150‑1 (MQ=255)                                                                                                tCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCg                                                   >  2:1179396/1‑150 (MQ=255)                                                                                                 ccGCGTGCGGGAGAAGTCGATGACGATGGGGGAATTGCTGATCGGGTCGCGCAAAGCCTACGGTCTGAATCTGATTGGTGGTGTTTGCCGCGATATTCTCAAAGAGCTACGTCTGCAAACGCTGAAACTGGTGCGCGAGATgcgcgcgg                                                   <  1:1224053/149‑3 (MQ=255)                                                                                                   gcgTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACgt                                               <  2:67274/151‑1 (MQ=255)                                                                                                   gcgTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACg                                                >  1:305137/1‑150 (MQ=255)                                                                                                       gCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTc                                            <  2:629844/150‑1 (MQ=255)                                                                                                         gggAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACGGGTGCGCGAGATGCGCGCCGACGTGTCGGa                                         <  1:80366/151‑1 (MQ=255)                                                                                                           gAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGc                                       >  2:587784/1‑151 (MQ=255)                                                                                                            aaaaGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCt                                      >  2:195430/1‑151 (MQ=255)                                                                                                             aaaGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTg                                     >  1:1059234/1‑151 (MQ=255)                                                                                                              aaGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTgg                                    >  2:586358/1‑151 (MQ=255)                                                                                                              aaGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTgg                                    <  1:890778/151‑1 (MQ=255)                                                                                                               aGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTTCGCGAGATGCGCGCCGACGTGTCGGAGCTGGt                                   >  1:1260729/1‑151 (MQ=255)                                                                                                                    aTGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAgaga                              <  2:879728/151‑1 (MQ=255)                                                                                                                     tGACGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTag                                 >  1:1133045/1‑147 (MQ=255)                                                                                                                       aCGATGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGAtgc                           >  2:460550/1‑151 (MQ=255)                                                                                                                          aTGGCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGAtgctgc                        >  2:288232/1‑151 (MQ=255)                                                                                                                            ggCGGAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATgctgct                       <  1:1256567/150‑1 (MQ=255)                                                                                                                                gAATTGCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCtgcttgct                   >  1:1264789/1‑150 (MQ=255)                                                                                                                                     gCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCg             <  2:959252/151‑1 (MQ=255)                                                                                                                                     gCTGATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCg             >  2:729404/1‑151 (MQ=255)                                                                                                                                        gATCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAat          >  1:951804/1‑151 (MQ=255)                                                                                                                                         aTCGGGTCGCGCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAat          >  2:661491/1‑150 (MQ=255)                                                                                                                                                cgcgCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGaaca  <  2:1038991/151‑1 (MQ=255)                                                                                                                                                cgcgCAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGaac   <  1:515721/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GATTCTGCTGGAAGTCGAACGGCTACACAGTCATTTGCTTAACCTTGGCCTCTCCTGCCATTTCGTTGGTTTTGATACCGGCTTTATGCAATTTTTCCGCGTGCGGGAAAAGTCGATGACGATGGCGGAATTGCTGATCGGGTCGCGTAAAACCTACGGTCTGAATCTGATTGGTGGTGTTCGCCGCGATATTCTCAAAGAGCAACGTCTGCAAACGCTGAAACTGGTGCGCGAGATGCGCGCCGACGTGTCGGAGCTGGTAGAGATGCTGCTTGCTACGCCGAATATGGAACA  >  NC_000913/2609311‑2609604

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