Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,917	0	A	C	46.3%	96.9 / 23.6	67	V43G (GTG→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (32/4);  new base C (0/31);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 5.85e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993770‑994064                                                                                                                                                       |                                                                                                                                                    aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccacc                                                                                                                                                   >  2:420308/1‑150 (MQ=255)              ccTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTg                                                                                                                                     >  1:824737/1‑151 (MQ=255)               cTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTg                                                                                                                                     >  2:299676/1‑150 (MQ=255)                tGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGAc                                                                                                                                   >  2:540337/1‑151 (MQ=255)                 gAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGAcc                                                                                                                                  >  1:784687/1‑151 (MQ=255)                  aaccatcttCCGCGTTTCTTCCTAATTTC‑CGGCCGGGGGGGGGGGCCCCGCAAACACATCGCCTGGGGTTGGCTTTC‑CCAGCCCTGCCAGCAGGCGCAGCAGGGTTTTTTTG‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGAc                                                                                                                                   <  1:943942/146‑1 (MQ=255)                   accatcttTCGCGTGTCTTCCTGAATTT‑CAGCCGGGGGTGTGGGCCCCGCTACCACATCGCCGGCGGTTGGCGTTT‑CCAGCCCGGCCAGCGGGCGCAGCAGGTTCTTTTTG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTg                                                                                                                                <  2:934580/149‑1 (MQ=255)                     catcatTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGcc                                                                                                                              >  2:647264/1‑151 (MQ=255)                     catcatTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGcc                                                                                                                              >  1:145002/1‑151 (MQ=255)                          ttCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAat                                                                                                                         >  1:663566/1‑151 (MQ=255)                             gcgTGTCTTCCTAAATTT‑CACCCGGGGGTGTGGTCCCCGCTAACACACCGCCTGGGGTTGGCTTTT‑CCAGCCCTGCCAGCGGGCGCAGCGGGGTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGa                                                                                                                           <  2:833633/146‑1 (MQ=255)                                tgtCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                    >  1:25980/1‑150 (MQ=255)                                    ttCCTGATTTT‑CAGCCAGGGGTGTGGTGCCCGCAAACACATCGCCGGGGGTGGGCGTTT‑CCAGCCCTCCCAGCAGGCGCAGCAGGGTATTTTTG‑CCCCCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                                 <  2:807176/149‑1 (MQ=255)                                    ttCCTGATTTC‑CAGCCAGGGGTGGGGTCCCCGTAAACACCCCGCCGGGGGTTGGCTTTC‑CCAGCCCTGCCAGCAGGCGAAGCGGGGTATTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGaa                                                                                                                          <  2:31252/140‑1 (MQ=255)                                     tCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                               >  1:491036/1‑150 (MQ=255)                                      ccgaaaTTC‑CAGCCGGGGGTGTGGTCCCCGCTACCCCCCCGCTGGGGGTTGGCTTTT‑CCACCCCGGCCAGCGGGCGCAGCAGGGTCCTTTTC‑CCCCCCCCGGGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                                <  1:444486/144‑1 (MQ=255)                                       cctaaTTTCCAGCCGGGGGTGTGGTGCCCGCTAACCCACCGCCGGCGGTTGGCTTTC‑CCAGCCCTGCCAGCAGGCGCGGCGGGGTATTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                               <  2:340434/146‑1 (MQ=255)                                          atttt‑CAGCCGGCGGTGTGGTCCCCGCTACCGCACCGCGGGCGGTTGGCTTTT‑CCAGCCCTGCCGGCGGGCGCAGCGGGGTCTTTTGG‑CCACCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                           <  2:849176/147‑1 (MQ=255)                                           aattccccccagcGGTGTGGCCCCCCCAAACCCACCCCCGGCGTTTGGCGTTC‑CCACCCCCCCCAGCAGGCGCGGAGGGTTTTTTTTG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                         <  1:180827/143‑1 (MQ=255)                                             tt‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                      >  2:496248/1‑151 (MQ=255)                                              t‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                      >  1:692885/1‑150 (MQ=255)                                               gcagccgggGGTGTGTCGCCCCCTAAACCCTCGCCGGGGTTGGGCTTTC‑CCAGCCCTGCCAGCAGGCGCAGGGGGTTCTTTTTC‑CCCCCCCCGCTGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCa                                                                                                     <  1:1023341/150‑1 (MQ=255)                                                   ccgggggggggggccccGCAAACCCATCGCCTGCGTTTGGCTTTC‑CCAGCCCTGCCAGCGGGCGCAGCAGGTTATTTTTG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGa                                                                                                  <  1:870127/136‑1 (MQ=255)                                                     agcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                                >  1:983638/1‑150 (MQ=255)                                                        ggtggggcccccGCAAACACATCGCCTGCGGTTGGCTTTC‑CAAGCCCTGCCAGCAGGCGCAGCAGGTTATTTTTG‑CCCCCACCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGata                                                                                             <  2:31821/141‑1 (MQ=255)                                                          tgtgGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                          >  2:378742/1‑151 (MQ=255)                                                             ggtcCCCGCAACCACATCGCCGGGGTTGGGCTTTC‑CCAGCCCTGCCAGCGGGCGCAGCGGGTTATTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                       <  2:42351/151‑1 (MQ=255)                                                                       aaCACACCGCCTGCGGTTGGCTTTC‑CCAGCCCTCCCAGCGGGCGCAGCAGGTTATTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt                                                                                <  2:646510/148‑1 (MQ=255)                                                                         cacaTCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                           >  1:762140/1‑151 (MQ=255)                                                                         cacaTCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                           >  1:498937/1‑151 (MQ=255)                                                                         cacaCCCCCTGCGGTGGGCTTTC‑CCAGCCCGCCCAGCGGGCGCGGCGGGTTCTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                           <  1:388691/151‑1 (MQ=255)                                                                            aTCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                        >  1:96646/1‑151 (MQ=255)                                                                             tCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                        >  2:715344/1‑150 (MQ=255)                                                                               gCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACt                                                                     >  1:256587/1‑151 (MQ=255)                                                                                 ccgggggtGGCCTTTCCCAGCCCGGCCGGCGGGCGCAGCAGGGTCTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                            <  2:528885/136‑1 (MQ=255)                                                                                   gggttgggCTTTC‑CCAGCCCGCCCAGCGGGCGCAGCGGGGTCTTTTTC‑CCACCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                        <  1:968596/138‑1 (MQ=255)                                                                                   gCGTTTGGCGTTC‑CCACCCCTGCCAGCAGGCGCAGAGGGGTACTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                   <  2:723691/149‑1 (MQ=255)                                                                                     gttgggCTTTC‑CCAGCCCGCCCAGCGGGCGCAGCGGGTTACTTTTC‑CCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                                   <  1:125263/143‑1 (MQ=255)                                                                                     ggTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                               >  1:157425/1‑151 (MQ=255)                                                                                     ggTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                               >  1:675936/1‑151 (MQ=255)                                                                                       tgggCTTTC‑CCAGCCCCCCCAGCAGGCGCAGCAGGTTTTTTTCC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                              <  1:658059/148‑1 (MQ=255)                                                                                         ggCTTTT‑CCAGCCCTGCCAGCGGGCGCAGCGGGGTATTTTTG‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                            <  2:712283/150‑1 (MQ=255)                                                                                          ggctttcccaccccTCCCAGCAGGCCCAGCGGGTTCCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                            <  1:438863/138‑1 (MQ=255)                                                                                                  cAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                                   >  2:356699/1‑151 (MQ=255)                                                                                                   aGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                                  >  1:900787/1‑151 (MQ=255)                                                                                                      ccTGCCAGCAGGCGCAGCGGGTTATTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                                <  2:794956/150‑1 (MQ=255)                                                                                                        tGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGccc                                             >  2:421149/1‑151 (MQ=255)                                                                                                         gCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                            >  1:404596/1‑151 (MQ=255)                                                                                                              cagGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                       >  2:781367/1‑151 (MQ=255)                                                                                                              cagGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                       >  2:832441/1‑151 (MQ=255)                                                                                                                ggCGCAGCAGGGTCCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                       <  2:836322/149‑1 (MQ=255)                                                                                                                gccgcAGCGGGGTCTTTTTG‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                      <  2:743034/148‑1 (MQ=255)                                                                                                                   gcagcaGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                  >  1:678779/1‑151 (MQ=255)                                                                                                                    cagcagGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                  >  1:146030/1‑150 (MQ=255)                                                                                                                      gcagGTTCCTTTTC‑CCCCCCCCGCTGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                                    <  2:45911/146‑1 (MQ=255)                                                                                                                       caggtttttttcG‑CCCCCCCCGCTGCGGCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                                    <  1:81301/137‑1 (MQ=255)                                                                                                                        agGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                               >  1:891356/1‑149 (MQ=255)                                                                                                                         ggttatttttG‑CCCCCCCAGGTCCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                              <  2:676220/143‑1 (MQ=255)                                                                                                                          ggTTCTTTTGCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                   <  2:1009489/144‑1 (MQ=255)                                                                                                                            tctttttC‑CCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa                           <  2:288299/146‑1 (MQ=255)                                                                                                                                tttgcccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCa                       <  2:971086/139‑1 (MQ=255)                                                                                                                                tttccccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                        <  1:713505/138‑1 (MQ=255)                                                                                                                                      cccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg                <  1:544802/146‑1 (MQ=255)                                                                                                                                        ccaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGc              >  1:594807/1‑151 (MQ=255)                                                                                                                                          cccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct            <  2:678844/150‑1 (MQ=255)                                                                                                                                          cccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct            <  1:61436/150‑1 (MQ=255)                                                                                                                                            cGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc          >  2:875869/1‑151 (MQ=255)                                                                                                                                             gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct         >  2:426005/1‑151 (MQ=255)                                                                                                                                                   cccccccccCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaa     <  1:1017043/148‑1 (MQ=255)                                                                                                                                                   cccccccccCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt   <  2:145011/150‑1 (MQ=255)                                                                                                                                                    ccccccccCGCCACAAACTGACCTGCCGGAATATGTAAATCCTGGTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa      <  1:469667/147‑1 (MQ=255)                                                                                                                                                    ccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATg  <  1:496206/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                    AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAATG  >  NC_000913/993770‑994064

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