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

CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993763‑994063                                                                                                                                                          |                                                                                                                                                        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                                      >  1:1139643/1‑150 (MQ=255)                                                                                                                     agGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      >  2:213174/1‑150 (MQ=255)                                                                                                                     agGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      >  1:82436/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)                                                                                                                                           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)                                                                                                                                                          |                                                                                                                                                        CCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACC‑GCTGCGGCCCACCA‑CCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993763‑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: ‑