Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,913	0	G	C	45.2%	78.0 / 27.8	62	G44G (GGC→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (31/3);  new base C (0/28);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 2.57e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCT  >  NC_000913/993764‑994057                                                                                                                                                        |                                                                                                                                                 cACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGc                                                                                                                                                 >  1:25568/1‑151 (MQ=255)       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)                                                                                                                                                        |                                                                                                                                                 CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCT  >  NC_000913/993764‑994057

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