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

GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA  >  NC_000913/993768‑994060                                                                                                                                                        |                                                                                                                                                gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCccac                                                                                                                                                >  1:590718/1‑151 (MQ=255) gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCccac                                                                                                                                                <  2:848230/151‑1 (MQ=255)   aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCccacc                                                                                                                                               <  1:4818/150‑1 (MQ=255)     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCaccacc                                                                                                                                            <  2:907553/151‑1 (MQ=255)     cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCaccacc                                                                                                                                            <  1:832305/151‑1 (MQ=255)                   gAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGAcc                                                                                                                              >  2:208384/1‑151 (MQ=255)                     accatcatCCGCTTTTCTTCCTAATTTCCACCCAGCGGTGGGGTCCCCGTAAACACACCCCCTGGGTTTGGCTTTT‑CCAGCCCGGCCAGCAGGCGCAGCAGGGTATTTTC‑CCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTg                                                                                                                            <  1:743092/149‑1 (MQ=255)                     aaCATCTTCCGCGTTTCTTCCGAAATTTCAGCCGGGGGGGGGGTCCCCGCAACCACACCG‑CGGGGGTTGGCCTTTCCCAGCCCGGCCAGCAGGCGCAGAGGGGTTCTTTT‑GCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTg                                                                                                                            <  2:797880/151‑1 (MQ=255)                       catcatCCGCGTGTCTTCCTGAATTTCAGCCGGCGGTGTGGGGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGCCCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGc                                                                                                                           <  1:716670/150‑1 (MQ=255)                        atcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                         >  1:779366/1‑151 (MQ=255)                         tcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGTTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                        >  2:417066/1‑151 (MQ=255)                           atTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGa                                                                                                                       >  1:280744/1‑150 (MQ=255)                                  tgtCTTCCTAAATTTCAGCCGGCGGTGTGGTGCCCGCAAACACACCGCCTGGGTTGGGCTTTC‑CCAGCCCTGCCGGCGGGCGCGGCGGGTTTCTTTT‑GCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                    <  1:395192/146‑1 (MQ=255)                                   gtCTTCCTGAATTTCAGCCAGCGGCGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                               >  1:206464/1‑150 (MQ=255)                                     cTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                            >  1:307987/1‑151 (MQ=255)                                      ttCCTAATTTCCAGCCGGGGGGGGGGGCCCCGCAACCCCATCCCCTGCGGTTGGCTTTC‑CCAGCCCTGCCGGCAGGCGCAGCAGGGTATTTTT‑GCCCCCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                <  2:290146/146‑1 (MQ=255)                                       tCCTGATTTCCACCCGGGGGTGTGGTCCCCGCTAACACCCCGCCGGGGGTTGGCTTTT‑CCAGCCCTCCCAGCGGGCGCAGCAGGGTATTTTT‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                             <  1:214025/148‑1 (MQ=255)                                         cTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                        >  2:587280/1‑151 (MQ=255)                                           gAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                      >  2:507469/1‑151 (MQ=255)                                              ttccacccgggggTGGGGTCCCCGCTACCCCACCGCCGGCGTTTGGCTTTC‑CCAGCCCTCCCAGCGGGCGCGGCGGGTTTTTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                      <  2:947571/137‑1 (MQ=255)                                                  agccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                               >  2:485608/1‑151 (MQ=255)                                                    ccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                             >  1:635591/1‑151 (MQ=255)                                                         ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                        >  1:408123/1‑151 (MQ=255)                                                         ggTGTGGTCCCCGCTAACACACCGCCGGGGGTTGGCGTTT‑CCAGACCGGCCAGCGGGCGCAGCGGGTTCTTTTT‑GCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGa                                                                                           <  1:541486/148‑1 (MQ=255)                                                           tgggggccccGCTAACACACCGCCTGCGGTGGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCGGGTTCCTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                       <  1:38503/143‑1 (MQ=255)                                                            gtggtcCCCGCTAACACCTCGCCTGGGTTTGGCTTTC‑CCACCCCTCCCAGCAGGCGCAGCGGGGTCCTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                       <  1:712174/149‑1 (MQ=255)                                                              ggtcCCCGCAAACCCACCGCCGGGGGTTGGCTTTC‑CCAGCCCTGCCAGCAGGCGCAGCGGGGTACTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTc                                                                                    <  1:226317/150‑1 (MQ=255)                                                              ggcgCCCGCAAACACACCGCCGGCGTTTGGCGTTC‑CCAGCCCCCCCAGCAGGCGCAGAAGGTTCTTTTT‑CCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                     <  2:1003275/146‑1 (MQ=255)                                                                  cccGCAAACACTCCGCTGGGGTTTGGCTTTC‑CCAGCCCGGCCAGCAGGCGCAGCAGGGTTTTTTT‑GCCCCCACCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg                                                                                <  1:353435/150‑1 (MQ=255)                                                                   ccGCTAACCCATCGCCTGCGGTTGGCTTTC‑CCAGCCCTGCCAGCAGGCGCAGCGGGTTATTTTT‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTa                                                                              <  1:802810/151‑1 (MQ=255)                                                                    cGCAAACCCACCCCCGGCGGTTGGCTTTC‑CCAGCCCGGCCAGCAGGCGCAGCGGGTTCTTTTT‑GCCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                             <  1:6673/151‑1 (MQ=255)                                                                        aaCACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                         >  1:145895/1‑151 (MQ=255)                                                                        aaCACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                         >  2:958899/1‑151 (MQ=255)                                                                             atcgcggggggTGGGCGTTC‑CCAGCCCTCCCGGCGGGCGCAGCAGGGTACTTTT‑GCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                      <  1:179179/140‑1 (MQ=255)                                                                             aTCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                    >  2:688194/1‑151 (MQ=255)                                                                                   ggctgTTGGCTTTC‑CCCGCCCCCCCCGCAGCCGCAGCAGTGTCTTTTC‑GCCGCCCCCGCGGCCCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                               <  2:743892/146‑1 (MQ=255)                                                                                    gCGGTTGGCTTTC‑CCAGCCCTGCCAGCGGGCGCAGCAGGGTACTTTT‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                <  1:748129/148‑1 (MQ=255)                                                                                      ggTTGGCTTTC‑CCACCCCGCCCAGCAGGCGCAGCGGGGTATTTTT‑CCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                             <  1:617210/149‑1 (MQ=255)                                                                                      ggTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                           >  1:557397/1‑151 (MQ=255)                                                                                      ggTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                           >  2:286871/1‑151 (MQ=255)                                                                                      ggTTGCCTTTC‑CCAGCCCCGCCAGGGGGCGCGGGAGGGTTCTTTT‑CCCCCCCCCGGGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                             <  2:786603/149‑1 (MQ=255)                                                                                       gTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                           >  1:415020/1‑150 (MQ=255)                                                                                              ttt‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaaca                                                    >  2:635972/1‑150 (MQ=255)                                                                                                   cAGCCCTCCCAGCAGGCGCAGCGGGGTATTTTT‑GCCCCCCCCGGTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                  <  2:824437/148‑1 (MQ=255)                                                                                                      cccTGCCGGCAGGCGCGGCGGGGTATTTTT‑GCCCCCACCGCGGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                <  1:844160/146‑1 (MQ=255)                                                                                                       ccTGCCAGCGGGCGCAGCAGGGTTTTTTC‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                            <  2:785814/150‑1 (MQ=255)                                                                                                                 ggCGCAGC‑GGGTTCTTTTCCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAg                                 <  1:962468/151‑1 (MQ=255)                                                                                                                   cgcAGCGGGTTTTTTTT‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                        <  1:320866/142‑1 (MQ=255)                                                                                                                   cgcAGCGGGGTTTTTTT‑CCCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                               <  2:280709/151‑1 (MQ=255)                                                                                                                     cggcgggttccTTTT‑CCCACCCCCTCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa                             <  1:988887/141‑1 (MQ=255)                                                                                                                      agcagGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                              >  2:786678/1‑149 (MQ=255)                                                                                                                        cagGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                           >  2:688756/1‑150 (MQ=255)                                                                                                                          ggtttttttt‑CCCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGa                             <  1:316772/140‑1 (MQ=255)                                                                                                                            gTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                      >  1:988932/1‑151 (MQ=255)                                                                                                                            gTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                      >  1:69799/1‑151 (MQ=255)                                                                                                                                tttt‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                  >  2:366788/1‑151 (MQ=255)                                                                                                                                tttt‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                    <  1:507515/149‑1 (MQ=255)                                                                                                                                 tttgcccaccgCCGTTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                    <  1:748210/144‑1 (MQ=255)                                                                                                                                 ttt‑CCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                  <  1:258285/150‑1 (MQ=255)                                                                                                                                      cccccccggCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATa                 <  2:206449/139‑1 (MQ=255)                                                                                                                                      ccccccccGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc               <  2:831542/143‑1 (MQ=255)                                                                                                                                      ccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt             <  2:630210/145‑1 (MQ=255)                                                                                                                                       cccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg            <  1:644483/146‑1 (MQ=255)                                                                                                                                       caccacCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg            >  1:669415/1‑151 (MQ=255)                                                                                                                                        accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg           >  2:14911/1‑151 (MQ=255)                                                                                                                                            ccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct        <  1:139531/150‑1 (MQ=255)                                                                                                                                            ccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctc       <  2:444007/151‑1 (MQ=255)                                                                                                                                             cgcgGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc      <  2:539109/151‑1 (MQ=255)                                                                                                                                               cggCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct     <  1:636023/148‑1 (MQ=255)                                                                                                                                               cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt    >  2:794140/1‑151 (MQ=255)                                                                                                                                                tGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc   >  1:402293/1‑151 (MQ=255)                                                                                                                                                  cgcccccccACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct     <  2:924006/144‑1 (MQ=255)                                                                                                                                                  cgccccACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa  <  2:280727/147‑1 (MQ=255)                                                                                                                                                  cGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa  >  2:274010/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑GCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCA  >  NC_000913/993768‑994060

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