Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 615,913 G→A S586N (AGT→AAT)  entF → apo‑serine activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913615,9130GA96.7% 203.5 / NA 62S586N (AGT→AAT) entFapo‑serine activating enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/32);  minor base T (0/2);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTGGGAGTTTTTCTGGCCGTTTATCGCAGGGGCAAAACTGGTGA  >  NC_000913/615669‑616152
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
ggggACAGCGTGGCGGTGGCTCTACCGCGCTCGGTCTGTTTGACCCTGGCACTTCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCGGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGcgcc                                                                                                                                                                                                                                           <  1:176084/251‑1 (MQ=255)
ggggACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGcgcc                                                                                                                                                                                                                                           <  1:59501/251‑1 (MQ=255)
                tGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTcaca                                                                                                                                                                                                                            >  1:282487/1‑250 (MQ=255)
                    aCTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACACGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGc                                                                                                                                                                                                                       >  1:319729/1‑251 (MQ=255)
                    aCTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGc                                                                                                                                                                                                                       >  2:166512/1‑251 (MQ=255)
                     cTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCATGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACCCCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGc                                                                                                                                                                                                                       >  1:13022/1‑250 (MQ=255)
                          gcgcTCGGTCTTTTTGAACCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATcac                                                                                                                                                                                                                  >  1:356592/1‑250 (MQ=255)
                           cgcTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCacac                                                                                                                                                                                                                <  1:442703/251‑1 (MQ=255)
                             cTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACgg                                                                                                                                                                                                              >  2:363997/1‑251 (MQ=255)
                               cGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTCCCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGc                                                                                                                                                                                                             <  2:348925/250‑1 (MQ=255)
                                     ttttGACCCTGGCATGCCAGGCGATAGGTGAAGCTGGAGCGGCATGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGGTGCTGGAAGATGCGCGTCCGTCGCTGTTAATTTCCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTatat                                                                                                                                                                                                       <  2:296690/250‑1 (MQ=255)
                                               ggCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACct                                                                                                                                                                                            <  1:271853/251‑1 (MQ=255)
                                                gCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACctc                                                                                                                                                                                           <  1:277277/251‑1 (MQ=255)
                                                  aCTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGcgccg                                                                                                                                                                                                                                          <  1:105346/202‑1 (MQ=255)
                                                  aCTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGcgccg                                                                                                                                                                                                                                          >  2:105346/1‑202 (MQ=255)
                                                  aCTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTg                                                                                                                                                                                         >  1:155015/1‑251 (MQ=255)
                                                    tCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCt                                                                                                                                                                                                                                        <  2:373801/202‑1 (MQ=255)
                                                    tCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCt                                                                                                                                                                                                                                        >  1:373801/1‑202 (MQ=255)
                                                     ccATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATGCCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTGCCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGc                                                                                                                                                                                       <  2:282487/250‑1 (MQ=255)
                                                     ccATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTgg                                                                                                                                                                                        <  1:373015/249‑1 (MQ=255)
                                                     ccATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTgg                                                                                                                                                                                        <  1:281476/249‑1 (MQ=255)
                                                     ccATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCt                                                                                                                                                                                      <  1:363997/251‑1 (MQ=255)
                                                       aTGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTcc                                                                                                                                                                                    <  1:166512/251‑1 (MQ=255)
                                                                  gAAGCTGGAGCGGCCTGGCTACCGCTGGATGCCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGcc                                                                                                                                                                         <  2:319729/251‑1 (MQ=255)
                                                                        ggAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGaaa                                                                                                                                                                     <  1:246833/249‑1 (MQ=255)
                                                                         gAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAggg                                                                                                                                                                  <  1:306062/251‑1 (MQ=255)
                                                                          aGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAgggg                                                                                                                                                                 <  2:339453/251‑1 (MQ=255)
                                                                              gCCTGGGTACCGCTGGATACCGGCTAGCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGGCGCTGTTAAGTACCACCGACGATCAACGGCCGCGCTTTAGCGATGTTCCCAATTCAACAAGCCTTTGCAATAACGCCCCGCTTACACCGCAGGGCACTGCGCCGCTGCAACGTTCACAACCGCATCACACGGCTTATGTCATCTTTACCTCTGGCTCCCCCGGCAGGCCGACAGGGGTAAt                                                                                                                                                             <  2:85330/251‑1 (MQ=255)
                                                                                  ggCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTc                                                                                                                                                         >  1:87302/1‑251 (MQ=255)
                                                                                      aCCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGc                                                                                                                                                     <  1:95780/251‑1 (MQ=255)
                                                                                             gatgccGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACGGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGc                                                                                                                                              <  1:89913/247‑1 (MQ=255)
                                                                                                aCCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTAt                                                                                                                                           <  1:62000/251‑1 (MQ=255)
                                                                                                         ccGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAAcc                                                                                                                                   <  1:456997/250‑1 (MQ=255)
                                                                                                               gATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCt                                                                                                                            >  1:339174/1‑251 (MQ=255)
                                                                                                                         aaaTGAGGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCaa                                                                                                                  <  1:433038/251‑1 (MQ=255)
                                                                                                                          aaTGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCaa                                                                                                                  >  1:447558/1‑250 (MQ=255)
                                                                                                                                       gATGCGCGTCCGTCTCTGTTAATTACCCCCGCCGATAAACTGCCGCGCTTTAGAGATGTTCCGAATTAAACAACCCTTTCCTATAACGCCCCGTTTACACCGCAGGGCATGGCGCCGCTGCACTTTTCACACCCGCATCCCAGGGTTGATATCATCTTTACCCCGGGCTCCCCCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCt                                                                                                    <  2:386933/251‑1 (MQ=255)
                                                                                                                                              gtccgtcGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGAGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAACGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAACTCATTATCCGCTTACAGGc                                                                                             >  1:290337/1‑251 (MQ=255)
                                                                                                                                                     gCTGTTAATTACCAGCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGACCCGCTTACACCGCAGGGCAATGCGCCGCCGCAACTTTCACAACCGAATAACACGGCTTATATCATCTTTACATCTGGCTCCACCGGCAGGCCGACAGGGGTAATGGTCGGGCAGCAGGCTATCGACAACAGCCTGCTATGGATGCAAAATCATTATCCGCTTACAGGCGACGATg                                                                                      >  1:296227/1‑251 (MQ=255)
                                                                                                                                                     gCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCTGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCCCCGGCAGGCCGAAAGGGGTAAGGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATg                                                                                      <  2:155719/251‑1 (MQ=255)
                                                                                                                                                       tGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTc                                                                                    >  1:163290/1‑251 (MQ=255)
                                                                                                                                                       tGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTACCACTTTAACAAGCCTTTGCTATAACGCACCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTAACTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTc                                                                                    >  1:423823/1‑251 (MQ=255)
                                                                                                                                                        gTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATTTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGt                                                                                     <  1:276986/249‑1 (MQ=255)
                                                                                                                                                         ttAATTACCACCGACGATCAACTCCCGCTCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACCCCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCCCCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCg                                                                                   <  1:53605/250‑1 (MQ=255)
                                                                                                                                                                      ccgaccAACTGCCGCGCTTTATCGATGTTCCCAATTTACCAGGCCTTTTCTATAACTCCCCGCTTACACCGCAGGGAAATGCGCCGCTGCAACTTTCACACCCGCATCAAGCGGCTTATATCATCTTTACCTCGGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACg                                                                     <  2:69627/246‑1 (MQ=255)
                                                                                                                                                                      acgaTCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAAcc                                                                     >  2:356705/1‑250 (MQ=255)
                                                                                                                                                                              cTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGc                                                               >  2:394490/1‑249 (MQ=255)
                                                                                                                                                                               tGCCGCGCTCTCGCGCTGTTCCCCATCAAACAAGCCTTTGCTATCACCCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACCCGGCTTATATCATCTTTACCTCTGGCCCCACCGTCAGGCCGAAAGGGGTAATGGTCGGGCCGACGGCTACAGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAAATGACGTTGCCCAAAAAACGCCGTGCAGt                                                            >  2:443681/1‑251 (MQ=255)
                                                                                                                                                                                        ttAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACCCCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTCCCTCTGGCTCCCCCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGt                                                    <  2:87302/250‑1 (MQ=255)
                                                                                                                                                                                        ttAGCGATGTCCCCAATTTAACAAGCCTTTCCTATAACGCCCCGCTTCCCCCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTCCCTCTGGCTCCACCGGCAGCCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGt                                                    <  2:339174/250‑1 (MQ=255)
                                                                                                                                                                                          aGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGnctc                                                 >  2:408414/1‑251 (MQ=255)
                                                                                                                                                                                               tttccccaatagaaCAAGTCTTTGCTATACCGCCCCGCTTACACCGCAGGCACATGCGACGCTGCAACTTTCACAACCGCCTCACCCGCATTATATTATCTTTCCCCTTGGCTCTACCGGGCGGCCTAAAGGGGTAATGGTCGGGCAGACGGTTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGgtg                                             <  2:453953/238‑1 (MQ=255)
                                                                                                                                                                                               tGTTCCCAATTTAACAAGACTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGgtgg                                            >  1:351420/1‑250 (MQ=255)
                                                                                                                                                                                                gTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTCCACAACCGCATCACACGGCTTATATCATCTTTCCCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGgtgt                                            <  1:132450/250‑1 (MQ=255)
                                                                                                                                                                                                  tCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGCTGTATCGGTGTggg                                         >  1:213344/1‑251 (MQ=255)
                                                                                                                                                                                                  tCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTggg                                         >  1:85984/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGCTATAACGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTGGGAGTTTTTCTGGCCGTTTAt                      >  2:378850/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaCGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGtc                                                   >  2:234754/1‑216 (MQ=255)
                                                                                                                                                                                                                           aaCGCCCCGCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGtc                                                   <  1:234754/216‑1 (MQ=255)
                                                                                                                                                                                                                                   gCTTACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGATCCACCGGCAGGCCGCAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTGGGAGTTTTTCTGGCCGTTTATCGCAGGGGCAAAAc        >  1:58552/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ttACACCGCAGGGCAATGCGCCGCTGCAACTTTCACAACCGCATCACAAGGCTTATATCATCTTTACCTCTGGCTCCACAGGCACGCCGAAAGGGGTAATGGTCGGGCAGATGGCTATCGTCAACCGCCTGCTTTGGCTGCAAAATCATTATCCGCTTACAGGCCAAGATGTCGTTGTCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTGGGAGCTTTTCTGGCCGTTTATCGCAGGGGCAACACTg      >  2:359555/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cccGCAGGGCATTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCGATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTGGGAGTTTTTCTGGCCGTTTATCGCAGGGGCAAAACTGGTGa  <  1:53818/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
GGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTGGGAGTTTTTCTGGCCGTTTATCGCAGGGGCAAAACTGGTGA  >  NC_000913/615669‑616152

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: