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

CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993763‑994063                                                                                                                                                              |                                                                                                                                                        ccACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCgg                                                                                                                                                         <  2:328332/151‑1 (MQ=255)    cGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGcc                                                                                                                                                       <  2:626608/150‑1 (MQ=255)       tAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccacc                                                                                                                                                   >  2:205876/1‑151 (MQ=255)        aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccacca                                                                                                                                                  <  1:877901/151‑1 (MQ=255)           aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCg                                                                                                                                              <  1:832207/151‑1 (MQ=255)               gcgcATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCcac                                                                                                                                          >  2:972578/1‑151 (MQ=255)                cgcATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCcac                                                                                                                                          >  1:415636/1‑150 (MQ=255)                 gcATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCcaca                                                                                                                                         >  2:179900/1‑150 (MQ=255)                  cATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACaa                                                                                                                                        >  1:142346/1‑150 (MQ=255)                   aTCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAAc                                                                                                                                      >  2:428236/1‑151 (MQ=255)                                 tccgcgTGTC‑TCCTTAA‑TTTCCGGCCGGCGGTGGGGGCCCCGCTACCACATCGCTGGCGTTTGCCTTTC‑CC‑AGC‑CCTG‑CCAGCGGGCGCAGCAGGTTATTTTTG‑CCCCCACCGCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCgg                                                                                                                           <  2:767546/146‑1 (MQ=255)                                  tCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAat                                                                                                                        >  1:259507/1‑150 (MQ=255)                                  tCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGCTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAat                                                                                                                        >  2:1162476/1‑150 (MQ=255)                                    gggtttCTTCCGAAT‑TTT‑CACCCCGGGGTGTGGTGCCCCCAAACACGCCGCCGGCGTTGGGCTTTC‑CA‑AGC‑CCTG‑CCGGCGGGCGCAGCGGGTTCTTTTTC‑CCCCCCCCGCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAata                                                                                                                       <  1:1126204/144‑1 (MQ=255)                                     cgTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                    >  2:533565/1‑151 (MQ=255)                                       tgtttcccTAAA‑TTC‑CAGCCAGGGGTGTGGTGCCCTCAAACACACCGCCTGGGGTGGGCGTTC‑CC‑AGC‑CCTG‑CCAGCGGGCGCAGAAGGTTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCAA‑CCGCCACAAACTGACCTGCCGGAatat                                                                                                                      <  2:979494/141‑1 (MQ=255)                                       tgtCTTCCTGAT‑TTC‑CACCCGGCGGTGTGGTCCCCGCTAACACACCCCCGGCGTTTGGCTTTC‑CC‑AGC‑CCCC‑CCAGCAGGCCCAGCAGGTTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAat                                                                                                                        <  2:39998/145‑1 (MQ=255)                                         tCTTCCTAAA‑TTC‑CAGCCGGGGGGGTGGTCCCCGCAAACACACCGCCTGGGGTTGGCTTTC‑CC‑AGC‑CCTG‑CCAGCGGGCGCAGCAGGGTATTTTTC‑CCCCCCCCGCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGaa                                                                                                                         <  1:449628/142‑1 (MQ=255)                                            tccctaa‑TTTCCAGCCGGCGGTGTGGCCCCCGCTACCAC‑TCCCCCGGCGTTGGCCTTTCCC‑AGC‑CCTC‑CCGGCAGGCGCAGCAGGTTTCTTTTG‑CCCCCCCCGCTGCGCCCCACCA‑CCGCCACAAACTGACCTGc                                                                                                                              <  1:439002/130‑1 (MQ=255)                                             ccTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                               >  1:217457/1‑148 (MQ=255)                                              cctaatttt‑CAGCCGGGGGTGTGGTGCCCCCTAACACCCCCCCGGGGGTTGGCGTTC‑CC‑AGC‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTC‑CCACCCCCGCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCa                                                                                                             <  2:826787/144‑1 (MQ=255)                                              cTGAATTTT‑CAGCCGGGGGGGTGGCGCCCGCAACCACAGCGCCTGGGTTGGGCTTTC‑CC‑AGC‑CCTG‑CCGGAAGCCGCAGGGGGGTACTTTGG‑CCCCCCCCGGGGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                              <  1:758961/149‑1 (MQ=255)                                                    ttt‑CAGCCAGGGGTGGGGTGCCCGCAAACCCACCGCCTGCGGTTGGCTTTC‑CC‑AGC‑CCTG‑CCAGCAGGCGCAGGGGGTTCTTTTTG‑CCCCCACCGCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                       <  2:598626/150‑1 (MQ=255)                                                     tt‑CAGCCGGGGGTGTGGGCCCCGCAAACACACCGCCGGCGGTTGGCTTTC‑CC‑AGC‑CCGG‑CCAGCGGGCGCAGCGGGGTACTTTTC‑CCCCCACCGCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                      <  2:482687/150‑1 (MQ=255)                                                        cagccaggGGTGTGGTGCCCGCTAACACAGCGCCTGCGGTTGGCTTTC‑CC‑AGC‑CCTC‑CCAGCAGGCGCAGCAGGGTACTTTTC‑CCCCCCCCGCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCa                                                                                                    <  2:807101/150‑1 (MQ=255)                                                        cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                     >  1:766971/1‑149 (MQ=255)                                                        cacccgggggTGTGGCCCCCGCTACCACCCCCCCGGGGGTTGGCTTTC‑CC‑AGC‑CCTC‑CCAGAAGGCGCAGCAGGGTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAg                                                                                                   <  1:506023/143‑1 (MQ=255)                                                             gggggTGGGGCGCCGGCAACCCCACCGCCGGCGGTGGGCTTTC‑CC‑ACC‑CCGC‑CCAGCAGGCGCAGCAGGGTATTTTCG‑CCACCCCCGCTGCGCCCCACCG‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAg                                                                                                   <  1:105497/143‑1 (MQ=255)                                                             agcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                               >  2:784643/1‑150 (MQ=255)                                                              ggggTGGGGTCCCCGCAACCACACCCCCGGGGGTGGGCTTTC‑CC‑AGC‑CCTG‑CCAGCGGGCGCAGCGGGTTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                                  <  2:419166/144‑1 (MQ=255)                                                              ggggTGGGGCGCCCCCAAACACCGCCCCGGGGGTGGGCTTTC‑CC‑AGC‑CCCG‑CCAGCGGGCGAAGCGGGTTTCTTTCC‑CCACCCCCGCTGCGGCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                                <  1:974325/146‑1 (MQ=255)                                                               cGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCACCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGCATATGTAAATCCAGTTGGTTCAGGACGat                                                                                             >  1:1073709/1‑150 (MQ=255)                                                                 gtgtGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                          >  2:891858/1‑151 (MQ=255)                                                                  tgtgGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                         >  2:1151525/1‑151 (MQ=255)                                                                       accccagcgccaccTCTCCTGCGGTTGGCTTTC‑CC‑CGT‑CTTG‑CCGCCGGGCCCGGAGGGGTAC‑TTTG‑CCCCCACCCCCGGCTCCCCCCAGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcg                                                                                     <  2:402225/140‑1 (MQ=255)                                                                        gCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg                                                                                   >  2:951435/1‑151 (MQ=255)                                                                         ccctcaaaCACACCGCCTGGGGTTGGCTTTC‑CC‑AGC‑CCGC‑CAAGCAGGCGCAGGAGGTTTCTTTGC‑CCACCACGGCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcg                                                                                     <  1:401900/142‑1 (MQ=255)                                                                         ccccCTAACCCCCCGCCGGGGTTTGGCTTTC‑CC‑AGC‑CCGC‑CCAGCAGGCGCAGCGGGTTTTTTTTC‑CCACCCCCGCTGCGCCCCCCCG‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                         <  1:808225/144‑1 (MQ=255)                                                                           cGCTACCACATCGCCTGGGGTGGGC‑TTT‑CC‑CGACCCTC‑CCAGCGGGCGCAGGAGGTTACTTTTC‑CCCCCCCCGCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                                <  2:687222/151‑1 (MQ=255)                                                                              aaaCACATCCCTGGGGGTTGGCTTTC‑CC‑AGA‑CCTC‑CCAGCGGGCACGGAGGGTTATTTTTC‑CCCCCCCCGCGGCGCCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                              <  2:706980/149‑1 (MQ=255)                                                                               aaCACATCGCCGGGGGTGGGCGTTC‑CC‑AGC‑CCTG‑CCAGCGGGCGCAGCAGGGTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                              <  2:997648/149‑1 (MQ=255)                                                                                 cacaCCGCCGGCGGTGGGCTTTC‑CC‑AGC‑CCTC‑CCGGCGGGCGCAGCGGGTTTTTTTTC‑CCCCCCCCGCTGGGCCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                           <  1:435605/150‑1 (MQ=255)                                                                                  acaTCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                         >  1:838011/1‑151 (MQ=255)                                                                                   caccGCCTGCGTTTGGCGTTC‑CC‑AGC‑CCGC‑CCAGAGGGCGCAGCAGGTTTCTTTTG‑CCACCCCCGCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                          <  2:416212/146‑1 (MQ=255)                                                                                     ccccccggggTTTGCCTTT‑CC‑CGG‑CCCGCCCAGCGGGCGCAGC‑GGGTTCTTTTGCCCCCCCCCGCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                       <  1:127519/143‑1 (MQ=255)                                                                                         cTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                  >  2:564092/1‑151 (MQ=255)                                                                                           gggtttGGCTTTC‑CC‑AGC‑CCGG‑CCAGCGGGCGCGGCAGGTTATTTTGC‑CCCCCCCCGCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                          <  2:302189/137‑1 (MQ=255)                                                                                              ggtGGGGTTT‑CCAAGA‑CCTG‑CCAGCAGGCGCAGCAGGGTTCTTTGG‑CCCCCCCCGCGGGGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                              <  2:1133714/149‑1 (MQ=255)                                                                                              gTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                              >  1:253941/1‑150 (MQ=255)                                                                                               ttGGCTTTC‑CC‑AGC‑CCTG‑CCAGCAGGCGCAGCAGGTTATTTTTC‑CCCCCCCCGCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCa                                                                 <  2:1178484/146‑1 (MQ=255)                                                                                                      tt‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                     >  1:342570/1‑151 (MQ=255)                                                                                                         cc‑AGC‑CCGG‑CCAGCAGGCGCAGCGGGGTACTTTTC‑CCCCCACCGCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg                                                    <  2:157614/150‑1 (MQ=255)                                                                                                              a‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                               >  2:870426/1‑151 (MQ=255)                                                                                                              a‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                               >  2:948594/1‑151 (MQ=255)                                                                                                              a‑CCTG‑CCAGCACCCGCAGCAGGGCACTTCCG‑CCACCACCGCCGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGCCGATATTTTCCGCGTAATGTTTGCTCACTGCATTGAGCAACAATGGCg                                                 >  1:337867/1‑149 (MQ=255)                                                                                                                ccgg‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCCCCCCCGCTGCGCCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGt                                                <  2:916123/146‑1 (MQ=255)                                                                                                                ccTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGc                                              >  2:751452/1‑151 (MQ=255)                                                                                                                   g‑CCAGCAGGCGCAGCAGGGTACTTTTG‑TCACCACCCCTGCGGCCCAACA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCAAGTTGGTTCAGGACGATATTTTCCGCGTAATTTTTGCTTCCTGCATTGTTCAACTATGGCGTGCCCt                                           >  1:614037/1‑151 (MQ=255)                                                                                                                      cagcagGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTgg                                         >  2:685348/1‑151 (MQ=255)                                                                                                                             cgcAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                  >  1:106780/1‑151 (MQ=255)                                                                                                                               cggcagGTTCCTTTCG‑CCCCCCCCGCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                 <  2:444017/148‑1 (MQ=255)                                                                                                                                agcagGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                 >  2:1033425/1‑149 (MQ=255)                                                                                                                                         cTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                       >  1:930829/1‑150 (MQ=255)                                                                                                                                                ccccccccGGTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg                 <  1:1042829/144‑1 (MQ=255)                                                                                                                                                ccccccccGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg                 <  2:177132/144‑1 (MQ=255)                                                                                                                                                ccaccaCCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg                 >  2:964904/1‑150 (MQ=255)                                                                                                                                                 caccacCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg               >  2:1038747/1‑151 (MQ=255)                                                                                                                                                   ccaccGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGc             >  2:449284/1‑151 (MQ=255)                                                                                                                                                    caccGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc            >  1:640879/1‑151 (MQ=255)                                                                                                                                                           gCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCTAGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACACTGGCCTGTCCTGGTTCAGACGCGCAGTATTCATACCGTTGCCTCCtt       >  1:178193/1‑149 (MQ=255)                                                                                                                                                              gCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  >  2:174221/1‑151 (MQ=255)                                                                                                                                                              |                                                                                                                                                        CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAA‑TTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CC‑AGA‑CCTG‑CCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993763‑994063

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