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

CGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCG  >  NC_000913/941564‑941861                                                                                                                                                      |                                                                                                                                                      cGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTg                                                                                                                                                      <  1:945735/151‑1 (MQ=255)           tccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTa                                                                                                                                              >  1:827694/1‑149 (MQ=255)            ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATc                                                                                                                                            <  2:589441/150‑1 (MQ=255)             cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTaattactt                                                                                                                                          >  2:857657/1‑144 (MQ=255)              gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATCtt                                                                                                                                          >  1:178875/1‑150 (MQ=255)              gcgcAAATTGCGGAAGGGTTGAAGTATACCTACGGCGGCTGGGCAGATGGGAAGCGGTCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAAGGAGTGGCGGTGGGGTGACTTAc                                                                                                                                                <  2:236451/144‑1 (MQ=255)               cgcAAATTGCGGAAGGTTTGAACTATGCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTCCTATCTTGa                                                                                                                                        <  2:1173531/151‑1 (MQ=255)               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGa                                                                                                                                        <  2:23186/151‑1 (MQ=255)               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTATTATCtt                                                                                                                                          >  1:73718/1‑149 (MQ=255)               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTTCTATCttt                                                                                                                                         >  2:604963/1‑149 (MQ=255)                   aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑Ag                                                                                                                                   <  2:362847/151‑1 (MQ=255)                    aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑Agg                                                                                                                                  <  1:907801/151‑1 (MQ=255)                     ttGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCAGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAc                                                                                                                                      <  2:519406/147‑1 (MQ=255)                       gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCa                                                                                                                                <  2:1032596/150‑1 (MQ=255)                             ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGa                                                                                                                         <  1:1248774/151‑1 (MQ=255)                                ttGAACTATAGCTACGGCGGCGGGGCAGATGGCAAGAGCCTGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAGCTTGAAC‑AGGTACGCCAGaaa                                                                                                                       <  2:400879/150‑1 (MQ=255)                                      tataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTATATCTTGAAA‑AGGCACGCCCGAAAACTaa                                                                                                                  >  1:929658/1‑149 (MQ=255)                                        taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGTTTACTATCTTTAAC‑AGGCACACCCGAAAACTAATTcc                                                                                                              >  1:496499/1‑151 (MQ=255)                                         aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGAGGCGGATGGGGGTCTTTTCATCTTTCAC‑CGGGACCCCAGAAATATAATgg                                                                                                               >  1:996326/1‑148 (MQ=255)                                             aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCa                                                                                                          <  2:696857/150‑1 (MQ=255)                                                     tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGACTTACTTTCCTGAAC‑AGGCACGCCCGAAAACTAATGCCCGCCCGatcatc                                                                                                  >  2:1181671/1‑150 (MQ=255)                                                      gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGatcatcat                                                                                                <  2:954742/151‑1 (MQ=255)                                                        gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAC‑AGGCACGCCCGAAAACTAATGCCCGCATGatcatcatc                                                                                               >  1:578863/1‑150 (MQ=255)                                                           gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc                                                                                           <  1:785693/151‑1 (MQ=255)                                                                   cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAA‑AGGGACGCCAGAAATCTAATGCCCCCATGATCATCATCCATCCccgc                                                                                      >  1:963060/1‑148 (MQ=255)                                                                    aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAA‑AAGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCcgc                                                                                        >  1:115853/1‑145 (MQ=255)                                                                    aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTACTATCTTGAAA‑AAGCACGCCCGAAATCTAATGCCCGCACGAACATCCACGGTACCCGCtat                                                                                   >  2:389640/1‑150 (MQ=255)                                                                         gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGac                                                                             <  1:766928/151‑1 (MQ=255)                                                                         gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATTTTTAAA‑AGGCACGCCAGAAAACTAATGCCCCCATGATCATCATCGATCCCCGCCATAcc                                                                                >  1:468334/1‑148 (MQ=255)                                                                          tCGGCTATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCCTGAAC‑CGGGCCCTCCTAAATCTAATTCCCGCACAAAAATCTCCGGTCCGCGCcataa                                                                                 >  2:1249739/1‑145 (MQ=255)                                                                           cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAC‑AGGCACGCCAGAAAACTAAAGCCCGCATGATCATCATCGATCCGCGCTATACCGaca                                                                            >  1:293106/1‑150 (MQ=255)                                                                            ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACAcc                                                                          <  1:118849/151‑1 (MQ=255)                                                                            ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTACTATCTTGAAC‑AAGCACGCCAGAAAACTAATGCCCCCCAGATCATCATCGATCCGCGCTATACCGCCAcc                                                                          >  2:943653/1‑151 (MQ=255)                                                                             gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGTTATGCCGACACCg                                                                         <  2:548413/151‑1 (MQ=255)                                                                              atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATATCTTTAAA‑AAGCACGCCCGAAATCTAAAGACCGCATGATCATCATCGATCCGCGCCATAACGACACCgg                                                                        >  2:1212698/1‑151 (MQ=255)                                                                                atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGTGTGACTTACTATCTTGATC‑AGGCACGCCAGAAAGCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTg                                                                      <  2:17716/151‑1 (MQ=255)                                                                                atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTGAAC‑AGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCgggg                                                                      >  2:189648/1‑149 (MQ=255)                                                                                atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGt                                                                       <  1:604971/150‑1 (MQ=255)                                                                                atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGt                                                                       <  2:914139/150‑1 (MQ=255)                                                                                      aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTTCTATCTTGAAC‑AGGCGCGCCAGAAATCTAATGCCCGCATGGTCATCATCGGTCCGCGCTATACCGACACCGGTGCCGGgc                                                                <  2:844420/151‑1 (MQ=255)                                                                                        cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTACTATCTTTAAC‑AAGCACCCCAGAAAACTAAAGCCCGCATAATAATCCTCGATCCCCGCTATAACCACACCCGTGCCggg                                                                 >  1:1246011/1‑148 (MQ=255)                                                                                         aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGcgcg                                                             <  1:917866/151‑1 (MQ=255)                                                                                           tAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAA‑AGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCCCGCTATACCGACACCGGTGCCGGGggcg                                                             >  1:1083418/1‑149 (MQ=255)                                                                                             aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGa                                                         <  2:1174068/151‑1 (MQ=255)                                                                                                    aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTg                                                   <  1:284518/150‑1 (MQ=255)                                                                                                      tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTTAAC‑AAGCACGCCAGAAATCTAATGGCCGCACGATCATCATCGATCCGCGCTATACCCACAACCGGGGCGGGCGCGGAGATGAggggg                                                 >  2:575000/1‑149 (MQ=255)                                                                                                       gCTGTTTGGTAATAACCCTGGCGAAACGCGAATTAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATc                                               <  2:1058485/151‑1 (MQ=255)                                                                                                         tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATTTTTAAC‑AAGCACGCCCGAAAACTAATGCCCGCCTGATCATCATCGACACGCGCGATAAAAAAACCCGGGCCGGGCGCGGAGAAGAGTgggaccc                                             >  2:812856/1‑146 (MQ=255)                                                                                                               gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAA‑AGGCACGCCAGAAATCTAAATCCCGCCTGATCATCATCGATCCGCGCTATACCCACACCCGGGCCGGGCGCGGAGATGAGTGGATCCCTAtt                                         >  1:9254/1‑149 (MQ=255)                                                                                                                taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTTTATTTTTAACAAGGCA‑GCCAGAAATTTAAATCCCCCGTCATCATCCTAGATCCGCCCCATAACGCCACCCGGGGCGGGGGCGGAGATGGGGGGCTCCCTAtt                                         >  1:704228/1‑148 (MQ=255)                                                                                                                 aataaCCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCTTTTCATCTTGTAG‑AGGCGCGCCAGAAAAATTATTCCCCCCTGATTCACATCGCTCCCCGCCATTACGACCACGGTGGCCGGCGGGGAGAAGAGGGGGGCCCTAtt                                         >  2:141690/1‑147 (MQ=255)                                                                                                                  ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTcc                                    <  2:455794/151‑1 (MQ=255)                                                                                                                    aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCg                                   <  2:649862/150‑1 (MQ=255)                                                                                                                         tGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAC‑AGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGGGCGAAGATGAGGGGGTCCCTATTCGTCCGGGGAc                              >  1:1043602/1‑150 (MQ=255)                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACa                             <  1:1142880/150‑1 (MQ=255)                                                                                                                             gAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAC‑AGGCCCGCCCGAAATCTAATTCCCCCATGATCATCATCGATCCGCGCTATAACGACACCGGGGCCGGGCGCGGAGATGAGTGGATCCCCATTCGGCCGGGaa                               >  1:555643/1‑145 (MQ=255)                                                                                                                                 cgcgAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGTGCCGGGCGCGAAGATGAGGGGATCCCTATTCGTCCGGGGACAGATGcc                       >  1:492243/1‑149 (MQ=255)                                                                                                                                      aTGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCACTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTgg                <  2:1004065/151‑1 (MQ=255)                                                                                                                                         aGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa             <  2:827680/151‑1 (MQ=255)                                                                                                                                              cGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt        >  1:863647/1‑151 (MQ=255)                                                                                                                                                   gggTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc   <  1:118687/151‑1 (MQ=255)                                                                                                                                                   gggTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc   <  2:1050190/151‑1 (MQ=255)                                                                                                                                                    ggTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCg  <  2:426448/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAC‑AGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCG  >  NC_000913/941564‑941861

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: ‑