Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	805,043	C→G	5.4%	T514S (ACT→AGT)	ybhJ →	aconitase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	805,043	0	C	G	5.4%	250.4 / 4.6	74	T514S (ACT→AGT)	ybhJ	aconitase family protein
Reads supporting (aligned to +/- strand):  ref base C (33/37);  new base G (2/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

GCTCGTTCTATCGCTGCGACTGCGGCAAACGGTGGCTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTC  >  NC_000913/804898‑805190                                                                                                                                                  |                                                                                                                                                    gCTCGTTCTATCGCTGCGACTGCGGCAAACGGTGGCTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCa                                                                                                                                                  <  1:775066/149‑1 (MQ=255)   tCGTTCTATCGCTGCGACTGCGGCAAACGGTGGCTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACCACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCCGc                                                                                                                                                >  2:276591/1‑149 (MQ=255)            cGCTGCGACTGCGGCAAACGGTGGCTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGa                                                                                                                                       <  1:1452708/149‑1 (MQ=255)             gCTGCGACTGCGGCAAACGGTGGCTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGAt                                                                                                                                      <  2:1098103/149‑1 (MQ=255)              cTGCGACTGCGGCAAACGGTGGCTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGAtt                                                                                                                                     >  2:867107/1‑149 (MQ=255)                  gACTGCGGCAAACGGTGGCTAGTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACg                                                                                                                                 <  2:642723/149‑1 (MQ=255)                   ccTGCGGCAAACGGTGGCTCTTAAACCTCTGCCAGCGAACTTGATAGCTGGGACAACGGGCCGGAGTACGCCTTCGATGTAACGCCGCATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACgg                                                                                                                                <  1:89531/148‑1 (MQ=255)                      gCGGCAAACGGTGGCTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGcc                                                                                                                             >  1:388416/1‑149 (MQ=255)                          cAAACGGTGGCTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAAc                                                                                                                         >  2:472610/1‑149 (MQ=255)                            aaCGGTGGCTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACAt                                                                                                                       >  2:1111485/1‑149 (MQ=255)                                gTGGCTATTTAACGTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCTGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTaaa                                                                                                                   <  2:322062/149‑1 (MQ=255)                                  ggCTATTTAACCTCTGCCAGCGAACTTGATTTCTGTGACAACGTGCCTGAGTACGCCTTCGATGTAACTCCGTATAAAAACAGTGTTTATCAGGGCTTTGTGAAAGTGGCAACTCAGCAACCGGTGATTTACGGTCCGAACATTAAAGa                                                                                                                 <  1:240122/149‑1 (MQ=255)                                  ggCTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGa                                                                                                                 <  2:518682/149‑1 (MQ=255)                                   gCTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCATATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGAc                                                                                                                <  2:570319/149‑1 (MQ=255)                                    cTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACt                                                                                                               >  2:1044418/1‑149 (MQ=255)                                    cTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACt                                                                                                               >  1:232520/1‑149 (MQ=255)                                     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                                                                                               >  2:410225/1‑149 (MQ=255)                                                    aGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTgc                                                                                               >  1:76389/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:1465770/1‑149 (MQ=255)                                                                          aCGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTg                                                                         >  2:1459875/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:313400/149‑1 (MQ=255)                                                                                                                                    tgtgAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCg               <  2:559222/149‑1 (MQ=255)                                                                                                                                    tgtgAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCg               <  1:452719/149‑1 (MQ=255)                                                                                                                                    tgtgAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCg               <  1:1111485/149‑1 (MQ=255)                                                                                                                                      tgAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGt             >  1:1491382/1‑149 (MQ=255)                                                                                                                                          aGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGaaa         <  1:1564688/149‑1 (MQ=255)                                                                                                                                          aGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGaaa         <  1:1465770/149‑1 (MQ=255)                                                                                                                                          aGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGaaa         <  1:1459875/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:1446467/149‑1 (MQ=255)                                                                                                                                               cAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACctct    >  1:1348211/1‑149 (MQ=255)                                                                                                                                                 aCTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGCCTGACAATAGCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCtcttc  <  1:782204/149‑1 (MQ=255)                                                                                                                                                 aCTCAGCAAACGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCtcttc  <  2:133161/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GCTCGTTCTATCGCTGCGACTGCGGCAAACGGTGGCTATTTAACCTCTGCCAGCGAACTTGATTGCTGGGACAACGTGCCGGAGTACGCCTTCGATGTAACGCCGTATAAAAACCGTGTTTATCAGGGCTTTGTGAAAGGGGCAACTCAGCAACCGCTGATTTACGGGCCGAACATTAAAGACTGGCCGGAATTGGGTGCGCTGACTGACAATATCGTCCTGAAAGTGTGCTCGAAGATCCTCGACGAAGTGACCACCACCGACGAACTGATTCCTTCCGGTGAAACCTCTTC  >  NC_000913/804898‑805190

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