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

GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993768‑994064                                                                                                                                                              |                                                                                                                                                    gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTG‑AATCTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCccac                                                                                                                                                    >  1:647462/1‑151 (MQ=255)     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCaccacc                                                                                                                                                >  2:486743/1‑151 (MQ=255)                     aaCATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTg                                                                                                                                >  2:850082/1‑151 (MQ=255)                         tcatTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                            >  1:603612/1‑151 (MQ=255)                          cttgggcgTGTCTCCGTG‑‑ATTTCCGCCGGGGG‑TGGGGGTCCCCGAAAACACACCCCCTGCGGTGGGCTTTC‑CC‑AGC‑CCGGCCGGCGGGGGCGGCAGGTTCCTTTG‑CCCCCCC‑CCGCTGGGCCCCCCCACCGCCACAAACTGACCTGCCgg                                                                                                                            <  2:230461/145‑1 (MQ=255)                           atTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCTCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGa                                                                                                                           >  2:44609/1‑150 (MQ=255)                             ccgcgTTTCTCCCTA‑AATTTCAGCGGGGGG‑T‑GGGGGGCCCTCAAACACACCCCGGGGGGTTGGCTTTC‑CA‑AGC‑CCTCCAAGAGGGCCCGGGAGGGTTCTTTC‑CCCACCA‑CGGGTCCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAat                                                                                                                         <  2:253240/149‑1 (MQ=255)                              tcgggtttctCCTGAAATTCCAGCCGGGGG‑T‑GTGGCCCCCCCAAACCCCTCGCC‑GGGGGTGGGGTTT‑CCAAGC‑CCGCCCGGCGGGCGCAGAGGGTTCCTTTT‑CCCCCCC‑CGGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAat                                                                                                                         <  2:770884/141‑1 (MQ=255)                              cgcgTGTCTTCCTG‑AATTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                       >  2:460787/1‑151 (MQ=255)                                 gtgtCTTCCTG‑AATTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                     >  2:175432/1‑150 (MQ=255)                                  tttCTTCCTG‑ATTTCCAGCCAGCGG‑T‑GTGGGCCCCCCTAACCCACCGCCTGGGGTGGGCTTTT‑CC‑GGC‑CCTGCCAGCGGGCGCAGCAGGGTACTTTT‑GCCACCC‑CCCCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                        <  1:926150/144‑1 (MQ=255)                                  tgtCTTCCTG‑AATTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                    >  1:617427/1‑150 (MQ=255)                                    tCTTCCTG‑AATTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                                 >  2:636679/1‑151 (MQ=255)                                        ccTG‑AATTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                               >  2:366414/1‑149 (MQ=255)                                        ccTG‑AATTCCAGCCAGGGG‑T‑GGGGTCCCCCCTAACCCCTCGCCGGCGGTTGGTGTTC‑CC‑AGC‑CCTGCCAGCAGGCGCGGCAGGGTTCTTTT‑GCCCCCC‑CCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAg                                                                                                             <  1:757337/151‑1 (MQ=255)                                        ccTA‑AATTTCGGCCGGCGG‑T‑GTGGTCCCCGCAAACACACCCCCGGGGTTTGGCGTTC‑CC‑AGC‑CCTGCCAGCAGGCGCAGCAGGGTATTTTG‑CCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTaaa                                                                                                                  <  1:66499/146‑1 (MQ=255)                                           g‑AATTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGC‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCTCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                            >  2:136489/1‑149 (MQ=255)                                            taaTTTCCACCCGCGGTT‑GTGGTCCCCCCTAACCCCCCGCCTGCGGTGGGCTTTC‑CC‑CGC‑CCTCCCCGCGGGCGCAGCGGGTTACTTTT‑CCCACCA‑CCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCa                                                                                                              <  2:374588/147‑1 (MQ=255)                                              aTTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                        >  1:452336/1‑151 (MQ=255)                                              aTTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                        >  1:950133/1‑151 (MQ=255)                                                 tCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGCACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATACGTAAATCTAGTTGGCTc                                                                                                      >  1:351247/1‑150 (MQ=255)                                                 ccagccgggGG‑T‑GTGGTGCCCTCTAACACACCGCCTGCGGTTGGCTTTC‑CC‑AGC‑CCTCCCAGCGGGCGCAGCGGGGTTTTTTT‑GCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                      <  1:717804/149‑1 (MQ=255)                                                  cagccggcGT‑T‑GTGGTCCCCGCAAACACACCGCCTGCGGTTGGCTTTT‑CC‑AGC‑CCTCCCAGCAGGCGCAGCGGGGTACTTTT‑CCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                      <  1:638864/149‑1 (MQ=255)                                                     ccggggg‑T‑GGGGTCCCGGCTACCCCACCCCCGGCGGTGGGCTTTC‑CC‑GCC‑CCTGCCAGCGGGCTCAGCAGGTTCCTTTC‑CCCCCCA‑CCGGTTCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCa                                                                                                     <  1:915239/142‑1 (MQ=255)                                                     ccggggg‑T‑GGGGTACCCGCAACCACACCCCGTGGGGTCGGC‑TTC‑CC‑AAACCCCGCCGGAAGGCGCAGGGGGGTATTTTT‑GCCCCCA‑CCGGTGCGGCCCACCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAg                                                                                                    <  2:538192/143‑1 (MQ=255)                                                        gggg‑T‑GTGGTGCCCGCTAACACACCGCCGGCGTTGGGCGTTC‑CC‑AGC‑CCTGCCAGCAGGCGCAGCGGGGTATTTTT‑CCCCCCC‑CCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                              <  1:258684/149‑1 (MQ=255)                                                           g‑t‑gtGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                           >  1:314697/1‑151 (MQ=255)                                                            gggggggtcCCCGCAACCACACCGCCTGGTTTTGGCTTTC‑CC‑ACC‑CCCGCCAGCAGCCGCAGCGGGTTTTTTTT‑GCCACCA‑CCGCTGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                       <  2:161381/135‑1 (MQ=255)                                                                 ggtgCCCGCTAACCCATCCCCTGGGGTTGGCTTTC‑CC‑AGC‑CCTGCCAGCAGGCGCAGCAGGTTACTTTT‑CCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                         <  1:765834/149‑1 (MQ=255)                                                                 ggtcCCCGCTAACACACCGCCTGGGGTGGGCTTTC‑CC‑AGC‑CCTCCCAGCGGGCGCAGCGGGTTATTTTG‑CCCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                       <  1:123037/151‑1 (MQ=255)                                                                          tAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                              >  2:726847/1‑151 (MQ=255)                                                                           accacaCCGCCTGCGGTTGGCTTTC‑CC‑GGC‑CCTCCCAGCAGGCGCGGCGGGTTATTTTT‑CCCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTa                                                                                  <  1:893450/144‑1 (MQ=255)                                                                            accacaccgccggggTTGGCCTTTCCC‑AGC‑CCTCCCAGCAGGCGCAGCGGGTTCTTTTT‑CCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                               <  2:743608/136‑1 (MQ=255)                                                                             cacaTCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                           >  2:888318/1‑151 (MQ=255)                                                                              acaTCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                             >  2:505769/1‑148 (MQ=255)                                                                                aTCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                        >  2:737756/1‑151 (MQ=255)                                                                                   ggcTGTGGTTGG‑GTTT‑CC‑CAC‑CCCGCCCGCGGGCGGCGCGGGGTTCTTTTCGCCCCCA‑CGGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                            <  2:882402/142‑1 (MQ=255)                                                                                     cTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                   >  1:259817/1‑151 (MQ=255)                                                                                         ggTTGGCTTTC‑CC‑AGC‑CCTGCCAGCGGGCGCAGCAGGTTATTTTT‑CCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                      <  1:700903/144‑1 (MQ=255)                                                                                         ggTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                               >  1:565902/1‑151 (MQ=255)                                                                                         ggTTGGCGTTC‑CC‑AGC‑CCTGCCGGCAGGCGCAGCAGGGTTTTTTT‑CCCCCCA‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCa                                                                  <  2:55098/148‑1 (MQ=255)                                                                                             ggCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                            >  1:880709/1‑150 (MQ=255)                                                                                                gTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                         >  2:954777/1‑150 (MQ=255)                                                                                                 gctttcccaccccTGCCAGCGGGCGCAGCAGGTTTTTTTT‑CCCCCCC‑CCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                                          <  2:730886/140‑1 (MQ=255)                                                                                                  tt‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                      >  2:896813/1‑151 (MQ=255)                                                                                                           cccTGCCGGCGGGCGCAGCGGGGTCTTTTT‑CCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                                <  1:777115/150‑1 (MQ=255)                                                                                                            ccggCCACCGGGCGCAGCGGGTTATTTTC‑GCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                  <  2:767428/145‑1 (MQ=255)                                                                                                                      ggCGCGGCAGGTTCCTTTT‑CCCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTgg                                          <  2:950001/146‑1 (MQ=255)                                                                                                                      ggCGCAGCGGGTTACTTTC‑CCCCCCC‑CCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                            <  1:574627/144‑1 (MQ=255)                                                                                                                      ggCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                       >  2:268340/1‑149 (MQ=255)                                                                                                                       gcgcAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                       >  2:8271/1‑148 (MQ=255)                                                                                                                             cagGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAg                                >  2:240122/1‑149 (MQ=255)                                                                                                                               ggttcttttt‑CCCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                            <  1:713602/145‑1 (MQ=255)                                                                                                                               ggggtttttt‑GCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                              <  1:207908/144‑1 (MQ=255)                                                                                                                                 gTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                          >  2:131870/1‑151 (MQ=255)                                                                                                                                           cccccctccGCTGCCGCCCCCCGCCGCCACAACCTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc                   <  2:234957/143‑1 (MQ=255)                                                                                                                                                 cccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct            <  1:788167/150‑1 (MQ=255)                                                                                                                                                    gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct         >  2:563272/1‑151 (MQ=255)                                                                                                                                                         ggCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaa     >  2:852816/1‑150 (MQ=255)                                                                                                                                                            ccaccaCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATg  >  1:295635/1‑150 (MQ=255)                                                                                                                                                              |                                                                                                                                                    GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGG‑T‑GTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993768‑994064

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