Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,164,581	A→G	6.1%	K163E (AAA→GAA)	nagZ →	beta N‑acetyl‑glucosaminidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,164,581	0	A	G	6.1%	199.2 / 2.6	65	K163E (AAA→GAA)	nagZ	beta N‑acetyl‑glucosaminidase
Reads supporting (aligned to +/- strand):  ref base A (34/27);  new base G (2/2);  total (36/29)
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.82e-01

GATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCC  >  NC_000913/1164434‑1164726                                                                                                                                                    |                                                                                                                                                  gATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGac                                                                                                                                                  >  8:50778/1‑148 (MQ=255)  atatCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCTGAATGAAAACGACCGGGaaa                                                                                                                                                 <  2:47338/149‑1 (MQ=255)  atatCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGaaa                                                                                                                                                 >  6:341862/1‑149 (MQ=255)    atCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAaca                                                                                                                                               >  7:118042/1‑149 (MQ=255)     tCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGaaatcg                                                                                                                                              >  3:126140/1‑146 (MQ=255)     tCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAacac                                                                                                                                              >  7:365668/1‑149 (MQ=255)                ccTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGaca                                                                                                                                   <  1:398470/149‑1 (MQ=255)                 cTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGacac                                                                                                                                  <  8:365668/149‑1 (MQ=255)                            gTCGGGCATATCAGCGCGGCGAGTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCCGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAAc                                                                                                                       >  5:266001/1‑149 (MQ=255)                             tCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGCTCCACAAAAAGCCCTGGCACTTGCCAGCCGGTTTATTGATGGTATGAATGAAGCCGGAATGAAAACGACCGGGGAACACTTCCCAGGACACGGTGCAGTAACg                                                                                                                      >  5:40462/1‑149 (MQ=255)                             tCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGGAACACTTCCCAGGACACGGTGCAGTAACg                                                                                                                      >  5:46204/1‑149 (MQ=255)                               gggCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACgg                                                                                                                     >  4:109850/1‑148 (MQ=255)                                ggCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCa                                                                                                                   <  2:308642/149‑1 (MQ=255)                                ggCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCa                                                                                                                   <  1:180044/149‑1 (MQ=255)                                    tatCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACt                                                                                                               <  5:92320/149‑1 (MQ=255)                                      tCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTca                                                                                                             <  3:109850/149‑1 (MQ=255)                                                gATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGaaa                                                                                                   <  5:341862/149‑1 (MQ=255)                                                  ttGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTAGGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGTAGTAACGGCAGACTCACACAAAGAAaca                                                                                                 <  1:234187/149‑1 (MQ=255)                                                    tgCGAGCGTTCTAATCATGCAGTTGGACTAAAAGCCCTGGTTATTGCCAGACGGTTTATTGATGGTATGCATGAAGCCTGAATGAAAACGACCGGGAAACACTTCACAGGACACGGTGCAGTAACGGCAGACTCACACTAAGAAacac                                                                                                <  4:137529/147‑1 (MQ=255)                                                     gCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCTGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGa                                                                                                     <  7:87523/142‑1 (MQ=255)                                                     gCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGAAACGGTGCAGTAACGGCAGACTCACACAAAGa                                                                                                     >  8:87523/1‑142 (MQ=255)                                                                 aTCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGt                                                                                  >  8:298219/1‑149 (MQ=255)                                                                  taaTGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGCAATGAAAACCACCGGGAAACACTTCACAGGACACGGTGCAGTAAGGGCAGACTCACACTAAGAACCACCGTGCGATCCACGTc                                                                                 >  2:400995/3‑149 (MQ=255)                                                                            ccACAAAAAGCCCTTGCAATTGCCAGCCTGTTTATTGTTTGTATGCATGAAGCCGTAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGga                                                                       <  5:269112/149‑1 (MQ=255)                                                                             cacaAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGga                                                                       <  1:104579/148‑1 (MQ=255)                                                                              acaAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCg                                                                         >  1:127094/1‑145 (MQ=255)                                                                              acaAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCg                                                                         >  8:214061/1‑145 (MQ=255)                                                                              acaAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCg                                                                         >  2:110399/1‑145 (MQ=255)                                                                              acaAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCg                                                                         >  5:271475/1‑145 (MQ=255)                                                                               caAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGCCCGGGAAACCCTTCCCAGGACACGGTGCAGTAACGGCAGCCTCCCACAAAGACACACCGTGCGATCCCCGTCCACAAGCGGAGAt                                                                    >  4:86955/1‑149 (MQ=255)                                                                                aaaaaGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCg                                                                         >  8:297403/1‑143 (MQ=255)                                                                                aaaaaGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCg                                                                         <  7:297403/143‑1 (MQ=255)                                                                                  aaaGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTcg                                                                 >  8:40408/1‑149 (MQ=255)                                                                                  aaaGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTcg                                                                 >  7:123373/1‑149 (MQ=255)                                                                                         tGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGTGATCCACGTCCACAAGCGGAGATTCGCGCTa                                                            >  8:87225/1‑147 (MQ=255)                                                                                         tGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTa                                                            >  4:4228/1‑147 (MQ=255)                                                                                           gCAATTGCCAGCCGGTTTATTCATGGTATGCATGAAGCCGTAATGAAAAAGACCGGGAAACAATTCCCAGGACACGGAGCAGTAACGGCAGACTCACACAAAGAAACACCGTGAGATCCACTTCCACAAGCGGAGATTCGCGCTaaa                                                          <  1:400995/147‑1 (MQ=255)                                                                                                  ccAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTc                                                  <  1:110399/148‑1 (MQ=255)                                                                                                  ccAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCg                                                 <  2:127094/149‑1 (MQ=255)                                                                                                       cGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCtt                                            <  6:271475/149‑1 (MQ=255)                                                                                                           ttATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATCCGCGCTAAAGATATGt                                                   <  8:366375/138‑1 (MQ=255)                                                                                                           ttATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATCCGCGCTAAAGATATGt                                                   >  7:366375/1‑138 (MQ=255)                                                                                                              ttGATTGTAGGCAGGACCCCGGAATGAAAACGCCCGGGGAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTc                                      <  6:40462/148‑1 (MQ=255)                                                                                                              ttGATGGTATGCATGAAGCCGGAATGAAATCGACCGGGGAACACTTCGCAGGACAGGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTc                                      <  6:46204/148‑1 (MQ=255)                                                                                                              ttGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGTGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTcc                                     <  7:87225/149‑1 (MQ=255)                                                                                                              ttGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATg                                                    >  8:14163/1‑134 (MQ=255)                                                                                                              ttGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATg                                                    <  7:14163/134‑1 (MQ=255)                                                                                                              ttGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTcc                                     >  7:390836/1‑149 (MQ=255)                                                                                                              ttGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAAAACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGGCCACACGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTcc                                     >  5:11645/1‑149 (MQ=255)                                                                                                                 aTGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCtt                                   >  3:276461/1‑148 (MQ=255)                                                                                                                    gTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCGCACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCt                                    <  4:129486/144‑1 (MQ=255)                                                                                                                    gTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCGCACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCt                                    >  3:129486/1‑144 (MQ=255)                                                                                                                    gTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATc                               <  3:4228/149‑1 (MQ=255)                                                                                                                      aTGCATGAAGCCGGAATGAAAACGACCGGGAAACCCTTCCCAGGCCACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCcg                             >  4:183875/1‑149 (MQ=255)                                                                                                                          aTGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGa                          >  1:407958/1‑148 (MQ=255)                                                                                                                           tGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGaa                         >  5:229901/1‑148 (MQ=255)                                                                                                                            gAAGCCGGAATGAAAACGACAGGGAAACACTTCCCAGGACACGGTGCAGTAACCGCAGACTCACACAAACAACCACCGTCCGATCCACATCCACAAGCGGAGATTCGCGCTAAATATATGTCGGTCTTCAGGTCCTTAATCAGCGaaaa                       >  7:326474/1‑149 (MQ=255)                                                                                                                              aGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGGCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAAt                      >  5:232130/1‑148 (MQ=255)                                                                                                                               gCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATaa                    <  2:407958/149‑1 (MQ=255)                                                                                                                                ccGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATaaa                   <  7:138661/149‑1 (MQ=255)                                                                                                                                  ggAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACt                 >  7:274713/1‑149 (MQ=255)                                                                                                                                        aaaaCGCCCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCCCGTCCCCAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTCATCCGCGAAAAGAAACTCGACGc           >  1:124251/1‑149 (MQ=255)                                                                                                                                          aaCGACCGGGAAACACTTGCCAGGACACGGTGCAGTTACGGCTGACTCACACAAAGTAACTCCGTGCGATCCACGTCCTGAAGCGGTGTTTGGCGCTTATGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCa         <  4:276461/149‑1 (MQ=255)                                                                                                                                          aaCGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGt                                        >  3:126141/1‑118 (MQ=255)                                                                                                                                          aaCGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGt                                        <  4:126141/118‑1 (MQ=255)                                                                                                                                             gACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTa      >  4:335363/1‑149 (MQ=255)                                                                                                                                             gACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTa      <  8:390836/149‑1 (MQ=255)                                                                                                                                               ccGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATg    >  5:341038/1‑149 (MQ=255)                                                                                                                                                 gggAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGCAACACCCTGCGATCCACGTCCACAAGCGGACATTCGCGCTAAAGATAAGTCGGTCTTCAGTTCCATAATCATCGAAAATAAACTCGAAGACATAATGcc  >  8:175689/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GATATCAGCTTTGCGCCTGTGCTGGATGTCGGGCATATCAGCGCGGCGATTGGCGAGCGTTCTTATCATGCCGATCCACAAAAAGCCCTGGCAATTGCCAGCCGGTTTATTGATGGTATGCATGAAGCCGGAATGAAAACGACCGGGAAACACTTCCCAGGACACGGTGCAGTAACGGCAGACTCACACAAAGAAACACCGTGCGATCCACGTCCACAAGCGGAGATTCGCGCTAAAGATATGTCGGTCTTCAGTTCCTTAATCCGCGAAAATAAACTCGACGCCATTATGCC  >  NC_000913/1164434‑1164726

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