Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,917	0	A	C	49.3%	97.1 / 32.0	75	V43G (GTG→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (33/5);  new base C (0/37);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 2.79e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993770‑994063                                                                                                                                                       |                                                                                                                                                    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                                      >  1:1139643/1‑150 (MQ=255)                                                                                                              agGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      >  1:82436/1‑150 (MQ=255)                                                                                                              agGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      >  2:213174/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                                 >  1:983395/1‑149 (MQ=255)                                                                                                                    agcagGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                 >  2:740960/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)                                                                                                                                    cccccccc‑GCTCCGCCCCCCCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc                  <  1:1158575/143‑1 (MQ=255)                                                                                                                                    ccaccaCC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg                 >  2:704945/1‑150 (MQ=255)                                                                                                                                     caccacC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTg               >  1:702220/1‑151 (MQ=255)                                                                                                                                               tGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc      >  1:6768/1‑151 (MQ=255)                                                                                                                                                   gCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  >  2:442838/1‑151 (MQ=255)                                                                                                                                                   cccccccA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa     <  1:452856/147‑1 (MQ=255)                                                                                                                                                    ccccccT‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa     <  2:1035327/147‑1 (MQ=255)                                                                                                                                                    cccACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaaa   <  1:796537/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                    AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993770‑994063

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