Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	826,942	0	A	C	42.4%	85.9 / 13.6	59	V104G (GTG→GGG)	ybhS	putative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (24/10);  new base C (0/25);  total (24/35)
Fisher's exact test for biased strand distribution p-value = 9.62e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCATTTGCCAGATCTGCCAGATCCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCC  >  NC_000913/826796‑827089                                                                                                                                                    |                                                                                                                                                      gctGCATTTGCCAGATCTGCCAGATCCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAa‑ccacc                                                                                                                                                     >  1:163531/1‑149 (MQ=255)          gccagatctgccagatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAAt                                                                                                                                         <  2:406665/151‑1 (MQ=255)             agatctgccagatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAAtttt                                                                                                                                      >  1:1233350/1‑151 (MQ=255)                   gccagatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCg                                                                                                                                <  1:180448/150‑1 (MQ=255)                   gccagatcCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAACAA‑CCACC‑AGACCGCGAATTTT‑CCCCg                                                                                                                                <  2:863057/150‑1 (MQ=255)                          cccTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCAtt                                                                                                                        >  1:379963/1‑151 (MQ=255)                           ccTTCGACATACCCCTTTACAAAGTTAGGGGTATTCGGCCCACTGCCGTCGGTAATCCCCTGAATCGGTGCGGTGGCGTTGGCGCCCCCCATCTGTTCCGCAAAACCCACCGAAAAAC‑CCCCC‑AGCCCGCAATTTTT‑CCCCGCCTGCAttt                                                                                                                       <  1:201428/151‑1 (MQ=255)                            cTTCGACATACCCCTGTACAAAGTTAGCGGGATTCGGCTAACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCAtttt                                                                                                                      >  2:1396736/1‑151 (MQ=255)                                  cccaCCCCGGAACAAGTTTCGCGGTTTTCGTCTCACTCCCGTGGGTGATCACCTAAACCGTTGCGGGGGCGTTGCCGCGCTTCTTCTCTCCCCAAAACTCACCCGAAATAG‑CCACC‑AGCCCTAGAATTTT‑CCCCGCCTGCATTTTGGCGAt                                                                                                                <  1:380090/148‑1 (MQ=255)                                     cccccTGTACAAAGTTAGCGGTATTCGGCTCACGGCCGTCGGTGATCACCTGAATCGGTGCGGGGGCGTTGGCGCGCTCCATCTGTTCCGCAAAACCCCCCGGAATAC‑CCCCC‑AGCCCGCGAATTTC‑CCCCGCCTGCATTTTGGCGATc                                                                                                               <  2:543679/148‑1 (MQ=255)                                      ccccTTTACAAAGTTAGCGGTATTCGGCTCCCTGCCGTCGGTAATCACCTGAATCGGGGGGGTGGCGTGGGCGCGCTCCATCTGTCCCAAAAAACCCCCGGGAATAA‑CCACC‑AGACCGCAATTTTT‑CCCCGCCTGCATTTTGGCGATCAGt                                                                                                            <  2:959947/151‑1 (MQ=255)                                      ccccTGAACAAAGTTGGGGTTATTCGGCTCACGGCCGTGGGTGACCCCCTGAATCGGGGGGGGGGCTTGGGCGCCCTCCATCTGTCCCGAAAAACCCCCGGAAATAC‑CCCCC‑AGACCGCGATTTTT‑CCCCGCCTGCATTTTGGCGATCa                                                                                                              <  1:1418694/149‑1 (MQ=255)                                       cccTGTACAAAGTTAGCGGTATTCGGCCCCCGGCCGCCGTTGATCCCCTGAAGCGGTGCGGTGGCGTTGGCGCGCCCCATCTGTCCCGCAAAATCCA‑CGGAAAAA‑CCCCCAAGCCCGCGATTTTT‑CCCCGCCTGCATTTTGGCGATCAGtt                                                                                                           <  2:243694/151‑1 (MQ=255)                                           gTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTg                                                                                                       >  2:624339/1‑151 (MQ=255)                                             aCAAATTTAGCGTTATGCGGCCCCCTGCCGTCGGTGATCCCCTGAACCGGTGCGGGGGCTTGGGCGCGCTCCATCTTTCCCACAAAACCCCCCGAAATAC‑CCCCC‑AGCCCGCGATTTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGAc                                                                                                     <  1:1241821/151‑1 (MQ=255)                                             aCAAAGTTGGCGGTATTCGGCTCCCTGCCGTCGTTGATCACCTGAATGGGTGGGGTGGCGTTGGCGCGCTCCATCTGTCCCGAAAAACCCACCGAAATAA‑CCCCC‑AGACCGCGATTTTT‑CCCCGCCTGCATTTTGGCGATCAGTTcc                                                                                                         <  2:759177/147‑1 (MQ=255)                                              cAAAGTTAGCGGCATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACg                                                                                                    >  2:885479/1‑151 (MQ=255)                                               aaaGTTGGGGGTATCCGCCTCCCTCCCGTCGGTAATCACCTAAATCGGTGCGGGGGCGTTGGCGCCCTCCATCTGTCCCAAAAAATCCACCGGAATAC‑CCCCC‑AGCCCGCGAATTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGa                                                                                                      <  1:919774/148‑1 (MQ=255)                                               aaaGTTGGCGGTATTCGCCTCACGGCCTTCGGTGACCACCTGGTTCGGGGGGGTGGCGTGGGCCCCCTCCATCTGTCCCGAAAAACCCACCGAAATAA‑CCACC‑AGCCCGCGATTTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACg                                                                                                    <  1:1349032/150‑1 (MQ=255)                                               aaaGTTAGCGTTTTTCGGCTCACTGCGGTCGGTGATCCCCTGAAGCGGGGCGGGGGGGTGGGCGCGCTCAATCTGTCCCGAAAAACCCCCCGGAATAC‑CCACC‑AGCCCGCGATTTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGa                                                                                                      <  2:1334022/148‑1 (MQ=255)                                               aaaGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACgg                                                                                                   >  1:1452717/1‑151 (MQ=255)                                               aaaGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACgg                                                                                                   >  1:475587/1‑151 (MQ=255)                                               aaaGTTAGCGGGATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGAGATCAGTTCCTGACgg                                                                                                   >  2:408909/1‑151 (MQ=255)                                                   ttAGCGTTATTCGGCTCCCTGCCGTCGGTGATCCCCTGAATCGGGGCGGTGGCGTTGGCGCGCTCCATCTGTCCCGCAAAACCCCCCGAAATAC‑CCCCC‑AGACCGCAATTTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTa                                                                                                <  2:569565/150‑1 (MQ=255)                                                   ttAGCGGTTTTCGGCTCCCTGCCGTGGGTGATCCCCGGAATCGGTGGGGGGGGGTTGGCGCGCCCCATCTGTCCCGAAAAACCCCCCGAAATAC‑CCCCC‑AGCCCGCGAATTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGtt                                                                                                 <  2:583115/149‑1 (MQ=255)                                                     aGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCACGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCg                                                                                             >  1:927723/1‑151 (MQ=255)                                                        gtttttCGGCTCCCTGCCGTCGGTGACCACCGAAATCGGTGCGGGGGGGTTGGCGCGCTCCATCTGTCCCGAAAAACCCCCCGAAATAC‑CCCCC‑AGCCCGCGATTTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTAt                                                                                               <  2:1162462/142‑1 (MQ=255)                                                        ggTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTg                                                                                          >  2:822358/1‑151 (MQ=255)                                                          tATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAt                                                                                        >  2:1312208/1‑151 (MQ=255)                                                               gcctcccggCCGTGGGTAATCACCTGATTGGGGGCGGTGGCGTGGGCGCGCCCCATCTGTCCCGAAAAACCCCCCGAAAAAC‑CCACC‑AGCCCGCGATTTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAtggtgg                                                                                   <  2:1233385/143‑1 (MQ=255)                                                                gCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGc                                                                                  >  1:1400934/1‑151 (MQ=255)                                                                 ctcccggCCGTGGGTGACCCCCTGATGGGGGGGGGTGGCGTGGGCGCCCCCCATCT‑TTCCCGAAAAACCCCCGGAAAAACCCCCC‑AGCCCGCAATTTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAtg                                                                                       <  2:766285/139‑1 (MQ=255)                                                                 ctcccgcccGTCGGTAACCCCCTGAATCGGTGGGGGGGCGTGGCCGCCCTCCTTTTGTCCCGAAAAATCCCCGGAAATAC‑CCCCC‑AGCCCGCAATTTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAtggtgg                                                                                   <  1:725146/142‑1 (MQ=255)                                                                 ctccccgCCGTGGGTGACCCCCTAAACGGGGGCGGTGGCGTTGGCGCGCTCCATCTTTCCCAAAAAACCCACCGAAATAC‑CCCCC‑AGCCCGCGATTTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGAtggtg                                                                                    <  1:499980/142‑1 (MQ=255)                                                                  tCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGt                                                                                >  2:670327/1‑151 (MQ=255)                                                                   cACTGCCGTCGGTGATCCCCTGAATCGGTGCGGTGGCGTTGGCGCCCTCCTTCTGTCCCGCAAACCCCCCGGAAATAA‑CCCCC‑AGCCCGCAATTTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTc                                                                               <  1:145197/151‑1 (MQ=255)                                                                    aCTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCg                                                                              >  1:307350/1‑151 (MQ=255)                                                                         cGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGt                                                                          >  1:442996/1‑150 (MQ=255)                                                                             gttGATCCCCGGAACCGGTGCGGTGGCGTTGGCCCGCTCCATTTGTCCCCAAAAACCCACGGGAATAC‑CCCCC‑AGCCCGCAATTTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAg                                                                        <  1:1040549/146‑1 (MQ=255)                                                                             ggTGATCCCCTGAATCGGTGCGGTGGCGTGGGCGCGCCCCATCTTTCCCGAAAAACCCCCCGAAATAA‑CCCCC‑AGCCCGCAATTTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGc                                                                     <  1:1437437/151‑1 (MQ=255)                                                                             ggTGATCCCCTGAACCGGGGGGGGGCCGTGGGCGCGCCCCTTCTTTCCCGAAAAACCCCCCGAAATAC‑CCCCC‑AGCCCGCGATTTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTa                                                                         <  1:768764/147‑1 (MQ=255)                                                                                   ccccTGAATCGGTGCGGTGGCGGTGGCGCGCTCCATCTTTCCCGAAAAATCCACCGGAATAC‑CCCCC‑AGCCCGAGATTTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAg                                                                        <  2:475600/140‑1 (MQ=255)                                                                                   ccccTGAACCGGTGCGGTGGCGTTGGCGCGCTCCTTCTGTTCCGCAAATCCCCCGGGAAAAC‑CCCCC‑AGCCCGCAATTTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAAcc                                                                <  1:574739/148‑1 (MQ=255)                                                                                    aCCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCgg                                                              >  2:283795/1‑151 (MQ=255)                                                                                     cctaattggttgcggtgGCGTGGGCGCGCTCCATCTGTTCCGAAAAACCCACCGAAATAC‑CCCCC‑AGCCCGCAATTTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGa                                                                   <  2:1218049/135‑1 (MQ=255)                                                                                      cTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGAGAACCGGt                                                             >  1:1094545/1‑150 (MQ=255)                                                                                         aaTCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATg                                                         >  2:865455/1‑151 (MQ=255)                                                                                                   gtggcggtggcgCGCTCCATTTGTCCCGCAAACTCCCCCGAAATAC‑CCCCC‑AGCCCGCTATTTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGa                                               <  1:381979/151‑1 (MQ=255)                                                                                                           ggCGCGCTCCTTCTG‑TCCCCAAAAACCCCCGGAAAAACCCCCC‑AGCCCG‑GAATTTTCCCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAAt                                            <  1:665244/146‑1 (MQ=255)                                                                                                            gcgcgcTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAgcg                                      >  1:1157747/1‑151 (MQ=255)                                                                                                                  tCCATCTGTTCCGCAAAATCCATCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCtt                                >  1:737075/1‑151 (MQ=255)                                                                                                                       cTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGCGAAATCCAGCGCCGCTTCGc                             >  2:648177/1‑149 (MQ=255)                                                                                                                        tGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTAc                          >  2:53911/1‑151 (MQ=255)                                                                                                                            ccccaaacccccccGAAATAC‑CCCCC‑GGCCCGGGATTTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCg                              <  2:268908/131‑1 (MQ=255)                                                                                                                                 aaaCCCACGGGAAAAC‑CCACC‑CCCCCGGGAATTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCa                 <  1:664647/151‑1 (MQ=255)                                                                                                                                     cccccGAAATAA‑CCCCC‑GGCCCGGGAATTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGt               <  1:421230/146‑1 (MQ=255)                                                                                                                                          ggAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATcc        >  2:337916/1‑151 (MQ=255)                                                                                                                                          gaaaTAT‑CCACC‑CGACCGGTAATTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTaaaa           <  1:207020/146‑1 (MQ=255)                                                                                                                                           gAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATccc       >  2:807096/1‑151 (MQ=255)                                                                                                                                            aaaaa‑CCCCC‑AGCCCGCAATTTTC‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCg      <  2:45282/148‑1 (MQ=255)                                                                                                                                                a‑ccacc‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGAccc  >  1:405914/1‑151 (MQ=255)                                                                                                                                                a‑ccacc‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGAccc  >  2:70338/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GCTGCATTTGCCAGATCTGCCAGATCCCTTCGACATACCCCTGTACAAAGTTAGCGGTATTCGGCTCACTGCCGTCGGTGATCACCTGAATCGGTGCGGTGGCGTTGGCGCGCTCCATCTGTTCCGCAAAATCCACCGGAATAA‑CCACC‑AGACCGCGAATTTT‑CCCCGCCTGCATTTTGGCGATCAGTTCCTGACGGTTATCGCTGATGGTGGCGTCGATGTAGGGCGAACCGGTCATGGTGTGGGTGAAATCCAGCGCCGCTTCGCTACGCTGTTCCAGTAAAATCCCGACCC  >  NC_000913/826796‑827089

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