Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	805,082	G→A	5.0%	W527* (TGG→TAG)	ybhJ →	aconitase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	805,082	0	G	A	5.0%	251.0 / 2.1	79	W527* (TGG→TAG)	ybhJ	aconitase family protein
Reads supporting (aligned to +/- strand):  ref base G (39/36);  new base A (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.30e-01

TATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTAC  >  NC_000913/804934‑805226                                                                                                                                                     |                                                                                                                                                 tATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTg                                                                                                                                                  >  2:1187413/1‑149 (MQ=255) tATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTg                                                                                                                                                  <  1:1366338/149‑1 (MQ=255)  aTTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTgg                                                                                                                                                 >  2:953758/1‑149 (MQ=255)       aCCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGa                                                                                                                                            >  2:1564688/1‑149 (MQ=255)        ccTCTGCCAGCGAACTTGATTGCTGGGACAACGTGTCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGaa                                                                                                                                           >  2:1001447/1‑149 (MQ=255)         ctctGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAAt                                                                                                                                          >  2:313400/1‑149 (MQ=255)                aGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTgc                                                                                                                                   >  1:76389/1‑149 (MQ=255)                aGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTgc                                                                                                                                   >  2:410225/1‑149 (MQ=255)                  cGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTgcgc                                                                                                                                 <  2:1010854/149‑1 (MQ=255)                   gAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGct                                                                                                                                <  2:388416/149‑1 (MQ=255)                     aCTTGATTTCAGGGACACCGTGCCGGAGTAGGCCTTCGATGTAACGACGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGctga                                                                                                                              <  2:598214/149‑1 (MQ=255)                       ttGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGctgact                                                                                                                            <  1:1187413/149‑1 (MQ=255)                         gATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAAGTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTAGCCGGAATTGGGTGCGctgactga                                                                                                                          >  1:19202/1‑149 (MQ=255)                         gATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAAGTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTAGCCGGAATTGGGTGCGctgactga                                                                                                                          >  1:13546/1‑149 (MQ=255)                          aTTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGTGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCtgactgac                                                                                                                         <  1:261567/149‑1 (MQ=255)                          aTTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCtgactgac                                                                                                                         <  1:876695/149‑1 (MQ=255)                           ttGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACa                                                                                                                        <  1:925371/149‑1 (MQ=255)                            tGCTGGGACAACGTGCCGGAGTACTCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACaa                                                                                                                       <  2:696520/149‑1 (MQ=255)                               tGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAatat                                                                                                                    >  2:1446467/1‑149 (MQ=255)                                 ggACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCg                                                                                                                  >  2:1413318/1‑149 (MQ=255)                                      aCGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTg                                                                                                             >  2:1459875/1‑149 (MQ=255)                                      aCGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTg                                                                                                             >  2:1465770/1‑149 (MQ=255)                                         tGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGaaa                                                                                                          <  2:76389/149‑1 (MQ=255)                                           ccGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGACGAACATAAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAgt                                                                                                        <  2:219503/149‑1 (MQ=255)                                            cGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAgtg                                                                                                       <  1:392395/149‑1 (MQ=255)                                               aGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGc                                                                                                    <  2:215543/149‑1 (MQ=255)                                                           aTGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTc                                                                                        >  2:1154443/1‑149 (MQ=255)                                                               aaCGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCTCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTcgacg                                                                                    >  1:220869/1‑149 (MQ=255)                                                                aCGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTcgacga                                                                                   >  2:1363661/1‑149 (MQ=255)                                                                  gCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAg                                                                                 >  1:997539/1‑149 (MQ=255)                                                                   ccGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGAGTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGt                                                                                <  1:867107/149‑1 (MQ=255)                                                                     gTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTgc                                                                              >  1:1359110/1‑148 (MQ=255)                                                                     gTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGa                                                                              <  1:953758/149‑1 (MQ=255)                                                                       ataAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGacc                                                                            >  1:1115251/1‑149 (MQ=255)                                                                          aaaaCCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGAAAATATAGTACTTAAAGTGTTCTCGAAGATCCTCGACGAAGTGaccacc                                                                         >  2:100266/1‑149 (MQ=255)                                                                           aaaCCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGaccacca                                                                        >  1:360331/1‑149 (MQ=255)                                                                                gtgtTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGTTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACcgac                                                                   <  2:1455189/149‑1 (MQ=255)                                                                                gtgtTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACcgac                                                                   >  1:1517528/1‑149 (MQ=255)                                                                                       tCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGa                                                            >  1:559222/1‑149 (MQ=255)                                                                                        cAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGAt                                                           >  1:1319965/1‑149 (MQ=255)                                                                                               ttGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGAttccttcc                                                    <  2:232520/149‑1 (MQ=255)                                                                                                tgtgAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCg                                                   <  1:1111485/149‑1 (MQ=255)                                                                                                tgtgAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCg                                                   <  2:559222/149‑1 (MQ=255)                                                                                                tgtgAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCg                                                   <  1:452719/149‑1 (MQ=255)                                                                                                tgtgAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCg                                                   <  1:313400/149‑1 (MQ=255)                                                                                                  tgAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGt                                                 >  1:1491382/1‑149 (MQ=255)                                                                                                      aGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGaaa                                             <  1:1459875/149‑1 (MQ=255)                                                                                                      aGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGaaa                                             <  1:1465770/149‑1 (MQ=255)                                                                                                      aGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGaaa                                             <  1:1564688/149‑1 (MQ=255)                                                                                                       ggggCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGCCTGACAATATCGTCCTGAAAGAGAGCTCGAAGATCATCGACGAAGTGAACACCAAAGACGAACTGATTACTTCCGGTGAAAc                                            >  1:211964/1‑149 (MQ=255)                                                                                                        gggCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAAcc                                           >  1:960953/1‑149 (MQ=255)                                                                                                          gCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACctc                                         <  1:1413318/149‑1 (MQ=255)                                                                                                           cAAGTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTAGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACctct                                        <  2:19202/149‑1 (MQ=255)                                                                                                           cAAGTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTAGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACctct                                        <  2:13546/149‑1 (MQ=255)                                                                                                           cAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACctct                                        >  1:1348211/1‑149 (MQ=255)                                                                                                           cAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACctct                                        <  1:1446467/149‑1 (MQ=255)                                                                                                             aCTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGCCTGACAATAGCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCtcttc                                      <  1:782204/149‑1 (MQ=255)                                                                                                             aCTCAGCAAACGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCtcttc                                      <  2:133161/149‑1 (MQ=255)                                                                                                                    aaCCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGt                               >  1:543039/1‑149 (MQ=255)                                                                                                                      ccGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTc                             >  1:1535576/1‑149 (MQ=255)                                                                                                                       cGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCa                            <  2:1115251/149‑1 (MQ=255)                                                                                                                       cGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCa                            >  1:672529/1‑149 (MQ=255)                                                                                                                        gCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCaa                           >  2:1555872/1‑149 (MQ=255)                                                                                                                        gCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCaa                           >  1:560353/1‑149 (MQ=255)                                                                                                                        gCTGATTTACGGGCCGAACATTAAAGACTGGCAGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCaa                           <  1:464521/149‑1 (MQ=255)                                                                                                                          tGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAAt                         <  1:1456216/149‑1 (MQ=255)                                                                                                                              ttACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGa                     <  1:1044418/149‑1 (MQ=255)                                                                                                                               tACGGGCCGAACATTTAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGGCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGAt                    <  1:276591/149‑1 (MQ=255)                                                                                                                                aCGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCATCCGGTGAAACCTCTTCTTATCGTTCAAATCCGAtt                   <  1:1001447/149‑1 (MQ=255)                                                                                                                                   ggCCGTACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGAGAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCAGCGACGAACTGATTCCTTACGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGt                <  2:220869/149‑1 (MQ=255)                                                                                                                                    gCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGac                                                                             >  1:997047/1‑87 (MQ=255)                                                                                                                                    gCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGac                                                                             <  2:997047/87‑1 (MQ=255)                                                                                                                                    gCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTc               >  1:851086/1‑149 (MQ=255)                                                                                                                                     ccGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCt              >  1:1366350/1‑149 (MQ=255)                                                                                                                                       gAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTgg            >  2:1164028/1‑149 (MQ=255)                                                                                                                                       gAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTCTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGAGTGGTCTgg            >  2:434604/1‑149 (MQ=255)                                                                                                                                         aCATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTGCAAATCCGATTGGTCTGGCg          >  1:817844/1‑149 (MQ=255)                                                                                                                                            ttAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCGTCTTATCGTTCAAATCCGATTGGTCTGGCGGAg       >  1:1235951/1‑149 (MQ=255)                                                                                                                                             tAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGt      >  2:1428880/1‑149 (MQ=255)                                                                                                                                              aaaGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGtt     >  1:60145/1‑149 (MQ=255)                                                                                                                                              aaaGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTAGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTACCGGTGAAACCACTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGtt     >  2:869443/1‑149 (MQ=255)                                                                                                                                                 gACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTAc  <  1:472610/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTCTTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTAC  >  NC_000913/804934‑805226

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑