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

CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993764‑994058                                                                                                                                                            |                                                                                                                                                  cACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGc                                                                                                                                                  >  1:684565/1‑151 (MQ=255)   cGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGccc                                                                                                                                                <  2:310055/151‑1 (MQ=255)     gTAGCAGACGCGCATCCTGGAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCccac                                                                                                                                              >  1:789752/1‑151 (MQ=255)      tAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCccacc                                                                                                                                             >  2:1072015/1‑151 (MQ=255)       aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCccacca                                                                                                                                            <  1:162730/151‑1 (MQ=255)                gcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCcaca                                                                                                                                    >  2:770971/1‑150 (MQ=255)                 cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACaa                                                                                                                                   >  2:139623/1‑150 (MQ=255)                   tCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACt                                                                                                                                >  1:237866/1‑151 (MQ=255)                           catcttCCGCGTTTTTTCTTGATTTCCAC‑CGGGCGGTGTGGTCCCCGCTAACCCCTCC‑CCGGGGTTGGGCTTTC‑CCACACCTGCC‑AGCGGGCGCAGCAGG‑TTCTTTTCG‑CCCCCCCCGCTGGGCCCCACCACCGCCACAAACTGACCTg                                                                                                                          <  2:26134/149‑1 (MQ=255)                                      tttCTTCCTAATTTTCAG‑CCAGGGGTGTGGTCCCCGTAAACTCCGCC‑CCTGGTTTGGGCTTTC‑CCAGCCCCGCC‑AGCAGGCGCAGCAGG‑TTTTTTTGC‑CCCCCCCCCCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATg                                                                                                                <  2:435109/146‑1 (MQ=255)                                      tgtCTTCC‑GAATTTCCGCCCAGGGGTGTGGCCCCCCATACCCCATCC‑CCGGGGGTTGGC‑TTT‑CCCGACCCGCCAAGCGGGCGCAGCAGG‑GTCTTTTTA‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                               <  2:468891/149‑1 (MQ=255)                                           tcttaattttCAG‑CCAGGGGTGTGGTGCCCGCTAACACATCG‑CCTGGGGTTGGCGTTC‑CCAGCCCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTC‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                           <  2:82900/141‑1 (MQ=255)                                             cTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                      >  2:106517/1‑151 (MQ=255)                                              tGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                     >  2:484641/1‑151 (MQ=255)                                               gAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    >  2:649572/1‑151 (MQ=255)                                                gaatcccccccagcGGTGGGGCCCCCCCAAACACCGCC‑CCTGGGGTGGGCTTTC‑CCAGCCCCGCC‑AGCGGCCGCAGAGGG‑TTATTTTCC‑CCACCCCCGCTGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                      <  1:384533/140‑1 (MQ=255)                                                attttCAG‑CCAGGGGTGTGGTCCCCGCAAACACACCG‑CCTGGGTTTGGCTTTC‑CCAGCCCGCCC‑AGCGGGCGCAGCGGG‑GTATTTTTC‑CCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    <  2:1027950/148‑1 (MQ=255)                                                aaTTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                   >  1:137410/1‑151 (MQ=255)                                                    tCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                >  1:639552/1‑150 (MQ=255)                                                         ccgggggTGGGGTCCCCCCAACCACCCCG‑CCGGCGGTGGCCTTTC‑CCAGCCCTGCA‑AGCGGGCGCGGAGGG‑GTCTTTTGG‑CCCCCCCCGCGCCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGa                                                                                            <  1:370412/145‑1 (MQ=255)                                                          cgggggTGTGGTCCCCACTACCCCACCG‑CCGGCGTTGGGCTTTT‑CCAGACCTCCC‑GGCAGGCGCGGCGGG‑GTATTTTTC‑CCCCCCCCGCCGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    <  2:401629/137‑1 (MQ=255)                                                          cagcGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                          >  2:29282/1‑151 (MQ=255)                                                              ggggtggggCCCCCCCACACACCGCCCTGGGGGTGGCTTTTCCCACCCCCGCC‑AGCAGGCCCA‑CAGGTGTACTTTCC‑CCGCCCCCGCCCCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                             <  1:198435/143‑1 (MQ=255)                                                              ggTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                      >  2:752746/1‑151 (MQ=255)                                                              ggTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                      >  1:154613/1‑151 (MQ=255)                                                                 gtggtgCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                   >  2:268535/1‑151 (MQ=255)                                                                   ggtgCCCGCTAACCCCCCG‑CCTGGGTTTGGCTTTC‑CCAGCCCCCCC‑AGCGGCCGCAGCGGG‑TTACTTTCC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                    <  2:359472/148‑1 (MQ=255)                                                                       cccGCTAACACATCG‑CGTGGGTTGGGCTTTC‑CCAGCCCTGCC‑AGCAGGCGCAGCAGG‑GTTCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg                                                                              <  2:249986/150‑1 (MQ=255)                                                                       cccGATAACCCACCG‑CCGGGGGTGGGCTTTC‑CCAGCCCTGCC‑AGCGGGCGCAGCGGG‑TTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc                                                                               <  2:740326/149‑1 (MQ=255)                                                                          gCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt                                                                          >  2:550134/1‑151 (MQ=255)                                                                             aaCACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                       >  1:24072/1‑151 (MQ=255)                                                                             aaCACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                       >  1:892392/1‑151 (MQ=255)                                                                              accccccccccggCGTTTGGCTTTC‑CCAGCCCTCCC‑AGCGGGCGCAGCGGG‑TTATTTTCC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                        <  1:47481/138‑1 (MQ=255)                                                                              acacATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                      >  1:394783/1‑151 (MQ=255)                                                                                  aTCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                  >  1:1056043/1‑151 (MQ=255)                                                                                        cTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                             >  2:448926/1‑151 (MQ=255)                                                                                          gCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                           >  2:671167/1‑151 (MQ=255)                                                                                            gttgggCTTTC‑CCAGCCCGGCC‑AGCAGGCGCAGCAGG‑GTACTTTTC‑CCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                  <  2:623119/138‑1 (MQ=255)                                                                                              tgggCTTTC‑CCAGCCCTGCC‑AGCGGGCGCGGAGGG‑GTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt                                                          <  2:404587/146‑1 (MQ=255)                                                                                                ggCTTTC‑CCAGCCCTGCC‑AGCGGGCGCAGCAGG‑GTTCTTTTC‑CCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                     <  2:389433/151‑1 (MQ=255)                                                                                                       cccggccccgCC‑AGGGGGCGCAGCGGG‑TTTCTTTCC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                <  1:576271/141‑1 (MQ=255)                                                                                                        ccagccccccc‑AGCGGGCGCAGAAGG‑GTTTTTTTC‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGCTTGCTTACTGCATTGAGCAACAAt                                                <  1:619433/140‑1 (MQ=255)                                                                                                        ccAGCCCGGCC‑AGCGGGGGCAGCGGG‑TTATTTGCG‑CCCCCCCGGCGGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGAAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                <  1:632896/149‑1 (MQ=255)                                                                                                          aGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                            >  1:34891/1‑151 (MQ=255)                                                                                                             ccgccc‑AGCGGGCGCGGCAGG‑TTCCTTTGC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                          <  2:140681/146‑1 (MQ=255)                                                                                                              ccgccccgcggccgcAGCGGG‑GTCCTTTTC‑CCCCCCCCCCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGc                                         <  2:175596/139‑1 (MQ=255)                                                                                                                 cc‑AGCGGGCGAAGCAGG‑TTCCTTTCC‑CCCCCCCCGCGGCCCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                      <  1:612804/150‑1 (MQ=255)                                                                                                                  ccggcgggCGCAGCGGG‑GTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                      <  1:253042/144‑1 (MQ=255)                                                                                                                    agcagGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                   >  1:804279/1‑151 (MQ=255)                                                                                                                     gcagGCGCAGCGGG‑GTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                   <  2:249009/150‑1 (MQ=255)                                                                                                                        ggcgcggaggg‑gTATTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                   <  1:845333/138‑1 (MQ=255)                                                                                                                        ggCGCAGCAGG‑TTTCTTTTC‑CCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                <  2:1057044/150‑1 (MQ=255)                                                                                                                            cagcagG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  2:852072/1‑150 (MQ=255)                                                                                                                            cagcagG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  1:14652/1‑150 (MQ=255)                                                                                                                             agcagG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  2:753749/1‑149 (MQ=255)                                                                                                                                    gttattttGCCCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa                     <  1:1091049/147‑1 (MQ=255)                                                                                                                                          ttccccccccccgccgCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                <  1:293568/140‑1 (MQ=255)                                                                                                                                          ttcccccccaccGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                   <  1:1083723/140‑1 (MQ=255)                                                                                                                                           tgcccacccCCGCTCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                  <  1:810143/144‑1 (MQ=255)                                                                                                                                             cccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg            <  2:645248/144‑1 (MQ=255)                                                                                                                                              ccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc             <  2:1046397/143‑1 (MQ=255)                                                                                                                                                accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg         >  1:73289/1‑151 (MQ=255)                                                                                                                                                accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg         >  2:507115/1‑151 (MQ=255)                                                                                                                                                  caccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc       >  1:985356/1‑151 (MQ=255)                                                                                                                                                    ccgccgCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct      <  1:211029/150‑1 (MQ=255)                                                                                                                                                    ccGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctc     <  1:385340/151‑1 (MQ=255)                                                                                                                                                     cGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc    >  1:97424/1‑151 (MQ=255)                                                                                                                                                      gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct   >  2:1048284/1‑151 (MQ=255)                                                                                                                                                       cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt  >  1:468583/1‑151 (MQ=255)                                                                                                                                                       cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt  >  1:762031/1‑151 (MQ=255)                                                                                                                                                            |                                                                                                                                                  CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993764‑994058

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