Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,575	0	A	C	49.2%	85.1 / 17.6	67	V555G (GTG→GGG)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (33/1);  new base C (0/33);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 2.46e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCG  >  NC_000913/501428‑501719                                                                                                                                                      |                                                                                                                                                 aaGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTcaccac                                                                                                                                                 >  1:797016/1‑149 (MQ=255)   gTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCATCATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCTTCCCCACCTCa                                                                                                                                             >  2:820865/1‑151 (MQ=255)             gttggtggcccAGATAGGGCGTCCCCGCCGCATCTGCCTTTTGGCTCTCGATGCCTGAGGCGACGCTGGCGCGTCTTATCATGCCTCCATATTTTTCATATTTTCCATCCGGCAA‑CCCCCTTTT‑CCCCCGTCCCCCCCTCCCCCGCCGGt                                                                                                                                    <  2:803778/141‑1 (MQ=255)             gTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTg                                                                                                                                   >  2:117440/1‑151 (MQ=255)               tCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGACACCACCTCACCCGCCGGTgg                                                                                                                                  >  2:452343/1‑150 (MQ=255)                         gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCccgca                                                                                                                       >  2:940437/1‑151 (MQ=255)                         gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                                       >  2:561995/1‑151 (MQ=255)                         gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagca                                                                                                                       >  1:448337/1‑151 (MQ=255)                              ggCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTc                                                                                                                   >  1:676203/1‑150 (MQ=255)                                 gtccccgccgcTTCGGGTTTTTGGCTCTCGATGCTTGATGCGCCGCGGGCGCGTTTTATCATCCCTACATTTTTTTCATTTTTTACATCCGGCAC‑CCCCCGTTT‑CCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                                                <  1:478710/145‑1 (MQ=255)                                      cgccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                          >  2:667466/1‑151 (MQ=255)                                      cgccgcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                          >  1:188777/1‑151 (MQ=255)                                          gcATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACgg                                                                                                      >  2:160726/1‑151 (MQ=255)                                                ggCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt                                                                                                >  1:758573/1‑151 (MQ=255)                                                ggCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt                                                                                                >  1:260275/1‑151 (MQ=255)                                                  cttttGCCTCTCGATGCCTGATGCGACGCTGGCCCTTCTTAACAGCCCTACATTTTTTTCATATTTTACATCCGGCAA‑CCCCCTTTT‑CCCCCCCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                                                                              <  1:172683/149‑1 (MQ=255)                                                   ttttGCCTCTCGTTGCCTGATGCGCCGCGGGCGCGTCTTATCATCCCTACATTTTTTTAATATTTTACATCCGGCAC‑CCACCTTTT‑CCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc                                                                                              <  2:648759/149‑1 (MQ=255)                                                      tgtctctGGATGCCGGTGGGGCCGCGGGCGCGTCTTATCTTCCCTACATTTTTTTCATTTTTTACTCCCGGAAC‑CCCCCTTTT‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGc                                                                                          <  2:835375/148‑1 (MQ=255)                                                       ggCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACg                                                                                           >  1:1041234/1‑149 (MQ=255)                                                        gCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                                         >  2:25226/1‑150 (MQ=255)                                                           ctcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                                      >  1:842770/1‑150 (MQ=255)                                                           ctcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                                      >  1:1134267/1‑150 (MQ=255)                                                                 gCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCa                                                                                >  2:242764/1‑150 (MQ=255)                                                                  cccgaggcgcCGCTGGCGCGTCTTCTCATCCCCCCTTTTTTTTAATATTTCCCCCCCGGCAC‑CCCCCTTTT‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                              <  1:824348/145‑1 (MQ=255)                                                                  ccTGATGCGCCGCGGGCGCTTTTTTTCATCCCTACATTTTTTTCATTTTTTCCACCCGGCAA‑CCCCCGTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                              <  1:117464/151‑1 (MQ=255)                                                                    tggtgcgccgggggCGCTTCTTCTCTTGCCTACATTTTTTCCATATTTTCCACCCGGCAA‑CCCCCGTTT‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTAc                                                                            <  1:205877/139‑1 (MQ=255)                                                                      atgccgccctcgcgcTTCTTCTCATCCCTATATATTTTTCATATTTTCCATCCGGCAC‑CCCCCTTTC‑CCCCCCCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTGcc                                                                           <  2:212457/139‑1 (MQ=255)                                                                       tGCGCCGCGGGCGCGTCTTATCATGCCTAAATTTTTTTCATATTTTACACCCGGCAC‑CCCCCTTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg                                                                         <  1:138460/151‑1 (MQ=255)                                                                         cgccgcgGGCGCGTTTTATCATGCCTACATTTTTTTCATATTTTACACCCGGCAA‑CCCCCGTTT‑CCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAt                                                                       <  1:566720/148‑1 (MQ=255)                                                                           ccGCTGGCGCGTCTTATCATGCCTACATTTTTTTCATATTTTACATCCGGCAC‑CCCCCGTTT‑CCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                                      <  1:744571/149‑1 (MQ=255)                                                                                     tttttATCACGCCTCCTTTTTTTTCATATTTTCCACCCGCAAC‑CCACCTTTC‑CCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                           <  2:512877/148‑1 (MQ=255)                                                                                     ttCTTATCATGCCTACATATTTTTCATTTTTTACACCCGGCAA‑CCACCTTTT‑CCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                                     <  2:328610/140‑1 (MQ=255)                                                                                     gTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTc                                                                     >  1:328686/1‑141 (MQ=255)                                                                                     gTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                           >  2:229466/1‑151 (MQ=255)                                                                                     gTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGt                                                            >  2:76557/1‑150 (MQ=255)                                                                                         tATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                                       >  2:133014/1‑151 (MQ=255)                                                                                           tCATCCCTAAATTTTTTTCATTTTTACCACCCGGCAC‑CCCCCTTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGa                                                     <  1:954023/151‑1 (MQ=255)                                                                                              tccctcattttttttAATTTTTTCCACCCGGCAC‑CCCCCTTTT‑CCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                  <  1:935397/141‑1 (MQ=255)                                                                                              tccccaCATTTTTTTCATTTTTCCCCCCCGGCCA‑CCCCCTTTT‑CCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                  <  1:433987/146‑1 (MQ=255)                                                                                              tGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                  >  1:704576/1‑151 (MQ=255)                                                                                               gCCTACATATTTTTAATATTTTACATCCTCCAA‑CCACCGTTT‑ACACCGTCACCACCTCACCAGCCGGTGGCGTTTCCAGCAGTTACAGCATGGTACCGGCGATTTCACGCTCGTCCATCACTACCTGATTCGCACCACGATCGGTGAtata                                                 >  2:1162543/1‑151 (MQ=255)                                                                                                ccTACATTTTTTTCATTTTTTCCATCCGGCAA‑CCCCCGTTT‑CCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                <  1:797248/151‑1 (MQ=255)                                                                                                ccTACATTTTTTTCATATTTTACATCCGGCAC‑CCCCCTTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                <  1:541071/151‑1 (MQ=255)                                                                                                ccTACATTTTTTTCATATTTTACACCCGCCCC‑CCCCCGTTC‑CCCCCCTCCCCCCCTCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                <  2:842436/151‑1 (MQ=255)                                                                                                ccTACATATTTTTCATATTTTACATCCGGCAA‑CCACCTTTT‑CCCCCGTCACCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                <  2:482663/151‑1 (MQ=255)                                                                                                ccTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGAtata                                                 >  2:648545/1‑150 (MQ=255)                                                                                                 ctccATATTTTTCATATTTTACACCCGGCAC‑CCCCCGTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                                               <  1:775573/148‑1 (MQ=255)                                                                                                 cTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACg                                               >  1:479246/1‑151 (MQ=255)                                                                                                  cactttttttcCATTTTTTCCCTCGGGAAC‑CCACCTTTT‑ACCCCCCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                              <  1:964200/145‑1 (MQ=255)                                                                                                  cactttttttacttttttcACAGCCGGCAA‑CCCCCTTTT‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc                                              <  2:1007387/136‑1 (MQ=255)                                                                                                   aCATTTTTTTCATTTTTTACACCCGGCAC‑CCCCCGTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGcc                                             <  1:977999/151‑1 (MQ=255)                                                                                                    cATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCa                                            >  2:449704/1‑151 (MQ=255)                                                                                                          tttCATTTTTTCCCCCCGGCAA‑CCCCCTTTT‑CCCCCCCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCAt                                      <  1:532024/151‑1 (MQ=255)                                                                                                            taattttttccccccGGCAA‑CCCCCGTTT‑CCCCCGCCCCCACCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCg                                    <  2:1200857/138‑1 (MQ=255)                                                                                                               tatTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCa                                 >  1:135666/1‑151 (MQ=255)                                                                                                                    tacccccGGCAC‑CCCCCTTTT‑CCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTCCGGTGATATACGCCACTTCATCGTCATAATg                            <  1:76575/146‑1 (MQ=255)                                                                                                                    tACATCCGGCAC‑CCCCCTTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg                            <  1:449850/151‑1 (MQ=255)                                                                                                                     cacaccGGCAACCCCCCGTTT‑CCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg                            <  2:164608/147‑1 (MQ=255)                                                                                                                      cATCCGGCAC‑CCCCCTTTT‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATggg                          <  1:260942/151‑1 (MQ=255)                                                                                                                         cccgCAA‑CCCCCCTTTCACCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg                        <  1:783925/148‑1 (MQ=255)                                                                                                                                    ccGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCCGCATGGTACGGGCGATTTCACGCTCGCCCGTCCCTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaatcat             >  2:98833/1‑151 (MQ=255)                                                                                                                                    ccGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat             >  2:395425/1‑151 (MQ=255)                                                                                                                                    ccGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat             >  1:675901/1‑151 (MQ=255)                                                                                                                                    ccGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaataat             >  2:213745/1‑151 (MQ=255)                                                                                                                                     ctttt‑CCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc            <  2:307364/149‑1 (MQ=255)                                                                                                                                      tttt‑CCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtct           <  1:285956/150‑1 (MQ=255)                                                                                                                                            ccccGTCACCACCTCACCCGCCGGCGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCCTCATAATGGGCGCGGGCAATAATCTCAata      >  1:58340/1‑151 (MQ=255)                                                                                                                                                gTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCg  >  2:970643/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                 AAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAA‑CCACCGTTT‑ACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCG  >  NC_000913/501428‑501719

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