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

GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941563‑941860                                                                                                                                                       |                                                                                                                                                     gCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                                                                                      <  1:834217/151‑1 (MQ=255)        ctcctcCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTa                                                                                                                                                <  1:161527/150‑1 (MQ=255)         tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTgccc                                                                                                                                                  >  1:1103584/1‑144 (MQ=255)             ccGCGCAAAGTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATc                                                                                                                                           <  2:854238/150‑1 (MQ=255)              cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑tt                                                                                                                                        <  2:177407/151‑1 (MQ=255)              cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑tt                                                                                                                                        <  1:336768/151‑1 (MQ=255)               gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTg                                                                                                                                       <  1:530065/151‑1 (MQ=255)               gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTATTATC‑tt                                                                                                                                        >  2:975761/1‑150 (MQ=255)                cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTATC‑tt                                                                                                                                        >  2:90906/1‑149 (MQ=255)                 gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGaa                                                                                                                                     <  2:45804/151‑1 (MQ=255)                        gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCAc                                                                                                                              <  2:824042/151‑1 (MQ=255)                        gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTTAACAAGCa                                                                                                                               >  1:849804/1‑150 (MQ=255)                         cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACg                                                                                                                             <  2:90441/151‑1 (MQ=255)                             aGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTCTTTC‑TTGGACAAGGACGcc                                                                                                                           >  2:853782/1‑149 (MQ=255)                                  tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAAt                                                                                                                     <  1:1169829/150‑1 (MQ=255)                                  tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTTAACAGGGACGCCAaaaa                                                                                                                      >  1:467835/1‑149 (MQ=255)                                    aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCt                                                                                                                   <  2:1006224/150‑1 (MQ=255)                                      cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACCTACTCTC‑TTTAAAAGGCCCGCCAGAAATCTaa                                                                                                                 >  1:167463/1‑150 (MQ=255)                                      cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATT‑TTTAACAAGGACGCCAGAAAACTaa                                                                                                                 >  1:443259/1‑150 (MQ=255)                                        ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAAt                                                                                                                <  2:236365/149‑1 (MQ=255)                                        ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAAt                                                                                                                <  2:252177/149‑1 (MQ=255)                                           ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCg                                                                                                           <  1:593211/151‑1 (MQ=255)                                               cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATg                                                                                                       <  1:9594/151‑1 (MQ=255)                                                ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGa                                                                                                      <  1:1116490/151‑1 (MQ=255)                                                  cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc                                                                                                    <  1:420649/151‑1 (MQ=255)                                                   ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCGTACTTTC‑TTTTAAAGGGGAGCCCGAAAAATAATACCCGCGTGaa                                                                                                     >  2:313478/1‑148 (MQ=255)                                                      tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTGAACAGGCACGCCCGAAAACTAATGccccc                                                                                                          >  2:88323/1‑139 (MQ=255)                                                       gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTTTGGGGGGTTACTATC‑TTTAAGAATGACGCCACGAAACTAATGACCggaagggcatc                                                                                                 >  2:1180870/1‑139 (MQ=255)                                                          cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGAATTATTATT‑TTGAACAGGCACGCCAGAAATCTAAAGGCCGCATGaatatca                                                                                                >  1:1016357/1‑147 (MQ=255)                                                            gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATT‑TTTAACAAGGACGCCAGAAATATAATGCCCGCATGATCATCATCGAt                                                                                           >  2:731536/1‑150 (MQ=255)                                                                  aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCt                                                                                    <  1:841336/151‑1 (MQ=255)                                                                  aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTTTATT‑TTGAAAAGGCAGGCCAGAAATCTCAAGCCCGCATGCACATCACCggttc                                                                                         >  1:1047739/1‑142 (MQ=255)                                                                   aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCta                                                                                   <  2:1123080/151‑1 (MQ=255)                                                                    cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATA‑TTTAACAGGCACGCCACAAAAATAATGCCCCCATGATCATAATCTGTCCccgc                                                                                     >  2:386653/1‑148 (MQ=255)                                                                        ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAACAGGCACGCCCGAAATCTAAAGCCCGCATGATCATCATCGGTCCGCGCTATAcc                                                                               >  1:629512/1‑150 (MQ=255)                                                                          gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGac                                                                            <  1:45711/151‑1 (MQ=255)                                                                          gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATA‑TTTAAAAAGCACGCCAGAAAACTAATGCCCGCATGATCATCCTTCATCCGCGCTATAcc                                                                               >  2:569674/1‑148 (MQ=255)                                                                                 atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTg                                                                     <  1:952113/151‑1 (MQ=255)                                                                                 atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTg                                                                     <  2:37965/151‑1 (MQ=255)                                                                                    gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCg                                                                  <  2:583119/151‑1 (MQ=255)                                                                                    gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATC‑TTTAAAAGGCACGCCAGAAAACTAATTCCCCCATGATCATCATCGATCCGCGCTATACCGACACCGGGGcc                                                                   >  2:258516/1‑150 (MQ=255)                                                                                    gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGTGGGTTACTATA‑TTTAAAAAGCAAGCCAGAAAACTAATTCTCCCCCGGTCATAATCGAGCCACGCTATACCCACCCCggggac                                                                   >  2:1123769/1‑146 (MQ=255)                                                                                                 tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTATC‑TTGAAAAAGCACGCCAGAAAACTTATGCCCGCCTAATTTTAATCCATCCGCGCTATATCCAGACCGGGGCCGGGGGCGGAAATg                                                      >  1:578273/1‑150 (MQ=255)                                                                                                 tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGAGGATTACTATCTTTGAACCAGCA‑GCCGGAAAACTAAAGCCCCCAAGGCCCTCATCGATCCCCGCTATACCGACACCGGGGGCCGGCGGGCAGATg                                                      >  1:603013/1‑150 (MQ=255)                                                                                                      gTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTATATC‑TTGAACAGGCACGCCCGAAATCTAAAGCCCGCCTGATCATCATTTATCCGCGCCATACCGACACCCGTGCCGGGCGGGAAGATGAGTgg                                                 >  2:601383/1‑150 (MQ=255)                                                                                                       tGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTACTATC‑TTTAAAAGGCACGCCCACAATCCCAAGACCGCAACATCAACAACGATTCCGGCTATAACGGCACCGGCGCCGGGCGCGGAagtgagtg                                                  >  2:257304/1‑148 (MQ=255)                                                                                                           gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATcc                                             <  1:975812/149‑1 (MQ=255)                                                                                                           gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAAAAAGGACGCCAGAAATCTAAAGCCCGCCTGATCATCATCGATCCCCGCTATACCGACACCGGGGGCGGGCGCGGAGATGAGGGGATccc                                            >  1:850003/1‑150 (MQ=255)                                                                                                              tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                         <  1:1089127/150‑1 (MQ=255)                                                                                                               ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGGGTGGATCCCTATTc                                       <  2:254833/151‑1 (MQ=255)                                                                                                                     aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACCTACTATC‑TTGAAAAGGCACGCCAGAAATCTAATGCCCGCATGAACATCATTGATCCGCGCTATACCCACACCCGTGCCGGGCGCGGAGATGAGTGGGTCCCTATTccgcc                                   >  2:528184/1‑145 (MQ=255)                                                                                                                      aCCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGACTTACTATT‑TTGAACAGGCACGCCCGAAAATTAATGCCCGCATGGACATTATCGATCCTCGCTATAACCACACCGGTGGCGGGGGCGAAGATGAGTGGATCCCTCTTCGTCCggg                                >  1:9822/1‑151 (MQ=255)                                                                                                                           ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAg                           <  2:767778/151‑1 (MQ=255)                                                                                                                           ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAg                           <  2:1103532/151‑1 (MQ=255)                                                                                                                           ggCGAAACGCGAATGAGTGGCGGAGGGGGGAGTTTCTATC‑CTTAACAAGCGCCCCAGAAAATTAATGCCCGCAAGATTATCATCCATCCGCCCTTTAACGCCAACCGGGGCGGGGGCGGAGATGAATGGGACCCTACTTCTCCggggaa                             >  2:759023/1‑146 (MQ=255)                                                                                                                            gCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTATTATT‑TTGAAAAAGCACGCCAGAAAATTAATGGCCGCCTGATCATCATCGATCCCCGCTATACCGACACCGGGGCCGGGGGCGGAGATGAGTGGATCCCCATTCTTCCgggga                              >  2:233769/1‑145 (MQ=255)                                                                                                                             cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                         <  1:88324/151‑1 (MQ=255)                                                                                                                               aaaCGCGAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAACACGGCCGCCAGAAATCTAATTCCCGCAAGATCATCATCGATCCGCCCCATACCTCCACCCGTGCCGGGCGCGGAGATGAGTGGATCCCCCTTCTTCCGGGAACa                            >  1:374994/1‑146 (MQ=255)                                                                                                                                  cgcgAATGAGTGGCGGTGGGGGGACTTACTATC‑TTTAAAAAGCACGCCCGAAATCTAATGCCCGCATGATCATCCTCGATCCGCCCCATACCCACACCCGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGc                    >  2:299976/1‑151 (MQ=255)                                                                                                                                             ggCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg         >  2:317427/1‑151 (MQ=255)                                                                                                                                             ggCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg         >  1:143241/1‑151 (MQ=255)                                                                                                                                                    gggTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc  >  2:710043/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                     GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATC‑TTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941563‑941860

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