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

CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA  >  NC_000913/993764‑994060                                                                                                                                                             |                                                                                                                                                    cACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑CCCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGc                                                                                                                                                    <  2:813994/151‑1 (MQ=255)        gCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GACGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCcaccac                                                                                                                                             <  1:429508/151‑1 (MQ=255)         cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCaccacc                                                                                                                                            >  2:637550/1‑151 (MQ=255)             cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCca                                                                                                                                        <  1:1142149/151‑1 (MQ=255)                  aTCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAAc                                                                                                                                   >  2:905525/1‑151 (MQ=255)                      tGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGAc                                                                                                                               >  1:62863/1‑151 (MQ=255)                        aaaaATCTTTGGCGTTTTTTCTTTAATTC‑CACCCG‑GGGGGGGGGTCCCCCCAAACAGCCCGCTTGGGGTG‑GGCTTTCCC‑AC‑CCCGCCCAGCGGCCGCAGAAGGGTTC‑TTTC‑GCCCCCACCGCGGCCCCCCACCACCGCCACAAACTGAc                                                                                                                               <  2:298473/149‑1 (MQ=255)                        aaaCATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCt                                                                                                                             >  1:555274/1‑151 (MQ=255)                         aaCATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTg                                                                                                                            >  1:413895/1‑151 (MQ=255)                           catcatTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGcc                                                                                                                          >  2:1083206/1‑151 (MQ=255)                              catTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                        >  1:417386/1‑150 (MQ=255)                                  cgcgTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                   >  2:1099747/1‑151 (MQ=255)                                   gcgTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                  >  2:790467/1‑151 (MQ=255)                                   gcgTGTCTCCCTAATTTC‑CGGCCG‑GGGGTGTGGCCCCCGCTAACACCCCGCCGGGGGTG‑GGCTTTCCC‑AG‑CCCTGCCAGCAGGCGCAGCGGGTTCT‑TTTT‑CCCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGaa                                                                                                                      <  2:31305/147‑1 (MQ=255)                                    gctttCTTTCTGAATTTCCAGCCG‑GGGGTGGGGCCCCCCCAAACCCACCGCCGGGGGTG‑GGCTTTCCC‑GG‑CCCGGCGGGCAGGCGCAGCGGGTTTT‑TTTT‑CCCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGaa                                                                                                                      <  2:1201537/143‑1 (MQ=255)                                    cgTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGCAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                 >  1:321105/1‑151 (MQ=255)                                      tgtCTTCCTGATTTC‑CACCCG‑GGGGTGGGGTCCCCGCAAACCCATCGCCGGCGTTT‑GGCTTTCCC‑AG‑CCCGCCCGGCGGGCGCAGCGGGGTCT‑TTTC‑CCCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                    <  2:640757/146‑1 (MQ=255)                                        tCTTCCTAATTTT‑CAGCCG‑GGGGGGTGGTCCCCCCAAACCCACCGCCGGCGGTT‑GGCTTTCCC‑AG‑CCCTGCCAGCAGGCGCAGCGGGGTAC‑TTTG‑CCCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                        <  1:1256929/140‑1 (MQ=255)                                        tCTTCCTAATTTC‑CACCCG‑GGGGTGTGGTCCCCGCAAACACCCCCCCTGGGTTT‑GGCTTTCCC‑AG‑CCCCCCCGGCAGGCGCAGCGGGTTAT‑TTTT‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                               <  1:507165/149‑1 (MQ=255)                                        tCTTCCTAAATTT‑CAGCCG‑GGGGTGGGGCGCCCGCAACCACATCGCTGGGGGTT‑GGCTTTCCC‑AG‑CCCTGCCAGCGGGCGCAGCGGGGTTC‑TTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                   <  1:1244218/145‑1 (MQ=255)                                          ttCCTAATTTC‑CAGCCG‑GGGGTGGGGTCCCCGCAAACACATCG‑CTGGGGTTGGGCTTTCCC‑AG‑CCCTGCCGGCGGGCGCAGGGGGGTTT‑TTTT‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                   <  1:1246879/143‑1 (MQ=255)                                            cctaatttt‑CAGCCG‑GCGGTGGGGTCCCCGCAACCCCACCGCCTGGGGTG‑GGCTTTCCC‑AG‑CCCGCCCAGCGGGCGCAGCGGGTTAT‑TTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                             <  1:497100/141‑1 (MQ=255)                                            ccTGATTTT‑CACCAG‑GGGGGGGGGGGCCCCCAAACACACCGCCGGCGGTG‑GGCTTTCCC‑AG‑CCCTGCCAGCAGGCGCAGAAGGGTTT‑TTTT‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                             <  1:552634/147‑1 (MQ=255)                                            ccTGAATTT‑CGCCCG‑GGGGTGTGGCCCCCGCAAACACATCGCCGGGGGTG‑GGTGTTCCC‑AC‑CCCTGCCAGCAGGCGCAGCGGGGTAC‑TTTC‑CCCACCACCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                           <  2:814263/149‑1 (MQ=255)                                             cgtgatttccgccca‑ggGGTGGGGTCCCCCCAAGCCCACCGCCGGGGGTT‑GGCGTCCCC‑CG‑CCCGGCAAGCGGGGGCCGCAGGTTCT‑TTTC‑GCCACCCCCGCTGCGCCCCCCAACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                           <  2:296742/137‑1 (MQ=255)                                               taaTTT‑CCGCCCGGCGGTGTGGTGCCCGCAAACACACCGCCGGGGGTT‑GGCTTTCCC‑AG‑CCCGGCCGGCGGGCGCAGCAGGGTAC‑TTTT‑GCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCa                                                                                                          <  1:956291/147‑1 (MQ=255)                                                atttt‑CAGCCG‑GGGGTGTGGCGCCCGCTACCACACCGCCGGGGGTG‑GGC‑TTTCC‑AGCCCCGGCCAGCGGGCGCAGGAGGGTCT‑TTTC‑CCCCCCACCTCGCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                           <  2:785106/143‑1 (MQ=255)                                                          ca‑gcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                            >  1:812567/1‑151 (MQ=255)                                                               ggTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                        >  2:652567/1‑151 (MQ=255)                                                                 tggggtcCCCGCTAACCCCCCGCCGGGGGTG‑GGCTTTCCC‑AC‑CCCGCCCAGCGGGCGCAGCGGGGTCT‑TTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGa                                                                                           <  1:928044/143‑1 (MQ=255)                                                                  gtggcccccccaaaCACACCGCTTGGGTTT‑GGCTTTCCC‑AG‑CCCTGCCAGGGGGCGCAG‑AGGGTTCTTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                      <  2:555269/138‑1 (MQ=255)                                                                    ggtgCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                   >  1:86690/1‑151 (MQ=255)                                                                    ggtcCCCGCAAACACACCCCCGGGGGTG‑GGCTTTCCC‑AG‑CCCTGCCAGCGGGCGCAGCGGGGTAC‑TTTT‑GCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                       <  1:241517/147‑1 (MQ=255)                                                                     gtgCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcg                                                                                  >  2:837903/1‑151 (MQ=255)                                                                        ccccctacccctccGCCGGCGTTT‑GGCTTTCCC‑AG‑CCCTGCCAGCGGCCGCAGCGGGTTTT‑TTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                               <  1:277276/138‑1 (MQ=255)                                                                        cccGCAACCACACCG‑CTGGGGTT‑GGGCTTTCCAAG‑CCCTGCCAGCAGGCGAAGCAGGGTAT‑TTTC‑CCCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc                                                                                 <  2:1169075/149‑1 (MQ=255)                                                                              accaccgcgcgTGCGTTG‑GGCGTCCCA‑AC‑ACCCGCCAGCGGGCGCAGAGGGTTTT‑TTTT‑CCCCCCCCCGCGGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                           <  2:896152/142‑1 (MQ=255)                                                                              aaCACACCGCCTGGGTTG‑GGCTTTCCC‑AG‑CCCTGCCAGCGGGCGCAGCGGGGTTT‑TTTT‑CCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                           <  1:828786/149‑1 (MQ=255)                                                                                 acaTCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                      >  1:866660/1‑151 (MQ=255)                                                                                  caTCGCGGGGGGTG‑GGCTTTCCC‑AG‑CCCTCCCGGCGGGCGCAGGGGGGTTT‑TTTT‑CCCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                       <  1:61755/149‑1 (MQ=255)                                                                                  caTCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCt                                                                     >  2:334152/1‑151 (MQ=255)                                                                                     cGCCTGCGGTC‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                  >  1:897215/1‑151 (MQ=255)                                                                                       ccggggttt‑GGCTTTCCC‑AC‑CCCGGCCGGCAGGCGAAGC‑GGGTTC‑TTTTCCCCCCCCCCGCGGCGCCCCACAACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                    <  1:553624/140‑1 (MQ=255)                                                                                       ccTGGGGTG‑GGCTTTCCC‑AG‑CCCGGCCAGCGGGCGCAGCAGGGTTC‑TTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                    <  1:765000/147‑1 (MQ=255)                                                                                          gCGTTT‑GGCTTTCCC‑AG‑CCCTGCCAGCGGGCGCAGCGGGGTAT‑TTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                       <  2:353131/141‑1 (MQ=255)                                                                                           gggttgggCGTTTCC‑AG‑CCCGCCGAGCAGGCGCAGCGGGTTTT‑TTTT‑CCCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                  <  1:775023/140‑1 (MQ=255)                                                                                              tt‑GGCGTTCCC‑AG‑CCCGGCCAGCAGGCGCAGAAGGTTAT‑TTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                               <  1:745403/145‑1 (MQ=255)                                                                                                  gCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                       >  1:287107/1‑150 (MQ=255)                                                                                                     tttCC‑AG‑CCCTGCCAGCAGGCGCAGCAGGGTTT‑TTTT‑CCCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaaca                                                    <  1:418902/150‑1 (MQ=255)                                                                                                     ttccc‑AG‑CCCTGCCAGCAGGCGCAGCGGGTTTT‑TTTC‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                                      <  1:1265661/145‑1 (MQ=255)                                                                                                       ttcccaccccTGCCAGCGGGCGCAGCGGGTTTT‑TTTC‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaaca                                                    <  2:771749/142‑1 (MQ=255)                                                                                                              cccgcccAGCGGGCGCAGCGGGTTAC‑TTTC‑GCCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                               <  2:379866/143‑1 (MQ=255)                                                                                                                  cccggcgggCGCAGCGGGTTCC‑TTTT‑CCCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                        <  2:491947/144‑1 (MQ=255)                                                                                                                       cgggCGCAGC‑GGGTTC‑TTTTCCCCCCCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                        <  1:1011544/144‑1 (MQ=255)                                                                                                                                cagGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                           >  2:308685/1‑150 (MQ=255)                                                                                                                                cagGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                           >  1:1220623/1‑150 (MQ=255)                                                                                                                                 agGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                           >  1:15964/1‑149 (MQ=255)                                                                                                                                    ggctcttttccccacccCCGCGGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                          <  1:1070358/138‑1 (MQ=255)                                                                                                                                         tttt‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                  >  1:1131684/1‑151 (MQ=255)                                                                                                                                           tt‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc                >  2:800566/1‑151 (MQ=255)                                                                                                                                               ccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc               <  2:256244/143‑1 (MQ=255)                                                                                                                                                caccacCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg            >  1:68876/1‑151 (MQ=255)                                                                                                                                                caccacCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg            >  2:176674/1‑151 (MQ=255)                                                                                                                                                 accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg           >  1:152215/1‑151 (MQ=255)                                                                                                                                                 accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg           >  1:37547/1‑151 (MQ=255)                                                                                                                                                  cccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGc          <  1:1028672/148‑1 (MQ=255)                                                                                                                                                        cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt    >  1:966351/1‑151 (MQ=255)                                                                                                                                                         tGCGGCCCAACACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGCTTCAGGACGATCTTTTCCGCGTAATGGTTACTTACGGCATTGAGCAACAATGGCGTGCTCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc   >  2:536223/1‑151 (MQ=255)                                                                                                                                                           cGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa  >  1:359189/1‑150 (MQ=255)                                                                                                                                                             |                                                                                                                                                    CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA  >  NC_000913/993764‑994060

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