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

ACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAA  >  NC_000913/993765‑994061                                                                                                                                                          |                                                                                                                                                     aCGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGcc                                                                                                                                                    >  2:388402/1‑151 (MQ=255)  cGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGcc                                                                                                                                                    <  1:1269877/150‑1 (MQ=255)  cGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCATTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGccc                                                                                                                                                   >  2:450452/1‑151 (MQ=255)   ggTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACTTCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCcca                                                                                                                                                  <  1:924767/151‑1 (MQ=255)   ggTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGccc                                                                                                                                                   >  2:93704/1‑150 (MQ=255)      aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccacc                                                                                                                                                >  2:962814/1‑150 (MQ=255)         aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCg                                                                                                                                            >  2:135263/1‑151 (MQ=255)           aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGcc                                                                                                                                          >  2:1230926/1‑151 (MQ=255)              ctcATCCTGAAACTTCATTCGCGTGTCTTCCTGAATTA‑CAG‑CCAGGGTTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTATTTTTG‑CCACCACCGCTGCGGCCCACCACCGcc                                                                                                                                          <  1:823705/146‑1 (MQ=255)                 aTCCTAAAACTTCTTCCGGGTTTCTTCCTGATTTC‑CAC‑CCGGGGTTGTGGTCCCCGCTAACACACCGCC‑GGCGTTTGGCGTTC‑CCACCCCCGCCAGCAGGCGCAGAAGGGTCTTTTTC‑CCCCCCCCGCTGCGGCCCACCACCGCCACaa                                                                                                                                      <  1:427267/149‑1 (MQ=255)                     tGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGAc                                                                                                                                >  1:468188/1‑151 (MQ=255)                              atTCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGa                                                                                                                        >  1:914191/1‑150 (MQ=255)                              atTCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGa                                                                                                                        >  1:866608/1‑150 (MQ=255)                                tCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAat                                                                                                                      >  1:165249/1‑150 (MQ=255)                                   cgTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                  >  1:1468631/1‑151 (MQ=255)                                         ttCCTGATTTT‑AAG‑CCAGGGGGGTGGTGCCCCATAACACACCGCC‑TGCGGTGGGCTTTT‑CCAGCCCGCCCAGCGGGCGCAGAGGGGTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                 <  2:422122/146‑1 (MQ=255)                                         ttCCTGATTTC‑CAG‑CCGGGGGTGGGGCCCCCGCAAACACACCGCG‑GGGGTTTGGTTTTC‑CCAGCCCGGCCAGCAGGCGCAGCGGGGTTTTTTGG‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGaa                                                                                                                       <  2:1462988/140‑1 (MQ=255)                                         ttCCTAATTTC‑CAC‑CCGGGGGTGGGGGCCCCCCTAACCCACCCCC‑GGGGTTTGGCTTTC‑CCAGCCCTCCCGGCGGGCGCAGCGGGGTCTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAat                                                                                                                      <  2:1448865/141‑1 (MQ=255)                                          tCCTTAATTT‑CCGCCCGGGGGTGGGGGCCCCCCAAACCCACCGCC‑TGGGTTGGGCTTTC‑CCAGCCCGGCCGGCAGGCGCAGCGGGTTTTTTTCC‑CCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGaa                                                                                                                       <  2:1341362/140‑1 (MQ=255)                                          tCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                            >  2:457992/1‑150 (MQ=255)                                           ccTGAATTT‑CGG‑CCGGGGGGGGGGTGCCCTCTACCACATCCCC‑TGGGTTGGGCTTTC‑CCAGCCCGCCCAGCGGGCGCAG‑AGGTTTCTTTTGCCCCCCCCCGCTGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                              <  2:523148/147‑1 (MQ=255)                                           ccTGAATTT‑CAG‑CCAGGGGTGCGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGATGCGGCCCACCACCGCCACAAACTGACCCGCCGGAATATGTAAATcc                                                                                                            >  1:495622/1‑149 (MQ=255)                                           ccTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                            >  1:384198/1‑149 (MQ=255)                                           ccTAAATTT‑CAG‑CCGGGGGTGGGGTGCCCGCTACCCCACCGCC‑GGCGGTGGGCTTTC‑CCAGCCCCGCCGGCGGGCGCAGCGGGTTCCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCa                                                                                                           <  2:1207695/150‑1 (MQ=255)                                            cTGAATTT‑CCG‑CCAGGGGTTGGGGGCCCCCTTAACCCCCCCC‑GGGGGTTGGCCTTTCCCCGCCCGCCCAGAGGGCGCAG‑GGGGTTCTTTGGCCCACCACCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                 <  1:653447/144‑1 (MQ=255)                                             tGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                        >  1:1558531/1‑151 (MQ=255)                                              taaTTTCCAG‑CCAGGGGTGTGGCCCCCGC‑AACCAATCCCCTTGGGGTGGGTTTTC‑CCACCCCTCCAAGCAGGCGCAGAGGGTTATTTTCC‑CCCCCCCCGCGGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAatat                                                                                                                    <  2:1544682/138‑1 (MQ=255)                                                  ttccag‑ccgggGGTGGGGTGCCCGCAAACACACCGCG‑TGGGGTTGGCTTTC‑CCACCCCTCCCAGCGGGCGCAGCAGGTTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                      <  2:1047577/146‑1 (MQ=255)                                                   ttccggccgggggTGTGGCCCCCGCAAACACACCGCG‑GGGGGTGGCCTTTC‑CCAGCCCGGCGAGGGGGCGCAGGGGGGTTCTTTCG‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                       <  1:718803/137‑1 (MQ=255)                                                     cag‑ccgggGGTGTGGTGCCCGCAAACACACCGCT‑GGGGTTTGGCTTTC‑CCAGCCCTGCCAGCGGGCGCAGCAGGTTCTTTTTC‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                   <  2:765139/149‑1 (MQ=255)                                                     cag‑ccagcGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                   >  2:912557/1‑149 (MQ=255)                                                       g‑ccagcGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                                >  1:1182756/1‑150 (MQ=255)                                                         ccagcGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                              >  1:1112783/1‑151 (MQ=255)                                                            ggggtgggggccccGCAAACCCACCGCT‑GGCGGTGGGCTTTC‑CCAGCCCTCCCAGCAGGCGCAGCGGGGTACTTTCC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                   <  2:87970/132‑1 (MQ=255)                                                               ggggggcccccgtaaaCACCTCGCT‑TGCGTTTGTCTTTC‑CCACCCCTGCCAGCGGCCGCAGCAGGTTTTTTTGC‑CCCCCCCCCCGGCGCCCCACCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                   <  2:384228/126‑1 (MQ=255)                                                                 gtggtgCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                      >  2:1364632/1‑151 (MQ=255)                                                                   ggtcCCCGCTACCCCCCCGCC‑TGGGGTTGGCTTTC‑CCAGCCCCCCCGGCGGGCGCAGCAGGTTTTTTTCC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                    <  2:843729/151‑1 (MQ=255)                                                                   gggccccGCAAACACACCGCG‑GGGGTTGGGCGTTC‑CCAGCCCGGCCAGCGGGCCCGGGAGGGTACTTTCC‑CCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                    <  2:156098/147‑1 (MQ=255)                                                                       cccGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                >  2:1175480/1‑151 (MQ=255)                                                                        ccGCAAACACTCCGCG‑GGGGGTTGGTTTTC‑CCAGCCCGCCCAGCGGGCGCAGGAGGTTTTTTTTC‑CCCCCCCGGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTa                                                                               <  1:421519/151‑1 (MQ=255)                                                                                 ccccccc‑TGTGTTTGGCTTTC‑CCAGCCCGCCCAGCCGGCGCAGCGGGTTTCTTTTC‑CCCCCCCCGCGGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                         <  1:225209/143‑1 (MQ=255)                                                                                 caccGCC‑GGGGGTTGGTTTTC‑CCAGCCCCGCCAGCAGGCGCAGCGGGTTATTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                           <  2:596979/143‑1 (MQ=255)                                                                                      cc‑TGCGGTTTGCGTTT‑CCAGCCCTTCCAGCAGTCGCAGCAGGGTACTTTGG‑CCACCACCGCTTCGGCCCACCACCGCCACAAAATGACCTGCTGGAATATGTCAATGCAGTTGGCTCAGGACGATATTTTCCGCGTAATGTTTGCTTACt                                                                  >  1:719073/1‑150 (MQ=255)                                                                                         tGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                >  2:154722/1‑150 (MQ=255)                                                                                          gCGGTTGGCTTTC‑CCACCCCGCCCAGCGGCCGCAGCGGGTTATTTTTC‑CCACCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTa                                                                    <  1:1125989/145‑1 (MQ=255)                                                                                          gCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGACTATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                              >  1:1332638/1‑151 (MQ=255)                                                                                          gCGGTGGGCTTTC‑CCAGCCCGGCCAGCAGGCGCAGCAGGTTTTTTTTG‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                     <  1:1314414/144‑1 (MQ=255)                                                                                            gttgggCTTTC‑CCAGCCCGCCCAGCGGGCGCAGAAGGTTATTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                     <  1:868264/138‑1 (MQ=255)                                                                                               tgggctttcccaccccTGCCACCGGGCGCAGAGGGTTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                 <  1:58686/130‑1 (MQ=255)                                                                                                ggCTTTC‑CCAGCCCGCCCAGCGGGCGCAGCGGGTTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                     <  2:552104/138‑1 (MQ=255)                                                                                                   gTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                      >  1:357287/1‑150 (MQ=255)                                                                                                          aGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                               >  1:587809/1‑151 (MQ=255)                                                                                                              cTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGcc                                           >  2:501502/1‑151 (MQ=255)                                                                                                                    gcagGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGtt                                     >  2:715116/1‑151 (MQ=255)                                                                                                                           cagcggTGTTCTTTTC‑CCCCCCCCGCGGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                <  1:456784/149‑1 (MQ=255)                                                                                                                           cagcagGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                >  1:1167679/1‑149 (MQ=255)                                                                                                                             gcagGGTACTTTTG‑CCACCACCGCTGCGGCCCACCCCCGCCACAACCTGACCTGCGGGAATTTGTAATTCCAGTTGGTTCGGGACGAAATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCACCAATGGCGTGCCCTGGTTAAGACggg                             >  1:1222045/1‑148 (MQ=255)                                                                                                                                      tttcG‑CCCCCCCCGCGGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                     <  1:16784/149‑1 (MQ=255)                                                                                                                                       ttcg‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                   <  1:846102/147‑1 (MQ=255)                                                                                                                                         tG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc                >  2:1039215/1‑151 (MQ=255)                                                                                                                                          g‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg               >  2:357839/1‑151 (MQ=255)                                                                                                                                             cccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg             <  1:1468147/146‑1 (MQ=255)                                                                                                                                              accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg            >  2:515573/1‑151 (MQ=255)                                                                                                                                                      tGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt     >  2:927230/1‑150 (MQ=255)                                                                                                                                                         ggCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaa  >  1:1300344/1‑150 (MQ=255)                                                                                                                                                         ggCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaa  >  1:1568235/1‑150 (MQ=255)                                                                                                                                                          |                                                                                                                                                     ACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCC‑TGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAA  >  NC_000913/993765‑994061

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