Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	941,713	0	T	G	45.2%	89.2 / 26.4	62	V252G (GTG→GGG)	dmsA	dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (7/27);  new base G (28/0);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 2.59e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGG  >  NC_000913/941569‑941859                                                                                                                                                 |                                                                                                                                                     actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                                 >  1:775283/1‑150 (MQ=255)    cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTa                                                                                                                                             <  1:174709/151‑1 (MQ=255)     ctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                                 >  2:318679/1‑146 (MQ=255)       ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAt                                                                                                                                            <  1:965323/149‑1 (MQ=255)       ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑t                                                                                                                                        <  1:738019/151‑1 (MQ=255)         gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTg                                                                                                                                      <  1:391100/151‑1 (MQ=255)           gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTAt‑t‑ttt                                                                                                                                      >  2:605168/1‑148 (MQ=255)            cAAATTGCGGAAGGTTTGAACTATACCTGCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTAT‑T‑TTGaa                                                                                                                                    >  1:5257/1‑150 (MQ=255)             aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTAT‑T‑TTGaa                                                                                                                                    >  2:339855/1‑149 (MQ=255)             aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTAT‑C‑TTTaa                                                                                                                                    >  1:32699/1‑149 (MQ=255)                ttGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTAT‑C‑TTGAACAgg                                                                                                                                >  1:786153/1‑150 (MQ=255)                   cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACg                                                                                                                            <  2:98851/151‑1 (MQ=255)                      aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTTTAT‑C‑TTTaaaaaggaaccc                                                                                                                          >  1:551993/1‑140 (MQ=255)                        ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGGTTAATAT‑C‑TTTAAAAAGGACGcccg                                                                                                                        >  2:218829/1‑148 (MQ=255)                             gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATc                                                                                                                   <  1:465874/150‑1 (MQ=255)                              aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTAT‑C‑TTGAAAAAGCACCCCCGaaaac                                                                                                                   >  2:280263/1‑147 (MQ=255)                                 tataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTAT‑C‑TTTAACCGGGACGCCCGAAATCTaa                                                                                                                >  1:945961/1‑149 (MQ=255)                                  ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAAt                                                                                                               <  1:458585/149‑1 (MQ=255)                                     ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTAT‑C‑TTGTAAAAGCACGCCCGAAAACTaaaagc                                                                                                            >  2:647874/1‑145 (MQ=255)                                          ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTAT‑C‑TTGAAAAAGCACGCCCGAAAACTAATGCCCGCCTg                                                                                                      >  2:227099/1‑150 (MQ=255)                                           gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTAT‑C‑TTTAACAAGGACGCCCGAAAACTAATGccc                                                                                                           >  1:568649/1‑144 (MQ=255)                                             ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTTT‑C‑TTTAACAAGCACGCCCGAAAACTAATGCCCGCATGat                                                                                                    >  2:918727/1‑149 (MQ=255)                                               cTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcat                                                                                                 <  2:525913/150‑1 (MQ=255)                                                tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTAT‑C‑TTGAACAGGCACGCCAGAAAACTAATGCCCGCATGatcatc                                                                                                >  1:400904/1‑150 (MQ=255)                                                   gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                            <  1:349913/151‑1 (MQ=255)                                                    cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTTT‑C‑TTTAAAAAGGACCCCCGAAATCTTAAGCCCGCCTGat                                                                                                    >  1:345905/1‑142 (MQ=255)                                                    cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTAT‑C‑TTGAACAGGCCCGCCAGAAAACTAATGCCCGCATGatc                                                                                                   >  1:769884/1‑143 (MQ=255)                                                       aTGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATcc                                                                                        <  1:393410/151‑1 (MQ=255)                                                           cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTAT‑C‑TTTAAAAAGCACGCCCGAAAACTAAAGCCCGCCTGATTATCATCGATTCgcgc                                                                                    >  2:591285/1‑151 (MQ=255)                                                             aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCta                                                                                  <  2:257609/151‑1 (MQ=255)                                                               aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtata                                                                                <  2:13608/151‑1 (MQ=255)                                                               aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTAT‑C‑TTTAACAAGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCtat                                                                                 >  1:66422/1‑150 (MQ=255)                                                                  ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCg                                                                             <  2:20108/151‑1 (MQ=255)                                                                         atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTAT‑C‑TTTAAAAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCCCGCTATACCGacac                                                                         >  1:619647/1‑148 (MQ=255)                                                                                aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATATCCCTGGCGAATCGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGATCAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCgg                                                                <  1:504018/150‑1 (MQ=255)                                                                                 aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgc                                                              <  2:66091/151‑1 (MQ=255)                                                                                    aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCGTTTCAT‑C‑CTGAAAAGGCACGCCAGAAATTTAAATCCCCCCTGATTATCCACCATCCCCGGCATACCCACACCCGTGCCggggg                                                             >  2:537103/1‑147 (MQ=255)                                                                                        aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGa                                                       <  2:780777/151‑1 (MQ=255)                                                                                          cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGATTTATTAT‑C‑TTTTACAGGGACCCCAGAAAAATAAATCCCGCAAGATCATCATCAAACCCACCCGTACCGACCACGGCGGCGGGCGCGAAGAt                                                      >  2:176855/1‑150 (MQ=255)                                                                                          cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTAT‑C‑TTTAAAAAGCACGCCCGAAATCTAATGCCCGCCAGATCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGGGCGgaga                                                       >  2:130366/1‑149 (MQ=255)                                                                                           tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTTCTAT‑C‑TTTTACAGGCACGACCGAAATCTAATGCCCGCATGAACATCATCTATCCTCGCTATACCGACACCGGGGCCcgg                                                               >  2:624782/1‑140 (MQ=255)                                                                                           tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTAT‑C‑TTGAAAAAGCACGCCAGAAATCTAATGCCCGCATGATCATCCTCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGAt                                                      >  2:666017/1‑149 (MQ=255)                                                                                                gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTAT‑C‑TTTAAAAGGCACGCCAGAAATCTAATGCCCGCCAGATCATCATCGATCCGCGCCATACCCACACCGGTGCCGGGGGCGGAGAAGAGTgg                                                >  1:979564/1‑150 (MQ=255)                                                                                                gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTAT‑C‑TTTAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGGAGATGAGTgg                                                >  1:188942/1‑150 (MQ=255)                                                                                                gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCTTCCTTT‑CTTTGAAACGGAAC‑CCAGAAAACTGATGCCCCCATGCACTTACTCCGTCCGCCCCTAAACAACAACGGTGGCGGGGGCGAAGATGAgggg                                                >  1:536054/1‑150 (MQ=255)                                                                                                       ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTAT‑C‑TTGAAAAAGCCCGCCCGAAATCTAAAGCCCGCATGAACATCATCGATCCCCGCTATACCGACACCCGGGCCGGGGGCGGAGATGAGTGGATccc                                           >  1:473505/1‑148 (MQ=255)                                                                                                          gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCg                                     <  1:255267/151‑1 (MQ=255)                                                                                                          gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCCATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCg                                     <  2:395848/151‑1 (MQ=255)                                                                                                            aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                                   <  1:512909/151‑1 (MQ=255)                                                                                                               aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCgg                                <  2:81709/151‑1 (MQ=255)                                                                                                               aaCCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGGCTTTCTTTCC‑TTAAAAAGGACCCCCGGAAAT‑TAATGCCCGCCAGATCACCAACGATCCGCGCTATACCCCCACCCGGGGCGGGGGGGGAAAAGAGTGGGTcaccttttgtcc                                  >  2:276990/1‑138 (MQ=255)                                                                                                                     ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                            <  1:812703/149‑1 (MQ=255)                                                                                                                     ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                            <  1:75427/149‑1 (MQ=255)                                                                                                                      gCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGa                         <  1:318680/151‑1 (MQ=255)                                                                                                                                   gAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGtt            <  1:394026/151‑1 (MQ=255)                                                                                                                                   gAGTGGCGGTGGGGGGACTTACTAT‑C‑TTTAACAGGCACGCCAGAAATCTAATGCCCGCAAGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTTGTCCGGGGACAGATGCCGCACTgggt            >  1:177020/1‑149 (MQ=255)                                                                                                                                      tGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc         >  2:333077/1‑151 (MQ=255)                                                                                                                                       ggCGGTGGGGTGACTTGCTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc         <  1:34748/150‑1 (MQ=255)                                                                                                                                          ggTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc     <  2:786148/151‑1 (MQ=255)                                                                                                                                            tGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg   >  1:691437/1‑151 (MQ=255)                                                                                                                                             ggggTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg  <  1:227100/151‑1 (MQ=255)                                                                                                                                             ggggTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg  >  1:261786/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                     ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAT‑C‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGG  >  NC_000913/941569‑941859

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