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

TACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGT  >  NC_000913/941568‑941855                                                                                                                                                  |                                                                                                                                               tactCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                           >  1:1324491/1‑151 (MQ=255)  actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTa                                                                                                                                          <  1:1099886/151‑1 (MQ=255)       tccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTa                                                                                                                                       >  2:1111502/1‑149 (MQ=255)        ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCt                                                                                                                                    <  1:37772/151‑1 (MQ=255)               aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCCTTCTATCTTTaaa                                                                                                                               >  2:1103027/1‑148 (MQ=255)                 ttGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTaatgtttttt                                                                                                                                 >  2:859521/1‑136 (MQ=255)                  tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGc                                                                                                                           <  2:659587/150‑1 (MQ=255)                  tGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTATCTTTAACAgg                                                                                                                            >  1:1260427/1‑149 (MQ=255)                   gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCAc                                                                                                                         <  1:813293/151‑1 (MQ=255)                     ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGc                                                                                                                       <  2:127823/151‑1 (MQ=255)                       aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCa                                                                                                                     <  1:246653/151‑1 (MQ=255)                          gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCTTATTATCTTTAACAGGCACGCCAGa                                                                                                                   >  1:446071/1‑150 (MQ=255)                             tGAACTATACCTACGGCGGCTGGGCAGATGGCCACAGCCCGTCGGATATCGCAAACAGTAAGCTGGTAGGGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTTCTAATTTTAACACACACGCCCGaaa                                                                                                                 >  2:408131/1‑149 (MQ=255)                             tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTTTCTTGTAAAAGCACGCCCGaaa                                                                                                                 >  1:837120/1‑149 (MQ=255)                               aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCCTACTATCTTGGAAAAGGACCCCAGAAAACt                                                                                                              >  1:227252/1‑150 (MQ=255)                                aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGTTATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATTTTGAACAAGCACGCCAGaa                                                                                                                  >  2:73151/1‑145 (MQ=255)                                aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCAACCCAGaaa                                                                                                                 >  1:302327/1‑146 (MQ=255)                                     aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTAACAGGGACGCCAGAAAACTAATGtcc                                                                                                       >  2:62490/1‑151 (MQ=255)                                      ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCg                                                                                                      <  1:327184/151‑1 (MQ=255)                                         aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                   <  1:1317973/151‑1 (MQ=255)                                           ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAAGCACCCCAGAAATCTAATGCCCGc                                                                                                     >  2:1068771/1‑147 (MQ=255)                                            gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGat                                                                                                <  2:885217/151‑1 (MQ=255)                                              ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                              <  2:903222/151‑1 (MQ=255)                                                  gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                           <  1:449552/150‑1 (MQ=255)                                                  gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATTTTTAAAAAGCACGCCAGAAAACTAATGCCCCCACGatca                                                                                              >  2:988241/1‑147 (MQ=255)                                                  gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGatcatc                                                                                            >  1:570093/1‑149 (MQ=255)                                                    gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                        <  1:778775/151‑1 (MQ=255)                                                     cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTTAACAGGCACGCCCGAAATCTAATGCCCGCCTGaacatca                                                                                           >  1:1231894/1‑147 (MQ=255)                                                      aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTTAACAGGCACGCCAGAAAACTAAAGCCCCCCAGATTATCATCggt                                                                                      >  1:874637/1‑149 (MQ=255)                                                       gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc                                                                                     <  2:885069/151‑1 (MQ=255)                                                       gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCAAGCCAGAAAACTAATGCCCGCATTGTCATCATCGAt                                                                                      >  2:323885/1‑150 (MQ=255)                                                       gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTACTATCTTTAACAGGGACGCCAGAAAACTAAAGCCCGCATGATCCACATCGAt                                                                                      >  2:904375/1‑150 (MQ=255)                                                                         gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCg                                                                   <  1:1253477/151‑1 (MQ=255)                                                                            atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATCTTTAACAAGCACGCCCGAAAACTAAAGCCCGCCCGATTATCATCGATACGCGCTATACCGACAACGGTg                                                                >  2:465795/1‑151 (MQ=255)                                                                              cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGcc                                                              <  1:370138/151‑1 (MQ=255)                                                                              cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGc                                                               <  2:733866/150‑1 (MQ=255)                                                                                  aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATATAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCggg                                                           >  1:1111298/1‑150 (MQ=255)                                                                                  aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCCTACTTTCTTTAACAGGCACGCCAGAAATCTAAAGCCCCCCTGAATATCATCGATACGCGCCCTACCTACACCgggg                                                                >  2:1047740/1‑143 (MQ=255)                                                                                             ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATg                                                 <  2:1135775/149‑1 (MQ=255)                                                                                              gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTAACAGGCCCGCCAGAAAATTAATGCCCCCCTGAACATCATCAATCCGCGCTATACCGGCACCGGTGGCGGGGGGGaa                                                     >  2:1072122/1‑144 (MQ=255)                                                                                                aGTGCTGGTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTTCTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTg                                             <  1:73234/150‑1 (MQ=255)                                                                                                    cTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCTGGGGGGCTTACTATCTTGAAAAAGCACGCCCGAAAACTTATGCCCGGACTAAAATCATTGATCCCCGCCATTACGCCACCGGGGGCGGGCGCGCAGAAGAGGGGTTcc                                        >  2:937028/1‑151 (MQ=255)                                                                                                      gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATccc                                       <  1:703757/150‑1 (MQ=255)                                                                                                         tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCCGAAAACTAATGCCCCCATGATCATCATCGATCCGCGCCATACCCACACCCGGGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                   >  1:1298410/1‑151 (MQ=255)                                                                                                           gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCg                                 <  2:1324461/151‑1 (MQ=255)                                                                                                             aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                               <  2:883664/151‑1 (MQ=255)                                                                                                             aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                               <  2:438887/151‑1 (MQ=255)                                                                                                              ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTc                               <  1:820462/150‑1 (MQ=255)                                                                                                                   ccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAATCTAATTCCCGCATGATCATCATCGATCCCCGCTATACCCACAACCGTGCCGGGCGGGGAGAAGGGGGGGTCCCTATTTGTCCgg                            >  2:1020238/1‑148 (MQ=255)                                                                                                                      ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                        <  2:574550/149‑1 (MQ=255)                                                                                                                      ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAg                      <  2:1260401/151‑1 (MQ=255)                                                                                                                         gAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGGACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCCGGGCCGGGCGCGAAGATGAGGGGATCCCTATTTGTCCggggaaag                      >  2:736841/1‑143 (MQ=255)                                                                                                                           aaCGCGAATGAGTGGCGGTGGGGGGACTTAATATCTTTAACAGGCACGCCAGAAAACTAAAGCCCCCATGATCATCATCGATCCGCGCCATACCGACACCGGGGCCGGGGGCGGAGATGAGGGGAACCCTATTCGTCCGGGGACAGAAGcc                 >  2:472398/1‑151 (MQ=255)                                                                                                                            aCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCg                <  1:1068792/151‑1 (MQ=255)                                                                                                                             cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGc               <  2:598324/151‑1 (MQ=255)                                                                                                                                    gAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGtt        <  1:904961/151‑1 (MQ=255)                                                                                                                                     aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa       <  2:4342/151‑1 (MQ=255)                                                                                                                                     aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa       <  1:1048499/151‑1 (MQ=255)                                                                                                                                        ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc     <  1:12035/150‑1 (MQ=255)                                                                                                                                        ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg    <  1:1010269/151‑1 (MQ=255)                                                                                                                                         gCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACgg   >  2:367429/1‑151 (MQ=255)                                                                                                                                          cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt  >  1:543643/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                               TACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGT  >  NC_000913/941568‑941855

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