Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	941,713	0	T	G	55.6%	49.6 / 96.7	63	V252G (GTG→GGG)	dmsA	dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (3/25);  new base G (35/0);  total (38/25)
Fisher's exact test for biased strand distribution p-value = 1.34e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-13

ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  NC_000913/941569‑941857                                                                                                                                                 |                                                                                                                                                 actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACtt                                                                                                                                             <  1:396335/150‑1 (MQ=255) actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGagt                                                                                                                                              >  1:395977/1‑147 (MQ=255)   tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTa                                                                                                                                            >  2:944204/1‑149 (MQ=255)       ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCt                                                                                                                                      <  1:277073/151‑1 (MQ=255)          cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTATTATAtt                                                                                                                                     >  2:600140/1‑149 (MQ=255)              aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAg                                                                                                                               <  1:302380/151‑1 (MQ=255)                ttGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTaa                                                                                                                                  >  2:471288/1‑146 (MQ=255)                  gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGGACAGGCa                                                                                                                            <  1:97668/150‑1 (MQ=255)                      aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGGAcccc                                                                                                                        >  2:114944/1‑150 (MQ=255)                              aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTa                                                                                                               <  2:641018/151‑1 (MQ=255)                               aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTaa                                                                                                              <  1:450664/151‑1 (MQ=255)                                    aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAGGGACGCCCGAAATCTAATGccc                                                                                                         >  2:840/1‑151 (MQ=255)                                    aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAAGCACGCCCGAAAACTAATGccc                                                                                                         >  2:111026/1‑151 (MQ=255)                                         cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTTTCTTGAAAAGGCACGCCCGAAAACTAATAccc                                                                                                         >  2:595331/1‑146 (MQ=255)                                           gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTAATATTTTGAACAGGCACGCCCGAAAAATAAAGccc                                                                                                         >  2:36352/1‑144 (MQ=255)                                               cTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCAAGCCAGAAAACTAATGCCCGCATGatca                                                                                                >  1:91686/1‑149 (MQ=255)                                                 gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                             <  1:668070/150‑1 (MQ=255)                                                 gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTAATATCTTTAACAAGCACGCCAGGAATCTAATTCCCGCATGa                                                                                                   >  1:497357/1‑144 (MQ=255)                                                 gggCAGATGGCAACAGCCCGTCGGATATCG‑AAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGTGGGGGGCGTCCTTTCTTGAAAACGCACGACAGAAAGCTaaagggcgcctg                                                                                                    >  1:703984/1‑132 (MQ=255)                                                   gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                          <  1:531694/151‑1 (MQ=255)                                                    cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAGTATCTTTAAAAGGGACGCCAGAGATCTAATAGCCGCATGctcaacttc                                                                                           >  1:332674/1‑144 (MQ=255)                                                    cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCTTACTATCTTTAACAAGCACCCCAGAAATCTAATGCCCGCCTGatcatcctc                                                                                           >  1:290847/1‑149 (MQ=255)                                                     aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGAt                                                                                        <  2:224530/151‑1 (MQ=255)                                                           cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCgcgc                                                                                  <  1:59788/151‑1 (MQ=255)                                                            aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCt                                                                                 <  1:329550/151‑1 (MQ=255)                                                             aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCta                                                                                <  1:600887/151‑1 (MQ=255)                                                             aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCta                                                                                <  1:224742/151‑1 (MQ=255)                                                             aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCgcgc                                                                                  >  2:873851/1‑149 (MQ=255)                                                               aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTATTATTTTTAACAGGCCCGCCAGAAAACTAAAGTCCGCATGATCATTATCGAACcgcg                                                                                   >  1:12257/1‑146 (MQ=255)                                                                    gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGac                                                                         <  1:486144/151‑1 (MQ=255)                                                                    gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATATCTTGAAAAGGCACCCCCGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATAccca                                                                          >  2:476938/1‑148 (MQ=255)                                                                     tCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCCCCTGATCCTCATCGATCCCCGCTATACCGac                                                                         >  2:193855/1‑150 (MQ=255)                                                                      cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTACTATCTTTAACAGGGACGCCAGAAATCTAATGCCCCCCTGATTATCATCGATCCGCGCTATACCGa                                                                          >  1:743985/1‑148 (MQ=255)                                                                        gagatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGATTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCCTCGATCCGCGCTATACCGACACCg                                                                     <  1:518681/148‑1 (MQ=255)                                                                        gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACAcc                                                                      >  2:618985/1‑150 (MQ=255)                                                                             cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAACGACACCgggggc                                                                >  1:782081/1‑147 (MQ=255)                                                                             cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATAACTTTGACAGGCACGCCAGAAATCTAATGCCCCCATGATTATCATCGAACCGCGCTATACCGACACCGGTggc                                                                >  1:934816/1‑149 (MQ=255)                                                                              gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTCTCTTGAAAAGGCACGCCAGAAAACTAATGCCCCCCAGATCCTCCTCGATCCGCGCTATACCCACACCCGGGcc                                                                >  1:186027/1‑150 (MQ=255)                                                                              gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAAACTAATGGCCGCCTGATCATCATCGATCCGCGCTATACCggcccccgggcc                                                                >  1:409673/1‑139 (MQ=255)                                                                               aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCgg                                                              <  2:938288/151‑1 (MQ=255)                                                                                aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTAATATCTTGAACAGGCACGCCAGAAATATAATGCCCGCATGATCATTATCGATCCGCGCTATACCCACCCCGGGGGCggg                                                             >  1:670085/1‑151 (MQ=255)                                                                                  aCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcg                                                           <  2:734739/151‑1 (MQ=255)                                                                                        aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGAATTACTATTTTGAAAAGGGACGCCAGAAAACTAAAGTCCGCATCATCATCATCAATCCGCGCTATAACCACAACGGTGCCGGGcgcg                                                         >  1:327152/1‑147 (MQ=255)                                                                                          cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGTGCGCGAAGAt                                                    <  1:917937/150‑1 (MQ=255)                                                                                             gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGGGCGCAGATGAg                                                 >  2:70631/1‑150 (MQ=255)                                                                                             gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGACCTATCTTTTTGAAAAGGGACCCCTGAAAATTAATTTCCGCGTGATCACCATCTCTCCGCGGTATACCGACACCGGTGCCCGGCGCGCAg                                                      >  2:633283/1‑145 (MQ=255)                                                                                                  gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATc                                           >  2:244377/1‑151 (MQ=255)                                                                                                     gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATcc                                          <  2:854588/149‑1 (MQ=255)                                                                                                     gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCTTGAAAAAGCACGCCCGAAATCTAATGCCCGCCTGATTCTCATTGATCCGCGCTATACCCAAACCCGGGGCGGGGGCGGAGATGAGGGGATgcc                                         >  1:477343/1‑150 (MQ=255)                                                                                                       ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTa                                       <  2:388778/150‑1 (MQ=255)                                                                                                            aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGt                                  <  2:820893/150‑1 (MQ=255)                                                                                                            aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                                 <  1:423875/151‑1 (MQ=255)                                                                                                               aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCgg                              <  1:985658/151‑1 (MQ=255)                                                                                                                  ccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTGAAAAGGCACGCCCGAAAACCAAAGCCCGCATGATCATCCTCGATCCCCCCTATACCGCCACCAGGGCCGGGGGGGGAGATGGTGGGGTCCCCATTTGTCCcgg                             >  1:525196/1‑149 (MQ=255)                                                                                                                  ccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCCTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCCCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGGGGATCCCTATTCGTCCggg                             >  2:128935/1‑149 (MQ=255)                                                                                                                  ccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGCTACTATATTGAACACGCACCCCCGAAATCTTATGCCCGCATGATCATCATCGAGCCCCGCTATACCCACACCCGTGCCGGGCGCGAAGAAGAGGGGAGCCCTAtt                                     >  2:718821/1‑141 (MQ=255)                                                                                                                    tGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTACTATCTTGAAAAGGCACGCCAGAAATCTAATGCCCGCATGAACATCATCGATCCGCGCTATAACGACACCGGTGCCGGGGGCGGAGATGAGGGGATCCCTATTCGTCCgggga                           >  2:785322/1‑147 (MQ=255)                                                                                                                      gCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATATATTGAAAACGGAAGGCAGAAAAATAATGGCCGGATGATCATTATCGATCCGCGCTAAAAAGACACCGGGGGCGGGCGGGAAGATGAGTGGAGTCCTTTTCTTCCGGGAACa                         >  1:295199/1‑149 (MQ=255)                                                                                                                       cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                      <  2:605707/151‑1 (MQ=255)                                                                                                                             gcgAATGAGTGGCGGTGGGGGGACTTTCTTTCTTGAAACGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGc                 >  1:422180/1‑150 (MQ=255)                                                                                                                                       ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc       <  2:300192/150‑1 (MQ=255)                                                                                                                                           gTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt  >  2:206738/1‑151 (MQ=255)                                                                                                                                           gTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATACCGCACTGGTTAACGGTCt  >  2:764619/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                 ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  NC_000913/941569‑941857

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