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

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAA  >  NC_000913/993770‑994062                                                                                                                                                       |                                                                                                                                                  aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCccacca                                                                                                                                                <  2:986934/151‑1 (MQ=255)   cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCcaccac                                                                                                                                               <  1:1036988/150‑1 (MQ=255)    aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCaccacc                                                                                                                                              >  2:575826/1‑150 (MQ=255)      aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGcc                                                                                                                                           <  2:51034/151‑1 (MQ=255)       cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGcc                                                                                                                                           >  2:480426/1‑150 (MQ=255)           cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACaa                                                                                                                                       >  1:517530/1‑150 (MQ=255)               cTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTg                                                                                                                                   >  2:1010212/1‑150 (MQ=255)                  aaaCATCTTTCGCTTGTTTTCCTGATTTCCAGCCGGCGG‑TGTGGTCCCCGCAAACCCACCG‑CTGGGGTTGGGCTTTC‑CCAGCCCGGCCAGCAGGCGCAGCAGGTTATTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCt                                                                                                                               <  2:892465/151‑1 (MQ=255)                  aaaCATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCt                                                                                                                               >  2:51839/1‑151 (MQ=255)                     catcatTCGCGTGACTTCCTGAATTTCAGCCAGCGG‑ATTGGCGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCTAGCAGCCGCAGCAGGGTACTTTTGCCACCACCGCCGCGGCCCACCACCGCCACAAACTGAccaacc                                                                                                                            >  2:296256/1‑147 (MQ=255)                          ttCGCGTGTCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGCACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAat                                                                                                                       >  1:183365/1‑151 (MQ=255)                            cgcgTGTCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAatat                                                                                                                     >  1:752156/1‑151 (MQ=255)                                 tttctCCTGAATTTCCCCCCGGGGTTGGGGTCCCCGCAACCCCACCCCCGGGGTTG‑GGCTTTC‑CCAGCCCTCCCGGCGGGCGCAGCGGGTTATTTTCCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTa                                                                                                                  <  2:510201/146‑1 (MQ=255)                                 gtCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                                 >  1:1027753/1‑150 (MQ=255)                                  tCTTCCTGATTTCCAGCCGGGGG‑GGTGGGGCCCGCAACCACACCGCCGGGGGTG‑GGCGTTT‑CCAGCCCTGCCAGCAGGCGCAGCGGGGTACTTTTGCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                      <  2:821438/144‑1 (MQ=255)                                  tCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                               >  1:110985/1‑151 (MQ=255)                                   cTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                              >  2:890276/1‑151 (MQ=255)                                    tcccTGAATTTCACCCGGGGG‑TGGGGTGCCCGATAACACACCCCCTGGGTTT‑GGCTTTC‑CCACCCCGCCCAGCAGGCGCAGCGGGTTTTTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                              <  2:721460/148‑1 (MQ=255)                                          atttccagccgggGG‑TGGGGTGCCCCCAAACCCACCGCCGGGGGTT‑GGCTTTC‑CAAGCCCGGCCAGCGGGCGCAGGGGGGTTTTTTCCCCACCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                       <  1:1113104/146‑1 (MQ=255)                                          atttccagccaggGG‑TGGGGTCCCCGCAACCACTGCCCCTGCGGTT‑GGCTTTC‑CCAGCCCTGCCAGCAGGCGCAGCGGGTTACTTTGCCCACCACCGCGGCGCCCCACCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                       <  1:96618/146‑1 (MQ=255)                                            tttCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                     >  2:803946/1‑151 (MQ=255)                                            ttgcagccgggGG‑TGTGGTCCCCGCTAACACAGCGCGTGGGTTT‑GGCTTTC‑CCAGCCCTCCCAGCGGGCGCAGCGGGGTCTTTTCCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                         <  2:383313/144‑1 (MQ=255)                                               cggccaggGG‑TGGGGTGCCCGCTAACACACCGCCTGCGTTT‑GGCTTTC‑CCAGCCCGCCCAGCGGGCGCAGCGGGTTATTTTTCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAg                                                                                                  <  2:559223/149‑1 (MQ=255)                                                   cagcGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                              >  2:1057462/1‑151 (MQ=255)                                                    ggggg‑TGGGGTCCCCGTAAACCCACCGCCGGGGTTT‑GGCTTTC‑CCAGCCCTCCCAGGGGGCGCGGGAGGTTCTTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                              <  2:935128/147‑1 (MQ=255)                                                      gggtgggggtgCCCCCTACCCCACCGCCTGCGGTG‑GGCTTTC‑CCAGCCCCGCCAGGAGGCCCAGCGGGGTCCTTTGCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                              <  2:138838/142‑1 (MQ=255)                                                      cGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                            >  1:930232/1‑150 (MQ=255)                                                       gg‑TGTGGTCCCCGCAAACACATCGCCTGGGGTT‑GGCGTTC‑CCAGCCCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                               <  2:665292/146‑1 (MQ=255)                                                              gggCCCGCTAACACCTCGCCTGGGTTT‑GGCTTTC‑CCGGCCCTCCCCGCAGGCGCAGCGGGGTCTTTTTGCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                     <  2:399131/148‑1 (MQ=255)                                                                 cccGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                 >  2:430331/1‑151 (MQ=255)                                                                 cccGCAAACCCACCGCCTGTGGTG‑GGTGTTC‑CCAGCCCTCCCAGCAGGCCCGACGGGTTATTTTCCCCACCCCCCCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                                 <  1:304966/151‑1 (MQ=255)                                                                   cGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTa                                                                                >  1:1084684/1‑150 (MQ=255)                                                                     cTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                             >  1:147764/1‑151 (MQ=255)                                                                       aaCACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                           >  2:936911/1‑151 (MQ=255)                                                                        acacATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                          >  1:439329/1‑151 (MQ=255)                                                                              ccccTGCGTTT‑TGCGTTC‑CGATCCCGGCCAGGGGGCGCAGCAGGTTATTTTTCCCCCCCCCTTCACGCCCCCTAACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                      <  2:551219/147‑1 (MQ=255)                                                                                ccTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                                  >  2:161514/1‑151 (MQ=255)                                                                                     ggTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                             >  1:930845/1‑151 (MQ=255)                                                                                        t‑GGCGTTTCCCGGCCCGGCCAGCAGGCGCAGCGGGGTCCTTTTCCCCCCCCCGCGGCGGCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                            <  1:308139/150‑1 (MQ=255)                                                                                            gcttccccAGACTTCCCAGGGGGCGCAACAGGGTATTTTGGCCCCCCCCTCTGGGCCCCACCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                           <  2:1116306/142‑1 (MQ=255)                                                                                              tttcccgccccggCCAGGAGGCGCAGCGGGGTCCTTTCGCCACCCCCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                                        <  1:594327/137‑1 (MQ=255)                                                                                              ttt‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaaca                                                      >  2:172807/1‑150 (MQ=255)                                                                                                       cccgCCAGAGGGCGCAGAGGGGTTCTTTTGCCCCCCCCGGGGGGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                              <  1:197804/147‑1 (MQ=255)                                                                                                            cagcagGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTgg                                        >  1:586613/1‑151 (MQ=255)                                                                                                             agcagGCGCGGCAGGTTTCTTTTCCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                                         <  1:532287/149‑1 (MQ=255)                                                                                                              gcgggCGCAGCAGGGTCTTTTCCCCCCCCCCTCGGCCGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTgg                                        <  1:587158/146‑1 (MQ=255)                                                                                                                 ggCGCAGCGGGTTTTTTTTCCCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                    <  2:942011/150‑1 (MQ=255)                                                                                                                    gcgccggggtttttttCCCCCCCCCGCAGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                 <  1:411300/138‑1 (MQ=255)                                                                                                                     cagcagGTTCCTTTTGCCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                 <  2:545488/149‑1 (MQ=255)                                                                                                                      agcagGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                >  2:471753/1‑149 (MQ=255)                                                                                                                          gggtttttttGCCCCCCCCCCCGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                 <  1:619084/138‑1 (MQ=255)                                                                                                                          gggTTCTTTCGCCCCCACCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                                  <  2:241835/143‑1 (MQ=255)                                                                                                                           ggTTCTTTTCCCCCCCCCCCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAg                                   <  2:564058/141‑1 (MQ=255)                                                                                                                            gTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                        >  2:1116530/1‑151 (MQ=255)                                                                                                                                tttccccccccccGCTCCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                      <  1:138880/138‑1 (MQ=255)                                                                                                                                            cGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc        >  1:616864/1‑151 (MQ=255)                                                                                                                                             gCTGGGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc        <  1:417515/150‑1 (MQ=255)                                                                                                                                              cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt      >  1:289992/1‑151 (MQ=255)                                                                                                                                                  ggCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaa   >  1:410060/1‑150 (MQ=255)                                                                                                                                                    cccACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCaaa  >  1:924444/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                  AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGG‑TGTGGTGCCCGCTAACACATCGCCTGCGGTT‑GGCGTTT‑CCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTGCCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAA  >  NC_000913/993770‑994062

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