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

CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCT  >  NC_000913/993764‑994057                                                                                                                                                         |                                                                                                                                                 cACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGc                                                                                                                                                 >  2:596332/1‑151 (MQ=255)       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                                                                                                                       >  1:145002/1‑151 (MQ=255)                           catcatTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGcc                                                                                                                       >  2:647264/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)                                                                                                                                                         |                                                                                                                                                 CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTT‑CAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCT  >  NC_000913/993764‑994057

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