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

ACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993765‑994058                                                                                                                                                       |                                                                                                                                                  aCGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGcc                                                                                                                                                 >  2:1108349/1‑151 (MQ=255)  cGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGccc                                                                                                                                                <  1:978936/151‑1 (MQ=255)    gTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCccac                                                                                                                                              <  2:282352/151‑1 (MQ=255)     tAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCccacc                                                                                                                                             >  1:682802/1‑151 (MQ=255)     tAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCccacc                                                                                                                                             >  1:216506/1‑151 (MQ=255)         aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCaccacc                                                                                                                                          >  1:1287685/1‑150 (MQ=255)            cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGcc                                                                                                                                       >  1:84243/1‑150 (MQ=255)            cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCca                                                                                                                                      <  2:305997/151‑1 (MQ=255)              cgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCcac                                                                                                                                     >  1:1389853/1‑150 (MQ=255)               gcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCcaca                                                                                                                                    >  2:423702/1‑150 (MQ=255)               gcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCcaca                                                                                                                                    >  2:892140/1‑150 (MQ=255)                 aTCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAAc                                                                                                                                 >  2:985376/1‑151 (MQ=255)                  tCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACt                                                                                                                                >  1:1128835/1‑151 (MQ=255)                      gAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGAcc                                                                                                                            >  1:313685/1‑151 (MQ=255)                       aaccatcttTCGGGTGTCTTCCTGATTTCCAGCCGGGGGTGGGGTCCCCGCTAACCCACCGCCGGGGTTGGGCTTCCCCAGCCCGGCCAGCGGGCGCAGCGGGTTCCTTTG‑C‑CCACCCCCGCGGCGCCCCACCACCGCCACAAACTGAc                                                                                                                             <  2:902337/146‑1 (MQ=255)                       aaccacctttCGCTTTTCTTCTTAAATTCCACCCGGGGGTGGGGGCCCCCCAAACACCCCGCCGGGGGTGGGCTTTCCCACCCCGCCCGGCAGGCGCAGCGGGGTCCTTTG‑C‑CCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCt                                                                                                                           <  2:673758/143‑1 (MQ=255)                          catcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGcc                                                                                                                        >  2:1299069/1‑151 (MQ=255)                             ctttCGCGTGTCTCCTTGATTTCCACCCGGGGGTGGGGTCCCCGCAAACACACCGCCGGGGTTTGGCTTTCCCAGCCCGCCCAGCAGGCCCAGCAGGGTTTTTTT‑C‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGcc                                                                                                                        <  1:620428/146‑1 (MQ=255)                              ttccgctTGTCTCCCTGATTTTCACCCGGGGTGGGGGCCCCCGCAACCACACCGCGTGGGTTTGGCGTTCCCACCCCCGCCAGCGGGCGCAGGGGGTTTCTTTC‑G‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTg                                                                                                                          <  2:701701/142‑1 (MQ=255)                                   cgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                >  2:637649/1‑150 (MQ=255)                                   cgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                               >  1:427698/1‑151 (MQ=255)                                     tgtCTTCCTTATTTCCGCCCGGGGGTGTGGCCCCCGCTACCCCTCCCCCGGGGGTGGGCTTTCCCCGCCCCCCCAGCAGCCGCAGCTGGGTATTTTC‑G‑CCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                              <  2:348962/150‑1 (MQ=255)                                      gtCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                             >  1:1203520/1‑150 (MQ=255)                                         ttCCTAAATTCCAGCGGGGGGGGGGGCCCCCGCAAACACCTCGCCTGGGGTTGGCGTTTCCACCCCTGCCAGCCGGCGCAGCGGGGTTTTTTT‑C‑CCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                  <  1:99644/142‑1 (MQ=255)                                          tctgaaaTTCCAGCCGGGGGTGGGGCCCCCGCAAACACCCCGCCGGCGGTTGGCTTTCCCAGCCCCCCCAGCGGGCGCGGCGGTTTTTTTTC‑C‑CCACCACCGGGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                <  2:405697/138‑1 (MQ=255)                                           ccggatttccagccgggGGTGGGGTCCCCGCAAACCCTCCGCCGGGGGTGGGCTTTCCAACCCCTGCCGGCGGGCGCGGCGGGGTCTTTTT‑C‑CCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                 <  1:229433/132‑1 (MQ=255)                                              gAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    >  2:1388796/1‑151 (MQ=255)                                               aaTTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                   >  1:483073/1‑151 (MQ=255)                                                  ttCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                >  1:366063/1‑151 (MQ=255)                                                   tCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                >  1:292686/1‑150 (MQ=255)                                                       ccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                           >  1:666493/1‑151 (MQ=255)                                                        cagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                          >  2:87269/1‑151 (MQ=255)                                                          ggggtgggggccccGCAAACCCACCGCCTGGGTTGGGCTTTCCCAGCCCGGCCGGGGGGCGCAGGAGGGTATTTTG‑C‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                           <  1:826763/137‑1 (MQ=255)                                                          ggggtgggggccccGATAACACACCGCCGGCGGTGGGCTTTCCCACCCCGCCCGGCAGGCGCAG‑AGGGTTTTTTT‑GCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                        <  2:567344/140‑1 (MQ=255)                                                           gggtggggcccccGCAAACCCACCGCGTGGGGTTGGCTTTCCCAGCCCGCCCAGCGGGCGCAGCAGGTTCTTTTC‑C‑CCCCCCCCCCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                           <  2:806258/137‑1 (MQ=255)                                                            ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                      >  2:156889/1‑151 (MQ=255)                                                              tggggcccccGCAAACCCCCCGCCGGCGGTTGGCTTTCCCAGCCCGCCCAGCGGGCGCGGAAGGGTCCTTTT‑C‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                     <  1:129671/143‑1 (MQ=255)                                                                 gggccccGCAACCACCGCGCGTGGGGTTGGCGTTCCAAGCCCTCCCAGAGGCCGCGGGGGGTTCTTTTT‑C‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                     <  1:378296/143‑1 (MQ=255)                                                                     cccGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                             >  2:54836/1‑151 (MQ=255)                                                                     cccGCTAACACATCGCCGGCGTTTGGCTTTCCCAGCCCGGCCAGCGGGCGCAGC‑GGGTTCTTTTCC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc                                                                               <  1:179892/149‑1 (MQ=255)                                                                     cccGCAACCACCCCGCCTGGGGTTGGCGTTTCCAGCCCTGCCAGCGGGCGCAGCGGGTTTCTTTT‑G‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTc                                                                                  <  1:182422/146‑1 (MQ=255)                                                                       cGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTa                                                                            >  1:219752/1‑150 (MQ=255)                                                                        cccacccccccGCCGGGGGTTGGCTTTCCCAGCCCGCCCAGCGGGCGCGGCGGGTTTCTTTC‑C‑CCCCCCCGGCCGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt                                                                          <  2:1388162/141‑1 (MQ=255)                                                                           aaCACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                       >  2:859195/1‑151 (MQ=255)                                                                             cacaCCGCCGGGGGTTGGCTTTCCCAGCCCCGCCAGCGGGCGCGGCGGGTTACTTTC‑C‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                     <  1:1426636/151‑1 (MQ=255)                                                                                 tCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                  >  1:1072796/1‑150 (MQ=255)                                                                                           ttggttttcccAGCCCTGCCAGCGGGCGCAGCGGGTTATTTTT‑C‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                        <  1:859217/141‑1 (MQ=255)                                                                                           tgggttttcccaccccTGCCAGCGGGCCCAGCAGGTTCCTTTC‑C‑CCCCCCCCGCTGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                        <  2:1389812/136‑1 (MQ=255)                                                                                             ggctttcccgccccTGCCGGCAGCCGAGGAAGGGTACTTTT‑C‑CCCCCGCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                     <  1:1001194/139‑1 (MQ=255)                                                                                             ggCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                      >  2:722644/1‑150 (MQ=255)                                                                                             ggCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                      >  1:857079/1‑150 (MQ=255)                                                                                               cGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                   >  1:165887/1‑151 (MQ=255)                                                                                                  ttCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                >  2:153282/1‑151 (MQ=255)                                                                                                   cccaccccTCCCAGGGGGCGCAGCGGGTTTTTTTT‑C‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                 <  1:184550/143‑1 (MQ=255)                                                                                                      accccccccggcagGCGCAGCGGGTTTTTTTC‑C‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                              <  2:429106/139‑1 (MQ=255)                                                                                                       cccccgcgagcggGCGCAGGGGTTTATTTTC‑G‑CCCCCCCGGCTGCCCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                                    <  2:1119114/134‑1 (MQ=255)                                                                                                        ccccgCCAGCAGGCGCAGCGGGGTTTTTTC‑G‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                <  2:847578/141‑1 (MQ=255)                                                                                                          cTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGcc                                        >  2:352097/1‑151 (MQ=255)                                                                                                             ccAGCAGGCGCAGCGGGTTTTTTTT‑G‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGccc                                       <  1:1367156/149‑1 (MQ=255)                                                                                                               agcagGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                   >  2:454866/1‑151 (MQ=255)                                                                                                                   ggCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                 >  1:1164677/1‑149 (MQ=255)                                                                                                                    gctcaggagGTTATTTTA‑C‑CCACCCCCGTTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAg                               <  1:568180/147‑1 (MQ=255)                                                                                                                         gcagGGTTTTTTC‑C‑CCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                              <  2:733130/146‑1 (MQ=255)                                                                                                                            ggttACTTTC‑C‑CCCCCCCCGCTGCCCCCCACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAg                          <  2:1231642/144‑1 (MQ=255)                                                                                                                            gggTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                      >  1:1055600/1‑151 (MQ=255)                                                                                                                             ggggatttt‑GCCCACCACGGCAGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                      <  1:433323/146‑1 (MQ=255)                                                                                                                                    tt‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc              >  2:293614/1‑151 (MQ=255)                                                                                                                                     ttccccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCa                 <  2:216503/139‑1 (MQ=255)                                                                                                                                        cccccccccccggcgcCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg            <  2:1128805/139‑1 (MQ=255)                                                                                                                                         ccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt           <  1:681606/145‑1 (MQ=255)                                                                                                                                         ccccccccGCGGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc             <  1:1078212/143‑1 (MQ=255)                                                                                                                                         ccaccaCCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt           >  1:776389/1‑151 (MQ=255)                                                                                                                                         ccaccaCCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt           >  2:369272/1‑151 (MQ=255)                                                                                                                                         ccaccaCCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt           >  1:803406/1‑151 (MQ=255)                                                                                                                                           ccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg         <  2:991826/147‑1 (MQ=255)                                                                                                                                            cccccGCTGCGGCCCCCAACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg            <  1:985401/144‑1 (MQ=255)                                                                                                                                              accGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct      >  1:401921/1‑151 (MQ=255)                                                                                                                                               ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctc     >  1:926104/1‑151 (MQ=255)                                                                                                                                               ccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct      <  1:90234/150‑1 (MQ=255)                                                                                                                                                  cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt  >  1:709133/1‑151 (MQ=255)                                                                                                                                                  cTGCGCCCCCCATCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCGTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc    <  2:934648/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                  ACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993765‑994058

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