Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,917	0	A	C	42.4%	78.0 / 22.9	59	V43G (GTG→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (28/6);  new base C (0/25);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 4.54e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993770‑994063                                                                                                                                                      |                                                                                                                                                   tgCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCccacc                                                                                                                                                  >  2:316791/2‑150 (MQ=255)  gCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCcaccac                                                                                                                                                >  1:753198/1‑151 (MQ=255)     gACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCATGGTAC‑TTTTG‑CCACCTCCGCTGCGGCCCACCACCg                                                                                                                                              <  2:10099/150‑1 (MQ=255)          gcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCcaca                                                                                                                                         >  2:876078/1‑150 (MQ=255)            aTCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAAc                                                                                                                                      >  1:704454/1‑151 (MQ=255)            aTCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAAc                                                                                                                                      >  2:275712/1‑151 (MQ=255)                   aaCATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTACAGACCTGCCAGCAGGCGCAGAAGGGTAA‑TTTTG‑CAACCACCGCTGCGGCCAACAACCGCCACAAACTGAc                                                                                                                                  >  2:739643/1‑148 (MQ=255)                      atcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                            >  2:179102/1‑151 (MQ=255)                            ccgttttcctCCTGATTTTCAGCCAGGGGGGTGGCCCCCCTTAACACACCCCCGGCGTTTGGCTTTTCAGGACCGCCCAGCGGGCGCAACGGGGTCT‑TTTTC‑CCACCCCCCCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                       <  1:276038/141‑1 (MQ=255)                              cgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                    >  1:78212/1‑151 (MQ=255)                              cgTGTCTCCCTAATTTTCACCCGGCGGGGGGGGGCCCGCTACCCCACCGCCTGGGGTTGGCGTTCCCAGCCCTCCCACCTGCCGCAGCGGGTTAC‑TTTTC‑CCACCTCCGTTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAat                                                                                                                        <  1:817954/147‑1 (MQ=255)                               ttgtCTCCCAGAATTCCAGCCGGGGGTGGGGCGCCCGCTAACCCATCGCGGGCGTTTGGTGTTCCCACCCCGGCCGGCGGGCCCAGCAGGTTTC‑TTGTG‑CCACCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGaa                                                                                                                         <  2:17387/144‑1 (MQ=255)                               gtgtCTTCGTAAATTCCCGCCGGGGGTGGGGTGCCCCCTAACACCCCGCCTGCGTTTGGCGTCCCCGCCCCCCCCAGCGGGCGCGGCGGGTTTT‑TTTTG‑CCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                       <  2:1052486/147‑1 (MQ=255)                                  tCTTCCTGAATTTCAGCCGGGGGTGGGGTGCCCGCAAACACACCGCCGGCGGTTGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCGGGGTCC‑TTTTG‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTaaa                                                                                                                 <  1:943868/150‑1 (MQ=255)                                  tCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                                >  1:200356/1‑151 (MQ=255)                                      cgtaattttCAGCCAGGGGTGTGGTCCCCGCTAACACAGCGCCGGCGGTGGGCGTTTCCAGCCCTGCCGGGGGGCGCAGCAGGGTTC‑TTTTG‑CCACCACCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAg                                                                                                            <  1:416472/145‑1 (MQ=255)                                      ccTGAATTTCAGCCGGGGGTGGGGGCCCCGCTAACACATCGCCGGGGGTGGCCTTTCCCAGCCCTCCCAGCGGGCGCAGAGGGGTCT‑TTTCC‑CCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                     <  2:98866/142‑1 (MQ=255)                                       cTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                           >  1:397598/1‑151 (MQ=255)                                       cTGAATTCCAGCCGGGGGTGGGGGGCCCGCTAACACACCGCCTGCGGTTGGCTTTCCCAGACCTGCCAGCAGGCGCAGCGGGTTCT‑TTTTC‑CCACCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                               <  1:674806/147‑1 (MQ=255)                                        tGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                          >  2:439977/1‑151 (MQ=255)                                         gatttccacccaggGGTGTGGTCCCCGAAAACCCAGCCCTGGGGGTGGGCGTGCCCAGCCCCGCCAGCGGGGGCAGCAGGTTAC‑TTTTG‑CCCCCCCCGGTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                           <  2:367883/140‑1 (MQ=255)                                           aTTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                       >  1:548927/1‑151 (MQ=255)                                            ttccagccgggtgtgtGGTCCCCGCAAACACACCCCCGGCGGTGGGCGTTCCCACCCCTGCCAGCAGGCGCAG‑AGGGTTCTTTTTC‑CCACCACCCCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                       <  2:99989/138‑1 (MQ=255)                                              tCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     >  2:420328/1‑150 (MQ=255)                                              tCAGCCAGCGGTGTGGTGCCCGCTAACACATCGACTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     >  1:1026751/1‑150 (MQ=255)                                               cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     >  2:1035024/1‑149 (MQ=255)                                                       ggtggggcccccGCAAACCCACCGCCTGGGGTTGCTTTTCCCAGCCCGCCCGCCAGGCGCAGCGGGTTCT‑TTTTC‑CCCCCCCCGTTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                           <  2:13163/142‑1 (MQ=255)                                                                 ccGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                  >  1:696566/1‑150 (MQ=255)                                                                      aaCACATCCCCTGTGTTTGGCGTTTCCACCCCTGCCAGCGGCGGCAGCAGGTTTT‑TTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                                <  1:631865/147‑1 (MQ=255)                                                                          caccGCCGGCGGTTGGCGTTCCCAGCCCTGCCAGCAGGCGCAGCGGGGTAT‑TTTTG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCt                                                                        <  1:217291/148‑1 (MQ=255)                                                                           accGCCTGGGGTTGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCAGGGTTC‑TTTTG‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                         <  1:220146/147‑1 (MQ=255)                                                                               ctggCGGTGGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCGGGTTTT‑TTTTG‑CCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTa                                                                      <  1:911838/145‑1 (MQ=255)                                                                               ccTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                   >  2:174631/1‑151 (MQ=255)                                                                                cTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                   >  1:405061/1‑150 (MQ=255)                                                                                   tggTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt                                                               >  2:946998/2‑151 (MQ=255)                                                                                   gggtgggtgTCTCCAGCCCTGACGGCAGGCCAGGCAGGGTAT‑TTTGC‑CCCCCTCCGTTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACt                                                                    <  1:773082/145‑1 (MQ=255)                                                                                      tgggCTTTCCCAGCCCTGCCAGCAGGCGCAGCGGGGTTT‑TTTGG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                             <  1:491351/148‑1 (MQ=255)                                                                                             ttCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                     >  1:121726/1‑151 (MQ=255)                                                                                             ttCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                     >  1:326464/1‑151 (MQ=255)                                                                                               cccgccccgCCGGCAGGCGCAGCAGGTTCC‑TTTTC‑CCCCCCCCGCTCGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                           <  2:1038062/135‑1 (MQ=255)                                                                                               ccacccccgCCAGCGGCCGCAGCAGGTTAT‑TTTTC‑CCCCCGCCGCGGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                   <  2:339101/143‑1 (MQ=255)                                                                                                 agcccccccAGCAGGCCCAGCAGGGTCC‑TTTCC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg                                                    <  1:901463/141‑1 (MQ=255)                                                                                                   cccgcccAGCAGGCGCAGC‑GGGTTT‑TTTTGCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                               <  1:844784/146‑1 (MQ=255)                                                                                                    ccTCCCAGCGGGCGCAGCGGGGTTC‑TTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                 <  1:782526/148‑1 (MQ=255)                                                                                                        ccAGCGGGCGCGGCAGGCTTG‑TTTCC‑CCCCCACCGTCGCCCCCCCCACCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                                          <  1:929402/151‑1 (MQ=255)                                                                                                        ccAGCAGGCGCAGCAGGGTAT‑TTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                           <  1:156715/150‑1 (MQ=255)                                                                                                          agcagGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                        >  1:870522/1‑151 (MQ=255)                                                                                                           gcagGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGtt                                       >  1:219816/1‑151 (MQ=255)                                                                                                                cgcAGCGGGTTTC‑TTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                  <  1:731558/151‑1 (MQ=255)                                                                                                                  cagcagGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                 >  2:237321/1‑150 (MQ=255)                                                                                                                   agcagGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGaag                                 >  2:448643/1‑147 (MQ=255)                                                                                                                     cagGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                              >  2:403651/1‑150 (MQ=255)                                                                                                                       gggTTC‑TTTTC‑CCCCCCCCTCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                           <  1:632419/151‑1 (MQ=255)                                                                                                                        ggttttttttC‑CCCCCCCCGCTGCCCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa                                <  2:871351/140‑1 (MQ=255)                                                                                                                        ggttatttttC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                                   <  1:213293/137‑1 (MQ=255)                                                                                                                        ggTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAATTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                           >  2:451196/1‑150 (MQ=255)                                                                                                                         gTAC‑TTTTT‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGATCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                         >  2:257543/1‑151 (MQ=255)                                                                                                                              ttttC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCa                      <  1:573892/150‑1 (MQ=255)                                                                                                                              ttttC‑CCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                       <  1:340937/149‑1 (MQ=255)                                                                                                                               tttG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                     >  2:830061/1‑150 (MQ=255)                                                                                                                                            gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct        >  2:87251/1‑151 (MQ=255)                                                                                                                                                  cccccccccCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  <  2:252508/150‑1 (MQ=255)                                                                                                                                                  cccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt       <  1:102514/145‑1 (MQ=255)                                                                                                                                                  cccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  <  2:576774/150‑1 (MQ=255)                                                                                                                                                   ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa     <  1:316815/147‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993770‑994063

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