Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,904	0	A	C	48.0%	128.0 / 17.0	77	G47G (GGT→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (32/8);  new base C (0/37);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 1.68e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTG  >  NC_000913/993756‑994050                                                                                                                                                       |                                                                                                                                                     ccGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑Ccaccac                                                                                                                                                     <  2:238441/150‑1 (MQ=255)  cGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑g                                                                                                                                                  <  2:373497/151‑1 (MQ=255)    aTTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCt                                                                                                                                                <  1:1109449/151‑1 (MQ=255)        ccACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTg                                                                                                                                               <  1:840049/148‑1 (MQ=255)        ccACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCgg                                                                                                                                            >  1:1139677/1‑151 (MQ=255)         cACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGc                                                                                                                                           >  1:44823/1‑151 (MQ=255)          aCGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGcc                                                                                                                                          <  2:1063476/151‑1 (MQ=255)               aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCccacc                                                                                                                                      <  2:128888/150‑1 (MQ=255)                 cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCacca‑cc                                                                                                                                  >  1:1171854/1‑151 (MQ=255)                    aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGcc                                                                                                                               >  2:452853/1‑151 (MQ=255)                       cgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCca                                                                                                                              >  1:963912/1‑149 (MQ=255)                            ccTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTg                                                                                                                       >  2:793793/1‑151 (MQ=255)                               gAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGAcc                                                                                                                    >  2:520643/1‑151 (MQ=255)                                      catCCGCGTGTCTCCCTAAATTCCAGCCGGGGGTGGGGTCCCCGCAAACACACCGCCGGGGGTGGGCTTTC‑CCACCCCGGCCAGCGGGCGCAGCGGGGTCTTTTCG‑CCCCCCCC‑GCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCgg                                                                                                              <  2:61806/150‑1 (MQ=255)                                        ttCGCGTTTCTCCTTAATTTTCACCCGGCGGTGGGGTCCCCGCAAACACACCGCCGGGGTTTGGCTTTC‑CCACCCCTGCCAGCGGGCGCAGCGGGGTCCTTTTG‑CCCCCCCC‑GCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGcc                                                                                                                <  1:1007197/146‑1 (MQ=255)                                           ggttttCTTCCTAATTTCCAGCCGGCGGTGGGGTGCCCGCAACCACTTCCCCGGGGGTGGGCTTTT‑CCACCCCGCCAAGGGGGCGCAGCAGGTTTTTTTTC‑CCCCCCCC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAat                                                                                                           <  1:825779/143‑1 (MQ=255)                                            cgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGt                                                                                                       >  1:453136/1‑151 (MQ=255)                                             gtgtCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGt                                                                                                       >  2:1124218/1‑150 (MQ=255)                                              tttctcCCTGAATTTCAGCCGGGGGGGGGGTCCCCGCAAACCCACCCCCGGGGGTGGGCTTTC‑CCAGCCCTCCCACCAGGCGCAGAAGGTTCTTTTTC‑CCCCCCCC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                      <  1:253834/148‑1 (MQ=255)                                              tgtctcccgaattttCAGCCGGCGGTGGGGCCCCCCAAAACCCCTCGCCGGCGGTGGGCTTTC‑CCGGCCCTCCCAGCAGGCGCAGAGGGTTACTTTCC‑CCACCCCC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAatat                                                                                                         <  2:31785/135‑1 (MQ=255)                                                tctcCCTAATTTCCAGCCGGGGTTGTGGTCCCCGCAAACACCCCGCCTGGGGTGGGCTTTC‑CCAGCCCGCCCAGCGGGCGCAGCGGGTTTCTTTGC‑CCCCCCCC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                     <  1:1064086/149‑1 (MQ=255)                                                  ttCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                 >  1:1016610/1‑151 (MQ=255)                                                   tCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GTTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                 >  1:893959/1‑150 (MQ=255)                                                   cccaaaattagagccagcGGTGGGGCCCCCGTTAACACGTCCTTGGGGGTTGGCTTTC‑CGAGCCCGCCCAGCAGGCGCAGCGGGTTTTTTTTC‑CCCCCACC‑GCGGGGCCCCCCCC‑CCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                 <  1:1167753/139‑1 (MQ=255)                                                           tccacccgggggTGTGGTCCCCCCAAACCCACCGCCTGGGTTTGCCTTTC‑CCAGCCCCGCCAGCAGGCGCAGCGGGTTCCTTTCC‑CCCCCACC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                           <  2:108366/138‑1 (MQ=255)                                                             cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                        >  1:914853/1‑149 (MQ=255)                                                               gccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                     >  2:196623/1‑150 (MQ=255)                                                                ccggcggggggggccccGCTACCACACCGCCGGGGGTGGGCTTTC‑CCGGCCCGGCCAGCGGGCGCGGCAGGTTCTTTTTG‑CCCCCCCC‑GCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAg                                                                                      <  2:1107725/134‑1 (MQ=255)                                                                     ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                              >  1:332257/1‑151 (MQ=255)                                                                     ggTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                              >  2:1179047/1‑151 (MQ=255)                                                                     ggTGTGGCCCCCCCAAACACCCCGCCTGCGGTTGGCTTTC‑CCGGCCCTGCCGGCGGGCGCAGGAGGTTTTTTTCC‑CCCCCCCC‑GCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                              <  1:446589/151‑1 (MQ=255)                                                                        gtggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                           >  2:582783/1‑151 (MQ=255)                                                                        gtggtcCCCGCTAACCCATCGCCTGCGGTGGGCTTTC‑CCAGCCCCGCCAGCAGGCGCGGCGGGTTTTTTTTC‑CCCCCCCC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                            <  1:985762/150‑1 (MQ=255)                                                                              cccGCAAACACATCGCCTGCGTTTGGTGTTC‑CCACCCCCGCCAGCGGGCGCAGCGGGTTCTTTTTC‑CCCCCCCC‑GCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc                                                                       <  2:984518/149‑1 (MQ=255)                                                                              cccGCAAACACAGCGCCGGGGGTGGGCTTTC‑CCAGCCCGCCCAGCGGGCGCGGAGGGTTTTTTTTC‑CCCCCCCC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc                                                                       <  1:453577/149‑1 (MQ=255)                                                                                cccaaaCACACCGCCGGGGTTGGGCTTTC‑CCAGCCCTCCCAGCAGGCGCAGCAGGTTTTTTTTG‑CCCCCCCC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTaa                                                                   <  1:500851/147‑1 (MQ=255)                                                                                cGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTa                                                                    >  1:649952/1‑150 (MQ=255)                                                                                   aacccccccGCCTGGGGTGGGCTTTC‑CCGGCCCGGCCGGCAGGCGCGGCGGGTTTCTTTTC‑CCCCCCCC‑GCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt                                                                  <  2:816918/141‑1 (MQ=255)                                                                                      cccaccGCCGGCGGTTGGCTTTT‑CCAGCCCTGCCAGCAGGCGCAGCAGGTTTCTTTTC‑CCACCCCC‑GCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                              <  1:440806/145‑1 (MQ=255)                                                                                       acaTCGCCTGCGTTTGGCTTTC‑CCAGCCCTGCCAGCGGGCGCAGCGGGGTACTTTTG‑CCCCCCCC‑GCGGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                              <  1:814563/149‑1 (MQ=255)                                                                                        ccccccggggggTTGCGGTTC‑CCAGGCCGGCAGGCAGGGGCAGAGGGGTGCTTTTC‑CCCCCCCC‑CCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                            <  2:534072/140‑1 (MQ=255)                                                                                        caccGCCTGGGGTGGGCTTTC‑CAAGCCCTCCCAGCGGGCGCAGCAGGTTTCTTTTC‑CCCCCCCC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                 <  2:252431/142‑1 (MQ=255)                                                                                        caccGCCGGGGGTTGGCGTTC‑CCAGCCCTGCCAGCGGGCGCAGGAGGTTTTTTTTC‑CCCCCCCC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCt                                                           <  1:653744/148‑1 (MQ=255)                                                                                         aTCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                          >  2:351960/1‑151 (MQ=255)                                                                                           cGCCTG‑GGTTGGGGTTTCCCGGCCCTCCCGGCGGGCGCAGCAGGTTTTTTTGC‑CCACCCCC‑GCGGGGCCCCCCCT‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTa                                                         <  1:235767/150‑1 (MQ=255)                                                                                             ccTGGGGTTGGCTTTT‑CCAGCCCTCCCAGCAGGCGCAGGGGGGTATTTTCC‑CCCCCCCC‑GCTGCGCCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTa                                                         <  1:340065/148‑1 (MQ=255)                                                                                               tGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                     >  1:936009/1‑150 (MQ=255)                                                                                               ggCGTTTGGTTTTC‑CCAGCCCCCCCAGCGGGCGCAGCGGGTTATTTTCC‑CCCCCCCC‑GCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                      <  1:685185/148‑1 (MQ=255)                                                                                                  gtttGGCGTTC‑CCACCCCGGCCAGCAGGCGCAGGAGGGTACTTTCG‑CCCCCCCC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCa                                                    <  1:1000900/146‑1 (MQ=255)                                                                                                  gttgggCTTTC‑CCGGACCGGCCAGAAGGCGCAGCAGGGTTTTTTCG‑CCCCCACC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                   <  2:387732/145‑1 (MQ=255)                                                                                                      ggCTTTC‑CCAGCCCGGCCAGCGGGCGCAGCGGGTTATTTTTC‑CCCCCCCC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                <  2:827589/148‑1 (MQ=255)                                                                                                       ggctttcccaccccTGCCGGCGGGCGCAGCGGGGTTTTTTTC‑CCCCCCCG‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                              <  2:44823/138‑1 (MQ=255)                                                                                                           tt‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                        >  1:893559/1‑151 (MQ=255)                                                                                                            tccgaccccccccAGCAGCCGCCGCAGTGTCCTTTGC‑CCCCCCCC‑GCAGGGCCCCCCCT‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                            <  2:1148226/136‑1 (MQ=255)                                                                                                            t‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg                                       >  1:536248/1‑151 (MQ=255)                                                                                                             cccAGCCCGGCCAGCAGGCCCAGCAGGGTACTTTCC‑CCCCCCCC‑GCGGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg                                       <  2:380126/150‑1 (MQ=255)                                                                                                                  cccggCCAGCGGGCGCAGGGGGTTTTTTTTC‑CCCCCCCG‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                           <  2:1166492/138‑1 (MQ=255)                                                                                                                  aCCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                  >  1:115101/1‑151 (MQ=255)                                                                                                                   ccgcccAGCGGGCGCGGAGGGGTA‑TTTTG‑CCCCCCCCGGGGGCGGCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGt                                   <  1:701758/145‑1 (MQ=255)                                                                                                                   ccTGCCAGCGGGCGCAGCGGGGTTTTTTTC‑CCCCCCCC‑GCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGt                                   <  1:797530/149‑1 (MQ=255)                                                                                                                     tGCCGGGAGGCGCAGCGGGGTTCTTTGC‑CCACCCCC‑GCGGCGCCCCCCAA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                    <  1:713886/146‑1 (MQ=255)                                                                                                                        cagcagGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTgg                            >  2:276049/1‑151 (MQ=255)                                                                                                                            gggCGCAGCGGGGTTCTTCCC‑CCCCCCCC‑GCGGCGCCCCCCCC‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                           <  1:1039921/147‑1 (MQ=255)                                                                                                                            agGGGCAGCAGGGTCCTTTTC‑CCCCCCCC‑GCTG‑GGCCCCCCAGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGtt                          <  1:631303/149‑1 (MQ=255)                                                                                                                            agGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                         >  2:213174/1‑150 (MQ=255)                                                                                                                            agGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                         >  1:82436/1‑150 (MQ=255)                                                                                                                            agGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                         >  1:1139643/1‑150 (MQ=255)                                                                                                                             ggCGCAGGGGGTTTCTTTTC‑CCCCCCCC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                             <  1:915826/145‑1 (MQ=255)                                                                                                                             ggCGCAGAAGGGTTTTTTTC‑CCCCCCCG‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                              <  1:1001775/144‑1 (MQ=255)                                                                                                                              ggggcagcgggTTCTTTTGCCCACCCCC‑GCGGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                           <  2:110757/139‑1 (MQ=255)                                                                                                                               cgcAGCGGGTTTCTTTTC‑CCCCCCCC‑GCTGCGGCCCCCCC‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                        <  2:1022124/148‑1 (MQ=255)                                                                                                                               cgcAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                     >  2:837450/1‑151 (MQ=255)                                                                                                                                  agcagGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                    >  2:740960/1‑149 (MQ=255)                                                                                                                                  agcagGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                    >  1:983395/1‑149 (MQ=255)                                                                                                                                    cagGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                 >  1:847201/1‑150 (MQ=255)                                                                                                                                         tACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt           >  1:883103/1‑151 (MQ=255)                                                                                                                                            ttttC‑CCCCCCCC‑GCTGCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt            <  1:849587/147‑1 (MQ=255)                                                                                                                                              ttG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc      >  1:856887/1‑151 (MQ=255)                                                                                                                                                  ccaccaCC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg    >  2:704945/1‑150 (MQ=255)                                                                                                                                                   caccacC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg  >  1:702220/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                     CCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTG  >  NC_000913/993756‑994050

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