Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,904	0	A	C	47.6%	96.1 / 12.8	63	G47G (GGT→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (24/9);  new base C (0/30);  total (24/39)
Fisher's exact test for biased strand distribution p-value = 2.50e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGG  >  NC_000913/993754‑994051                                                                                                                                                              |                                                                                                                                                    caCCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑Gccacca                                                                                                                                                     <  2:486584/151‑1 (MQ=255)   ccGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCaccacc                                                                                                                                                   <  1:254822/151‑1 (MQ=255)    cGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCg                                                                                                                                                  >  2:455754/1‑151 (MQ=255)         tCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGc                                                                                                                                              <  2:631430/150‑1 (MQ=255)          cctcGGGAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACTCATCGCCTGCGGTT‑GGCGTTGCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCg                                                                                                                                             <  2:796457/147‑1 (MQ=255)          ccACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAGCACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACGACCGCTGCg                                                                                                                                             <  2:163696/150‑1 (MQ=255)           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                  >  1:1220623/1‑150 (MQ=255)                                                                                                                                          cagGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                  >  2:308685/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)                                                                                                                                                          caccacCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg   >  1:68876/1‑151 (MQ=255)                                                                                                                                                          caccacCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg   >  2:176674/1‑151 (MQ=255)                                                                                                                                                           accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg  >  1:37547/1‑151 (MQ=255)                                                                                                                                                           accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg  >  1:152215/1‑151 (MQ=255)                                                                                                                                                              |                                                                                                                                                    CACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCA‑GCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTTCC‑AG‑ACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGG  >  NC_000913/993754‑994051

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