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

ACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTC  >  NC_000913/941566‑941856                                                                                                                                                    |                                                                                                                                                 actactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGGTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGAc                                                                                                                                               >  2:1205895/1‑151 (MQ=255)    actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                                                                                              <  1:656063/149‑1 (MQ=255)    actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                             >  2:195669/1‑150 (MQ=255)       cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTAc                                                                                                                                           <  1:660085/149‑1 (MQ=255)       cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTa                                                                                                                                         <  2:801544/151‑1 (MQ=255)             cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTAttttt                                                                                                                                    >  1:354911/1‑149 (MQ=255)               cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAc                                                                                                                                 <  1:485711/151‑1 (MQ=255)                   ttGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGGTACTATTTTGaaaa                                                                                                                                >  1:566828/1‑146 (MQ=255)                        gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGcc                                                                                                                        <  2:777458/151‑1 (MQ=255)                           ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGa                                                                                                                     <  2:492032/151‑1 (MQ=255)                           ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGa                                                                                                                     <  1:503086/151‑1 (MQ=255)                           ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGa                                                                                                                     <  2:1039069/151‑1 (MQ=255)                            gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaa                                                                                                                    <  2:317226/151‑1 (MQ=255)                            gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaa                                                                                                                    <  1:406067/151‑1 (MQ=255)                             tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTCTCTTTTACAAGCCCGcc                                                                                                                        >  1:740396/1‑146 (MQ=255)                                gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATc                                                                                                                 <  1:625880/150‑1 (MQ=255)                                gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATc                                                                                                                 <  1:158600/150‑1 (MQ=255)                                gAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACCATCTTTAACAAGCACGCCAGAAATc                                                                                                                 >  1:513293/1‑150 (MQ=255)                                        ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGGTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCg                                                                                                       <  1:1205835/151‑1 (MQ=255)                                          tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTTTATATTTTAAAAGCACGCCAGaaaaataaaaccccc                                                                                                       >  1:377993/1‑138 (MQ=255)                                           aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCa                                                                                                     <  1:439038/150‑1 (MQ=255)                                             ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCAGAAAACTAATG‑CCCGc                                                                                                      >  2:765508/1‑147 (MQ=255)                                             ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGATCTAACATCTTTAAGAAGCACGCCAGAAAACTAATA‑GCCGc                                                                                                      >  2:454726/1‑147 (MQ=255)                                             ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATCTTGAAAACGGACGCCCGAAAACTTATG‑ccc                                                                                                        >  1:927424/1‑145 (MQ=255)                                                ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTTGGGGTACTTAATTTTTTTAACCCGGATCTCAGAAATCTAATG‑ccc                                                                                                        >  2:362652/1‑142 (MQ=255)                                                ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGatca                                                                                               <  1:682255/151‑1 (MQ=255)                                                ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTTTCTTGAACAAGCACGCCCGAAAACTAATG‑cccccct                                                                                                    >  2:393486/1‑142 (MQ=255)                                                ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCATATTATCTTTAACAGGGAGGCCCAAAAAATAAAG‑CCCGcgaaa                                                                                                  >  2:51412/1‑144 (MQ=255)                                                     ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGatcatcatc                                                                                          <  1:819905/151‑1 (MQ=255)                                                      gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCg                                                                                         <  2:675091/151‑1 (MQ=255)                                                        aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGAt                                                                                       <  2:157395/151‑1 (MQ=255)                                                         gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATc                                                                                      <  1:744758/151‑1 (MQ=255)                                                                    cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATAcc                                                                           <  1:509277/151‑1 (MQ=255)                                                                     ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGCCCGCCAGAAAACTAATG‑CCCGCATGATCATCATCGATCCCCGCTATAcc                                                                           >  1:333425/1‑150 (MQ=255)                                                                      cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCg                                                                          <  1:695691/150‑1 (MQ=255)                                                                        tCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATTTTGAACAGGCACGCCAGAAAACTAAAG‑CCCGCCTGATCATCATCGATCCGCGCTATACCCac                                                                        >  2:250871/1‑150 (MQ=255)                                                                          ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACAcc                                                                     <  1:325562/151‑1 (MQ=255)                                                                               tCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCCATCTTTAACAGGGACGCCAGAAAACTAATG‑CCCGCATGATCATCATCGATCCCCGCTATACCCACACCgggg                                                                 >  2:883354/1‑148 (MQ=255)                                                                                cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGcc                                                               <  2:776608/151‑1 (MQ=255)                                                                                cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGc                                                                <  1:71852/150‑1 (MQ=255)                                                                                cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGc                                                                <  2:83663/150‑1 (MQ=255)                                                                                         tAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTTAACAAGCACGCCACGAAACTAAAG‑GCCGCCTGATCATCATCGAACCGCGCCATACCCACAAccgggccggg                                                            >  2:1292554/1‑145 (MQ=255)                                                                                          aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGAATTTATATCCTGAACAAGGCCGGCCGAAATAGAAAA‑GCCGGAAGAGTATCATCAATCCGCGCCATAACGAAACCCGGGCCgcg                                                            >  1:428360/1‑142 (MQ=255)                                                                                            gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                   <  2:66704/151‑1 (MQ=255)                                                                                            gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCCCGCCCGAAAACTAATG‑CCCCCATGATCATCCTCCATCCCCGCTATACCGACACCCGGGCCGGGCGCGGAg                                                     >  2:684157/1‑149 (MQ=255)                                                                                               ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCCTAATATCTTTTACAGGGACCCCAGAAAACTTATT‑ACCCCCAAATTCTCATCGGACCGCGCCATAAACAAACCGGGGCCGGGcgcgga                                                      >  2:1100638/1‑143 (MQ=255)                                                                                               ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTTAACAAGCACCCCAGAAATCTAAAG‑CCCGCAAGATCATCCACGATCCGCGCTATACCGACACCCGTGCCggg                                                            >  2:863928/1‑139 (MQ=255)                                                                                                   gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACGCCCGAAATCTAATG‑CCCGCATGATCATCATCGATCCCCGGTCTACCCCCAACGGGGCCGGGGGCGAAGATGAg                                                >  2:1148309/1‑147 (MQ=255)                                                                                                     gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGAt                                           <  1:454716/150‑1 (MQ=255)                                                                                                     gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGAt                                           <  2:369216/150‑1 (MQ=255)                                                                                                           tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGTGTTCTATCTTGAACAGGGACGCCAGAAATCTAATG‑CCCGGCAGATAAAACTCGAGCCGCGATTTACCCACACCCGTGGCGGGCGTGGAAATGAggag                                             >  1:485590/1‑139 (MQ=255)                                                                                                             gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATCAGTGGATCCCTATTCg                                  <  1:195668/151‑1 (MQ=255)                                                                                                                  aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTTAAAAAGCACGCCCGAAAACTAATG‑GCCCCATGATCCTCATCTATTCCCGCCATACCGACACCGGGGCCGGGGGCGGAGATGAGGGGATCCCCATTCCTCCgg                             >  2:1140850/1‑151 (MQ=255)                                                                                                                  aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAAGCACGCCAGAAAACTAAAG‑ACCCCCTGATCATCATCCGTCCCCGCCATAACCACACCGGGGCCGGGGGCCAAGGTGAGTGGATCCCTATTCGTcccg                             >  2:146777/1‑149 (MQ=255)                                                                                                                  aaCCCTGGCGAAACGCGAATGAGTGGCGGACGGGGGGCTTAATATTTTGAAAAGGGCCGCCAGAAATATAATG‑GCCGCATGGTAATCAACCAGCCTCGCTATAACGTCACCCGGGCCGCGGGCGCAGATGAGGGGGCGCCTTTTCGTCCgg                             >  1:226986/1‑151 (MQ=255)                                                                                                                   aCCCTGGCGAAACGCGAATGAGTGGCGGTTGGGGTAGTTAGTATCCTGAAACGACACTCCCCAAAACTAAAGCCCCGCA‑GAACATCCCTCCTACACCCTATACACACACCCGGGCCGGGCGCGGAGATGATGGGGACCCTAtttgctccgg                            >  2:67900/1‑143 (MQ=255)                                                                                                                     ccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGaa                          <  2:130780/151‑1 (MQ=255)                                                                                                                       tGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAAACTAATG‑CCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGGGCGGAGATGAGTGGGTCCCTATTTGTCCggg                            >  2:188671/1‑147 (MQ=255)                                                                                                                        ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                         <  2:912425/149‑1 (MQ=255)                                                                                                                          cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                     <  2:652901/151‑1 (MQ=255)                                                                                                                                gcgAATGAGTGGCGGTGGGGGGGCTTACTATTTTTAACAGGCACGCCAGAAAACTAATG‑GCCGCCTGATCATCATCGATTCGCCCTATACCGACACCGGTGCCGGGGGCGAAGATGAGTGGGTCCCTATTCGTCCGGGAACAGAAGCCGc                >  1:1132570/1‑150 (MQ=255)                                                                                                                                  gAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACt             <  1:766455/151‑1 (MQ=255)                                                                                                                                   aaTGAGTGGCGGTGGGGGGACTTACTTTCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCCTGATCATCATCGATCCGCGCTATACCCACACCCGTGCCGGGCGCGAAGATGAGTGGGTCCCTATTTGTCCGGGAACAGATACCGCACTg            >  2:469375/1‑151 (MQ=255)                                                                                                                                    aTGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAATCTAATG‑CCCCCATGATCATCATCGATCCGCGCTATACCGACACCAGGGCCGGGCGCGCAGATGAGTGGGTCCCTATTCGTCCCGGGACAGATGCCGCACTgg           >  2:1244261/1‑151 (MQ=255)                                                                                                                                       aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa        <  2:1160591/151‑1 (MQ=255)                                                                                                                                           gCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACgg    >  2:61023/1‑151 (MQ=255)                                                                                                                                             ggTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc  >  2:932188/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                 ACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTC  >  NC_000913/941566‑941856

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