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

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993770‑994058                                                                                                                                                      |                                                                                                                                              aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCccacca                                                                                                                                            <  2:863422/151‑1 (MQ=255)      aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGcc                                                                                                                                       >  2:991646/1‑151 (MQ=255)      aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCAGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGcc                                                                                                                                       >  2:529396/1‑151 (MQ=255)        gcgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCcac                                                                                                                                     >  1:588776/1‑151 (MQ=255)         cgcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCcac                                                                                                                                     >  2:646246/1‑150 (MQ=255)          gcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCcaca                                                                                                                                    >  2:789430/1‑150 (MQ=255)           cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACaa                                                                                                                                   >  1:379550/1‑150 (MQ=255)           cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACaa                                                                                                                                   >  2:589682/1‑150 (MQ=255)                  aaaCATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCt                                                                                                                           >  1:91710/1‑151 (MQ=255)                     catcatTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGcc                                                                                                                        >  2:294565/1‑151 (MQ=255)                       tcatTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                      >  1:895943/1‑151 (MQ=255)                         tttCGCGTTTCTTCCGGATTTTCAC‑CCGGGGGTGGGGCCCCCCCTACCCCCCCGCCGGCGGTGGGCTTTCCCGCCCCGGCCGGCGGCCCCCGCGGGGTTC‑TTTTCCCCCCCCCGCTGCGCCCCCCCTCCGCCACAAACTGACCTGCCg                                                                                                                       <  1:366947/147‑1 (MQ=255)                         atTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGa                                                                                                                     >  2:919022/1‑150 (MQ=255)                             gcgTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                >  1:351485/1‑151 (MQ=255)                              cgTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                               >  2:335091/1‑151 (MQ=255)                                tgtCTCCCTAATTTCCAG‑CCAGGGGTGTGGCCCCCGCAAACACATCGCCGGCGGTTGGCTTTCCCAGCCCCCCCAGCGGGCGCAGCGGGTTAT‑TTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                  <  2:468499/146‑1 (MQ=255)                                    ttCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                         >  2:712159/1‑151 (MQ=255)                                    ttCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                         >  1:538682/1‑151 (MQ=255)                                     tCCTTAATTTCCGCCCGGGGGTGGGGGCCCCGCAAACACACCGCCTGGGGTGGGCTTTCCCAGCCCTGCCGGCGGGCGCAGCGGGGTTC‑TTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                             <  2:130224/147‑1 (MQ=255)                                      ccTGATTTTCAG‑CGGGGGGTGGGGCCCCCGCTAACACATCGCCGGGGGTTGGCTTTCCCAGCCCTGCCGGCGGGCGCAGCAGGTTTT‑TTTTGCCCCCACCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAg                                                                                                       <  1:314470/151‑1 (MQ=255)                                      ccTGAATTTCAG‑CCAGGGGTGTGGTCCCCGCTAACCCACCGCCTGCGGTTGGCTTTCCCAGCCCTCCCAGCAGGCGCAGCGGGGTCC‑TTTTGCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                           <  2:695448/147‑1 (MQ=255)                                         gAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    >  2:929062/1‑151 (MQ=255)                                            tttCAG‑CCAGCGGTGTGGTCCCCCCTACCACTCCGCCTGCGGTTGGCTTTCCCACCCCTGCCAGCAGGCGCAGCGGGTTAT‑TTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                     <  2:820468/147‑1 (MQ=255)                                             ttccacccgggggggtggtgCCCGCTACCCCCCCGCCTGGGTTTGGCGTTCCCAGCCCGGCCAGCAGGCCCAGCAGGGTCT‑TTTTCCCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                  <  1:315762/136‑1 (MQ=255)                                                   ccgggggTGTGTTCCCCGCTACCACCCCGCCGGGGGTGGGCTTTCCCAGCCCCGCCAGAAGGCCCAGCGGGGTAC‑TTTGCCCCCCCCGGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                             <  2:46682/144‑1 (MQ=255)                                                        ggTGGGGTCCCCGCTAACACATCGCCGGGGGTTGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCGGGTTTT‑TTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGata                                                                                       <  2:588838/150‑1 (MQ=255)                                                             ggtcCCCGCTAACACACCGCCGGGGTTTGGCGTTTCCAGCCCTGCCAGCGGGCGCGGCAGGGTAC‑TTTTCCCCCCCCCGCGGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTc                                                                                  <  2:356664/150‑1 (MQ=255)                                                               tgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc                                                                               >  1:81435/1‑151 (MQ=255)                                                                 cccGCTAACCCACCGCTGGGGTTTGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCAGGGTTT‑TTTCCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                   <  2:946578/145‑1 (MQ=255)                                                                 cccGCAAACACCCCCCCGGGGGTGGGCTTTCCCGGCCCTGCCAGCGGGCGCGGCAGGGTTC‑TTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc                                                                               <  1:538061/149‑1 (MQ=255)                                                                         cacaTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                     >  1:417272/1‑151 (MQ=255)                                                                          acacCGCCGGGGGTTGGCTTTCCCAGCCCGCCCAGCGGGCGCAG‑AGGGTTCTTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                     <  2:567958/150‑1 (MQ=255)                                                                          acaTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                    >  1:601379/1‑151 (MQ=255)                                                                           caccgccggggttgggCTTTCCCAGCCCTGCCAGCAGGCGCAGCAGGTTTT‑TTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                      <  2:954006/134‑1 (MQ=255)                                                                                cgggggttgggCTTTCCCAGCCCGCCCGGCAGGCCCGGAAGGGTTC‑TTTGCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTa                                                                 <  2:1718/139‑1 (MQ=255)                                                                                ccTGCGTTTGGCTTTCCCAGCCCTGCCAGCAGGCGCAGCAGGTTCT‑TTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACt                                                               <  2:544618/150‑1 (MQ=255)                                                                                  tGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                             >  1:483755/1‑150 (MQ=255)                                                                                   gCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                           >  2:491209/1‑151 (MQ=255)                                                                                    cGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt                                                          >  1:477614/1‑151 (MQ=255)                                                                                        gggttttcccAGCCCTCCCGGCGGGGGCGGCGGGGTTT‑TTTTGCCCCCCCCGGGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                  <  2:634824/131‑1 (MQ=255)                                                                                         ggCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                      >  1:496811/1‑150 (MQ=255)                                                                                           ctttcccacccccgCCAGCAGGCGGAGCGGGGTCT‑TTTTCCCCCCCCCGCTGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                              <  1:583427/127‑1 (MQ=255)                                                                                            gTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                   >  1:615388/1‑150 (MQ=255)                                                                                             ttcccAGCCCTGCCAGCGGGCGCAGCAGGGTTC‑TTTTGCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                     <  2:425260/144‑1 (MQ=255)                                                                                              ttCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                >  2:522741/1‑151 (MQ=255)                                                                                              tccaagcccccccAGCAGGCGCAGGAGGTTCT‑TTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                 <  2:422890/139‑1 (MQ=255)                                                                                                ccAGCCCTGCCAGCAGGCGCAGCAGGTTTT‑TTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                <  2:114537/149‑1 (MQ=255)                                                                                                   gACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGt                                           >  1:394198/1‑151 (MQ=255)                                                                                                        gCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                      >  1:309013/1‑151 (MQ=255)                                                                                                         ccAGCAGGCGCAGCAGGGTAC‑TTTTGCCCCCCCCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGcc                                        <  2:994296/148‑1 (MQ=255)                                                                                                          cagcagGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTgg                                    >  2:521145/1‑151 (MQ=255)                                                                                                               ggCGCAGCGGGGTAC‑TTTTGCCCCCCCCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                <  2:972466/150‑1 (MQ=255)                                                                                                               ggCGCAGCAGGTTAC‑TTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAg                               <  2:86427/151‑1 (MQ=255)                                                                                                               ggCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                 >  1:828721/1‑149 (MQ=255)                                                                                                               ggCGCAGAAGGTTTT‑TTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTgg                                    <  1:110087/146‑1 (MQ=255)                                                                                                                  gcagcaGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  2:578257/1‑151 (MQ=255)                                                                                                                       agGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                         >  2:515423/1‑149 (MQ=255)                                                                                                                       agGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                         >  2:43550/1‑149 (MQ=255)                                                                                                                        gggTCC‑TTTGCCCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                <  2:423037/141‑1 (MQ=255)                                                                                                                        gggTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                      >  1:781467/1‑151 (MQ=255)                                                                                                                               ttttGCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                  <  1:990466/149‑1 (MQ=255)                                                                                                                                    ccccccccGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt           <  1:679794/145‑1 (MQ=255)                                                                                                                                    cccccaccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc             <  1:881758/146‑1 (MQ=255)                                                                                                                                    ccaccaCCGCTGCGGCCCACCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt           <  2:968949/151‑1 (MQ=255)                                                                                                                                    ccaccaCCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg            >  2:423091/1‑150 (MQ=255)                                                                                                                                        caccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc       >  1:301348/1‑151 (MQ=255)                                                                                                                                             cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt  >  2:361958/1‑151 (MQ=255)                                                                                                                                             cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt  >  1:431639/1‑151 (MQ=255)                                                                                                                                             cTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctc     <  2:670728/148‑1 (MQ=255)                                                                                                                                                cGGCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct   <  1:812267/147‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                              AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTAC‑TTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993770‑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: ‑