Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	501,575	0	A	C	40.0%	151.5 / 14.8	75	V555G (GTG→GGG)	ybaL	inner membrane putative NAD(P)‑binding transporter
Reads supporting (aligned to +/- strand):  ref base A (43/2);  new base C (0/30);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 6.34e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTT  >  NC_000913/501427‑501725                                                                                                                                                     |                                                                                                                                                       tAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCaccacc                                                                                                                                                      >  2:736650/1‑151 (MQ=255)    gTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTc                                                                                                                                                    >  2:17605/1‑150 (MQ=255)     ttATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCAc                                                                                                                                                  >  1:416627/1‑151 (MQ=255)     ttATCGTCGGTTCCTAGGCCAGATAAAGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCa                                                                                                                                                   >  1:682210/1‑150 (MQ=255)        tcgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccg                                                                                                                                               >  2:689681/1‑151 (MQ=255)         cgtcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccgc                                                                                                                                              >  1:264642/1‑151 (MQ=255)           tcgGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccgccg                                                                                                                                            >  1:145647/1‑151 (MQ=255)               ttCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTgg                                                                                                                                        >  2:1218910/1‑151 (MQ=255)                     ggCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTc                                                                                                                                  >  2:1072008/1‑151 (MQ=255)                       ccAGATAGGGCGTTCCCCCCGCATCTGGCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCTTCTTTTCATCCCTACATTTTTTTCATTTTTTACATCCGGCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTc                                                                                                                                  <  2:639899/149‑1 (MQ=255)                        cAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTcc                                                                                                                                 >  2:127266/1‑149 (MQ=255)                          gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCcagc                                                                                                                              >  2:739814/1‑150 (MQ=255)                             agggcttccccgccgcATCTGGCTTTTGGCTCTCGATGCCTGATGCGCCGCGGGCGCTTCTTATCAGCCCTACATATTTTTCATTTTTTACACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCTCACCCGCCGGTGGCGTTTCcagca                                                                                                                             <  1:963363/138‑1 (MQ=255)                             agggCGTTCCCCCCCCATCTGGCTTTTGCCTCTCGATGCCTGAGGCGACCCTGGCGCTTCTTATAATGCTTAAATATTTTTCATTTTTTCAATCCGCCACCCCCCTTTTCCCCCGTCCCCACCCCCCCCGCCGGTGGCGTTTCCAGCAGt                                                                                                                           <  2:1155496/148‑1 (MQ=255)                               ggCTTTCCCGCCGCATCTGGAATTTGGCTCTCAATGCCTGATGCGACGCTGGCGCGTTTTATCATCCCTACTTTTTTTTCTTATTTTCCACCCGGCAACCACCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCcag                                                                                                                               <  1:480128/144‑1 (MQ=25)                                   tacACGCCCAATCTGGTATTTGGCTCTCGTTGCCTGATGCGCCGCTGGGACTTTTTATCATGCTTACATATTTTTCATATTTTACACCCGGCACCCACCTTTTCCCCCGCCCCCCCCGCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAg                                                                                                                      <  1:1235561/147‑1 (MQ=255)                                         ccgcTTCTGGTTTTTGTCCCTCGTTGCCTGATGCGCCGCTGGCGCTTCTTTTTTTCCCTACTTTTTTTTAATTTTTTACACCCGGCAACCACCGTTTCCCCCGTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCa                                                                                                                    <  1:1086617/145‑1 (MQ=255)                                          cgcTTCTGGCTTTTGCCTCTCGATGCCTGTTGCGCCGCTGGCGCTTCTTATCATGCCTACATTTTTTTAATATTTTACATCCGGAAACCCCCTTTTCCCCCGTCACCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTa                                                                                                               <  2:491922/149‑1 (MQ=255)                                            cATCTGGTTTTGGGCTCTTTATCCCGGAGGCGCCTCGGGCGCGTTTTACCATCCATATATTTTTTACATATTTTACTCCCGGCACCCCCCTTTTCCCCCTTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGt                                                                                                                <  1:708187/146‑1 (MQ=255)                                             atctggcatctggcTCTCGATGCCTGATGCCACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACCTCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTCGCGTTTCCAGCAGCTCCAGCTTGGTACggg                                                                                                           >  1:303356/1‑150 (MQ=39)                                                tGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGTGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAt                                                                                                       >  2:676612/1‑151 (MQ=255)                                                 gtcttttGCCTCTCGAGGCCTGATGCGCCGCTGGCGCGTTTTATCATGCCTACATTTTTTTCATATTTTCCACCCGGCACCCCCCTTTTCCCCCGTCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGAtt                                                                                                      <  1:644726/147‑1 (MQ=255)                                                 ggCTTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTTTTATCATGCCTACATTTTTTTCATATTTTACATCCGGCAACCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCg                                                                                                         <  2:950525/148‑1 (MQ=38)                                                 ggCTTTTGGCTCTCGATCCCTGATGCGCCGCGGGCGCTTTTTATCATGCCTACATATTTTTCATATTTTACATCCGCCACCCCCCGTTTCCCCCGCCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACgg                                                                                                            <  2:94428/145‑1 (MQ=255)                                                        gtctctCGATGCCTGAGGCGACGCTGCCGCGTTTTATCTTCCCAACATTTTTTTCATTTTTTACATCCGCCAACCCCTTTTTCCCCCTCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCt                                                                                               <  1:595367/149‑1 (MQ=255)                                                           tctcGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGc                                                                                            >  1:267422/1‑151 (MQ=255)                                                                    cTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                                    >  1:582199/1‑150 (MQ=255)                                                                       aTGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCt                                                                                >  2:930070/1‑151 (MQ=255)                                                                        tGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg                                                                               >  2:958914/1‑151 (MQ=255)                                                                         gctaCGCTGGCGCGTCTTATCGTCCCTACATTTTTTTCATTTTTTCCACCCGGCACCCACCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACt                                                                                    <  1:722447/142‑1 (MQ=255)                                                                           gccGCTGGCGCGTTTTATCATCCTTACATATTTTTAATATTTTACTCCCGCCACCCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg                                                                               <  2:654233/146‑1 (MQ=255)                                                                           gACGCTGGCGCGTCTTATCATGCCTACATATTTTTTATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                                            >  1:1208980/1‑151 (MQ=255)                                                                             cGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCg                                                                          >  2:277888/1‑151 (MQ=255)                                                                               cTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGc                                                                         >  1:31674/1‑150 (MQ=255)                                                                                 ggCGCGTCTTATCATCCCTACATTTTTTTCATTTTTTACACCCGCCAACCCCCTTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGAtt                                                                            <  1:19341/145‑1 (MQ=255)                                                                                    gcgTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg                                                                   >  2:179988/1‑151 (MQ=255)                                                                                      gTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGtt                                                                 >  2:881190/1‑151 (MQ=255)                                                                                        cTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                                               >  1:713672/1‑151 (MQ=255)                                                                                        cTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                                               >  1:243096/1‑151 (MQ=255)                                                                                          tATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGt                                                             >  2:1238644/1‑151 (MQ=255)                                                                                          tATCATCCCTACATTTTTTTCATTTTTTACATCCGGCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCg                                                               <  2:947219/149‑1 (MQ=255)                                                                                           aTCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTg                                                            >  1:257785/1‑151 (MQ=255)                                                                                               tGCCTACATTTTTTTCATTTTTTCCACCCGAAACCCCCTGTTCCCCCCTCCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGa                                                           <  2:516174/148‑1 (MQ=255)                                                                                               tGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                        >  1:502666/1‑151 (MQ=255)                                                                                               tGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                        >  1:1485/1‑151 (MQ=255)                                                                                                gCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGatat                                                        >  1:1140847/1‑150 (MQ=255)                                                                                                 ccTACATTTTTTTCATTTTTTACATCCGGCAACCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGAtata                                                       <  2:535384/150‑1 (MQ=255)                                                                                                 ccTACATTTTTTTCATTTTTAACAGCGGCCACCCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGAtata                                                       <  1:433759/150‑1 (MQ=255)                                                                                                 ccTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATAc                                                      >  2:70027/1‑151 (MQ=255)                                                                                                 ccTACATATTTTTCATATTTTACACCCGGCAACCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGAtata                                                       <  2:772410/150‑1 (MQ=255)                                                                                                       tatTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACt                                                >  2:487414/1‑151 (MQ=255)                                                                                                         tttttCATATTTTACACCCGGCAACCCCCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCAc                                                 <  1:564731/148‑1 (MQ=255)                                                                                                           tttCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCAt                                            >  2:255110/1‑151 (MQ=255)                                                                                                             tCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCg                                          >  1:779510/1‑151 (MQ=255)                                                                                                                tatTTTACACCCGCCACCCCCCTTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATc                                           <  2:76671/147‑1 (MQ=255)                                                                                                                  ttttACATCCGGCAACCACCGTTTACCCCCTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCAGCGTCAt                                      >  2:538224/1‑150 (MQ=255)                                                                                                                       cgcccGGCACCCCCCTTTTCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATg                                  <  2:904187/146‑1 (MQ=255)                                                                                                                       cATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATggg                                >  2:563376/1‑151 (MQ=255)                                                                                                                          ccGGCACCCCCCTTTTCCCCCGCCCCCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATgg                                 <  2:1077077/147‑1 (MQ=255)                                                                                                                          ccGGCACCCACCGTTTCCCCCGTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc                             <  2:598379/151‑1 (MQ=255)                                                                                                                          ccGGCAACCCCCTTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc                             <  2:427734/151‑1 (MQ=255)                                                                                                                          ccGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcgc                             >  2:458115/1‑151 (MQ=255)                                                                                                                            ggCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGgcg                              <  1:753077/148‑1 (MQ=255)                                                                                                                               aaCCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGc                         >  1:287251/1‑150 (MQ=255)                                                                                                                                 cccccGTTTCCCCCGTCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCgg                           <  2:921747/143‑1 (MQ=255)                                                                                                                                 cccccGTTCCCCCCTTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCa                        <  1:499004/146‑1 (MQ=255)                                                                                                                                 ccaccTTTTCCCCCTCCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGc                         <  1:915205/148‑1 (MQ=255)                                                                                                                                   accGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCaata                     >  1:163249/1‑150 (MQ=255)                                                                                                                                      gTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtct                 >  2:293364/1‑151 (MQ=255)                                                                                                                                       tttCCCCCGCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtct                 <  1:326504/150‑1 (MQ=255)                                                                                                                                             ccGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATc          >  1:1282710/1‑151 (MQ=255)                                                                                                                                             ccGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTCCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATc          >  1:1213730/1‑151 (MQ=255)                                                                                                                                              cGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATcc         >  2:917039/1‑151 (MQ=255)                                                                                                                                              cGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATc          >  2:404686/1‑150 (MQ=255)                                                                                                                                               gTCCCCCCCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGACATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCg        <  1:1218870/151‑1 (MQ=255)                                                                                                                                                tCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCgg       >  2:104680/1‑151 (MQ=255)                                                                                                                                                     accTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGAtttt  >  1:1203417/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                       TAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTT  >  NC_000913/501427‑501725

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