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

CGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941578‑941862                                                                                                                                        |                                                                                                                                                      cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATATCtt                                                                                                                                          >  2:684531/1‑149 (MQ=255)         gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCa                                                                                                                                 <  1:216178/150‑1 (MQ=255)          cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTACTATCTTTAACAgg                                                                                                                                   >  2:972180/1‑147 (MQ=255)               ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGa                                                                                                                          <  2:143144/151‑1 (MQ=255)               ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTATCTTTaaaaaggaagcc                                                                                                                             >  1:653029/1‑138 (MQ=255)                         ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAAt                                                                                                                  <  2:1122454/149‑1 (MQ=255)                          taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGGAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCCTAATATTTTTAAAAGGGAccccccaaaactcaaggc                                                                                                               >  1:271866/1‑137 (MQ=255)                            ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCg                                                                                                             <  1:386360/151‑1 (MQ=255)                            ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTTATATTTTTAAAAGGGCCGCCCGAAAAATAAAGccc                                                                                                              >  1:499824/1‑150 (MQ=255)                              tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTATATCTTTAACACGGCCGCCCGAAAACTaaaagccgc                                                                                                            >  1:268277/1‑143 (MQ=255)                               aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGc                                                                                                            <  2:1103144/149‑1 (MQ=255)                                 ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGa                                                                                                        <  2:1079895/151‑1 (MQ=255)                                 ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAAGCACGCCCGAAATCTAATGccccc                                                                                                            >  2:294383/1‑145 (MQ=255)                                  gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGat                                                                                                       <  2:299941/151‑1 (MQ=255)                                   cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc                                                                                                      <  1:990813/151‑1 (MQ=255)                                        gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                                  <  2:326686/150‑1 (MQ=255)                                        gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatc                                                                                                   <  2:572203/149‑1 (MQ=255)                                        gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatc                                                                                                   <  2:1354/149‑1 (MQ=255)                                         ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                                <  1:116230/151‑1 (MQ=255)                                          gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                               <  2:21624/151‑1 (MQ=255)                                          gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTACTATCTTGAACAGGGACGCCAAAAAACTTAAACCCGCATGGACATTATCg                                                                                               >  2:768589/1‑151 (MQ=255)                                           cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGa                                                                                              <  2:1048961/151‑1 (MQ=255)                                            aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTTAACAGGCACGCCCGAAAACTAATGCCCGCCTGATCATCATCGAt                                                                                             >  2:492880/1‑151 (MQ=255)                                            aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGACTTACTCTTTTTAACAGGCACGCCAGAAAACTTATGCCCGCCAGatc                                                                                                      >  1:1100627/1‑142 (MQ=255)                                             gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc                                                                                            <  1:686184/151‑1 (MQ=255)                                                  cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCgcgc                                                                                       <  1:837852/151‑1 (MQ=255)                                                     cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTATTATCTTTAAAACGGCCGCCAGAAAACTAATGACCCCATGAAAATCATCCAATCCCGCtat                                                                                    >  2:1011446/1‑151 (MQ=255)                                                     cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGGCAAAAAACTAATTCCCGCATGATCATCATCGAACCccgc                                                                                       >  2:691123/1‑148 (MQ=255)                                                      aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAAACTAATTCCCGCATGATCATCATCGATCCGCGCcat                                                                                    >  2:399899/1‑150 (MQ=255)                                                          cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTTCTATCTTTAACAGGCACGCCCGAAATCTAAAGCCCCCAAGATCATCATCGAACCGCGCCATAcc                                                                                 >  1:512957/1‑149 (MQ=255)                                                           gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTTTCTTGAACAGGCACGCCAGAAATCTAATTCCCGCCTGATTATCATCGAACCGCGCTATACCGa                                                                               >  2:20824/1‑150 (MQ=255)                                                           gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAATCTAATGGCCGCCTCATCATCATTGATCCGCGCTATACCGac                                                                              >  2:817323/1‑151 (MQ=255)                                                              ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACAcc                                                                           <  1:397315/151‑1 (MQ=255)                                                               gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTATACCGACAcc                                                                           >  2:825701/1‑150 (MQ=255)                                                                  atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGt                                                                        <  2:106279/150‑1 (MQ=255)                                                                  atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTg                                                                       <  2:80871/151‑1 (MQ=255)                                                                      aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCCGAAAACTAAAGCCCCCATGATCATCATCGATCCGCGCTATACCGACACCCGGGCCgg                                                                   >  2:1171829/1‑151 (MQ=255)                                                                        aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgc                                                                 <  2:196322/151‑1 (MQ=255)                                                                          cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcgc                                                               <  1:359777/151‑1 (MQ=255)                                                                                 cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                         >  2:35064/1‑150 (MQ=255)                                                                                  tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCCCGCTATACCGACACCGGGGCCGGGGGCGAAGAt                                                         >  2:557770/1‑149 (MQ=255)                                                                                   ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAg                                                      <  1:319419/151‑1 (MQ=255)                                                                                      aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTgg                                                   <  2:225463/151‑1 (MQ=255)                                                                                       gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGa                                                  <  2:909356/151‑1 (MQ=255)                                                                                           tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGCGCGAAGATGAGTGGATccc                                              >  2:282070/1‑151 (MQ=255)                                                                                            gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                             <  1:319159/151‑1 (MQ=255)                                                                                             tttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATAACCACACCGGTGCCGGGCGCGAAGATGAGTGGATccc                                              >  2:1195757/1‑149 (MQ=255)                                                                                                ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATGCCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                         <  1:612204/151‑1 (MQ=255)                                                                                                   aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                                      <  2:194402/151‑1 (MQ=255)                                                                                                    ataaCCCTGGCGAAACGCGAATGAGTGGCGGAGGGGTGACTTAATATCTTTAACAAGGCCGCCAGGAAAATTATGCCCGCCTGATCATCATCGATCCCCGCTATATCGACACCCGGGCCGGGGGCGGAGATGAGTGGGTTCCCATTCGTcc                                     >  1:245250/1‑151 (MQ=255)                                                                                                      aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCg                                    <  2:1151686/150‑1 (MQ=255)                                                                                                       aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCCGAAATCTAATGCCCCCATGATCATCATCGATCCCCGCCATACCCACACCCGGGCCGGGGGCGGAAAAGAGGGGATCCCTATTCGTCCggg                                  >  2:577985/1‑151 (MQ=255)                                                                                                       aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTAATATCTTTAACAGGCACGCCAGAAATCTAATGCCCCCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGAAGATGAGTGGAACCCCATTCGTCCggg                                  >  2:10399/1‑151 (MQ=255)                                                                                                        cccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGa                                 <  1:289248/151‑1 (MQ=255)                                                                                                        cccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAAACTAATGCCCCCCTGGTCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGCGCGGAGATGAGGGGATCCCCCTTCGTCCggg                                  >  1:758679/1‑150 (MQ=255)                                                                                                         ccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGGACGCCAGAAATCTAATGGCCGCGTGATCATCATCGATTCGCGCTATACCGACACCGGGGGCCGGGGCGGAGAAGAGTGGATCCCTATTCGTcc                                     >  1:290606/1‑146 (MQ=255)                                                                                                           tGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGAACCGCGCTATACCCACACCGGGGCCGGGCGCGAAGAAGAGGGGATCCCTATTTGTCCgggga                                >  2:412010/1‑147 (MQ=255)                                                                                                               gAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGGAGATGAGTGGATCCCTATTCGTCCGGGGACa                              >  2:110766/1‑147 (MQ=255)                                                                                                                aaaCGCGAATGAGCGGCGGTGGGGTGACTTACTATCTTGATCAGGCCCGCCAGAAATCTAACGCCCGCATGATCATCATCGATCCGCGCTGTACCGCCACCGGTGCCGGGCGCGAAGATGAGTGGATCACTAGTCGTCCGGGAACAGATGc                         <  1:186478/151‑1 (MQ=255)                                                                                                                 aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGcc                        <  1:131575/151‑1 (MQ=255)                                                                                                                 aaCGCGAATGAGTGGCGGTGGGGGGGATTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCCCGCTATACCGACACCGGGGCCGGGGGCGAAGATGAGGGGGTCCCTATTTGTCCGGGGACAGATGcc                        >  2:108816/1‑151 (MQ=255)                                                                                                                    gcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCa                     <  2:111696/151‑1 (MQ=255)                                                                                                                     cgAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGGAGATGAGGGGATCCCTATTCGTCCGGGGACAGAAGCCGccc                    >  1:536479/1‑149 (MQ=255)                                                                                                                           aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa              <  2:1114105/151‑1 (MQ=255)                                                                                                                              ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg           <  2:1113575/151‑1 (MQ=255)                                                                                                                                  gTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt       <  1:1043080/151‑1 (MQ=255)                                                                                                                                   tGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg      >  1:631902/1‑151 (MQ=255)                                                                                                                                      ggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCg   <  2:658108/151‑1 (MQ=255)                                                                                                                                       gTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  1:52854/1‑151 (MQ=255)                                                                                                                                        tGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  1:1119847/1‑150 (MQ=255)                                                                                                                                        |                                                                                                                                                      CGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941578‑941862

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