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

GGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑GTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTGGGAGTTTTTCTGGCCGTTTATCGCAGGGGCAAAACTGGTGA  >  NC_000913/615669‑616152
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
ggggACAGCGTGGCGGTGGCTCTACCGCGCTCGGTCTGTTTGACCCTGGCACTTCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCGGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGcgcc                                                                                                                                                                                                                                           <  1:176084/251‑1 (MQ=255)
ggggACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGcgcc                                                                                                                                                                                                                                           <  1:59501/251‑1 (MQ=255)
                tGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTcaca                                                                                                                                                                                                                            >  1:282487/1‑250 (MQ=255)
                    aCTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACACGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGc                                                                                                                                                                                                                       >  1:319729/1‑251 (MQ=255)
                    aCTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGc                                                                                                                                                                                                                       >  2:166512/1‑251 (MQ=255)
                     cTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCATGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACCCCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGc                                                                                                                                                                                                                       >  1:13022/1‑250 (MQ=255)
                          gcgcTCGGTCTTTTTGAACCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATcac                                                                                                                                                                                                                  >  1:356592/1‑250 (MQ=255)
                           cgcTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCacac                                                                                                                                                                                                                <  1:442703/251‑1 (MQ=255)
                             cTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACgg                                                                                                                                                                                                              >  2:363997/1‑251 (MQ=255)
                               cGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTCCCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGc                                                                                                                                                                                                             <  2:348925/250‑1 (MQ=255)
                                     ttttGACCCTGGCATGCCAGGCGATAGGTGAAGCTGGAGCGGCATGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGGTGCTGGAAGATGCGCGTCCGTCGCTGTTAATTTCCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTatat                                                                                                                                                                                                       <  2:296690/250‑1 (MQ=255)
                                               ggCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACct                                                                                                                                                                                            <  1:271853/251‑1 (MQ=255)
                                                gCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACctc                                                                                                                                                                                           <  1:277277/251‑1 (MQ=255)
                                                  aCTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGcgccg                                                                                                                                                                                                                                          <  1:105346/202‑1 (MQ=255)
                                                  aCTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGcgccg                                                                                                                                                                                                                                          >  2:105346/1‑202 (MQ=255)
                                                  aCTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTg                                                                                                                                                                                         >  1:155015/1‑251 (MQ=255)
                                                    tCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCt                                                                                                                                                                                                                                        <  2:373801/202‑1 (MQ=255)
                                                    tCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCt                                                                                                                                                                                                                                        >  1:373801/1‑202 (MQ=255)
                                                     ccATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATGCCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTGCCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGc                                                                                                                                                                                       <  2:282487/250‑1 (MQ=255)
                                                     ccATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTgg                                                                                                                                                                                        <  1:373015/249‑1 (MQ=255)
                                                     ccATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTgg                                                                                                                                                                                        <  1:281476/249‑1 (MQ=255)
                                                     ccATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCt                                                                                                                                                                                      <  1:363997/251‑1 (MQ=255)
                                                       aTGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTcc                                                                                                                                                                                    <  1:166512/251‑1 (MQ=255)
                                                                  gAAGCTGGAGCGGCCTGGCTACCGCTGGATGCCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGcc                                                                                                                                                                         <  2:319729/251‑1 (MQ=255)
                                                                        ggAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGaaa                                                                                                                                                                     <  1:246833/249‑1 (MQ=255)
                                                                         gAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAggg                                                                                                                                                                  <  1:306062/251‑1 (MQ=255)
                                                                          aGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAgggg                                                                                                                                                                 <  2:339453/251‑1 (MQ=255)
                                                                              gCCTGGGTACCGCTGGATACCGGCTAGCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGGCGCTGTTAAGTACCACCGACGATCAACGGCCGCGCTTTAGCGATGTTCCCAATTCAACAAGCCTTTGCAATAACGCCCCGCTTACACCGCAGGGCA‑CTGCGCCGCTGCAACGTTCACAACCGCATCACACGGCTTATGTCATCTTTACCTCTGGCTCCCCCGGCAGGCCGACAGGGGTAAt                                                                                                                                                             <  2:85330/251‑1 (MQ=255)
                                                                                  ggCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTc                                                                                                                                                         >  1:87302/1‑251 (MQ=255)
                                                                                      aCCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGc                                                                                                                                                     <  1:95780/251‑1 (MQ=255)
                                                                                             gatgccGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACGGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGc                                                                                                                                              <  1:89913/247‑1 (MQ=255)
                                                                                                aCCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTAt                                                                                                                                           <  1:62000/251‑1 (MQ=255)
                                                                                                         ccGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAAcc                                                                                                                                   <  1:456997/250‑1 (MQ=255)
                                                                                                               gATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCt                                                                                                                            >  1:339174/1‑251 (MQ=255)
                                                                                                                         aaaTGAGGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCaa                                                                                                                  <  1:433038/251‑1 (MQ=255)
                                                                                                                          aaTGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCaa                                                                                                                  >  1:447558/1‑250 (MQ=255)
                                                                                                                                       gATGCGCGTCCGTCTCTGTTAATTACCCCCGCCGATAAACTGCCGCGCTTTAGAGATGTTCCGAATTAAACAACCCTTTCCTATAACGCCCCGTTTACACCGCAGGGCA‑TGGCGCCGCTGCACTTTTCACACCCGCATCCCAGGGTTGATATCATCTTTACCCCGGGCTCCCCCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCt                                                                                                    <  2:386933/251‑1 (MQ=255)
                                                                                                                                              gtccgtcGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGAGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAACGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAACTCATTATCCGCTTACAGGc                                                                                             >  1:290337/1‑251 (MQ=255)
                                                                                                                                                     gCTGTTAATTACCAGCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGACCCGCTTACACCGCAGGGCA‑ATGCGCCGCCGCAACTTTCACAACCGAATAACACGGCTTATATCATCTTTACATCTGGCTCCACCGGCAGGCCGACAGGGGTAATGGTCGGGCAGCAGGCTATCGACAACAGCCTGCTATGGATGCAAAATCATTATCCGCTTACAGGCGACGATg                                                                                      >  1:296227/1‑251 (MQ=255)
                                                                                                                                                     gCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCTGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCCCCGGCAGGCCGAAAGGGGTAAGGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATg                                                                                      <  2:155719/251‑1 (MQ=255)
                                                                                                                                                       tGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTc                                                                                    >  1:163290/1‑251 (MQ=255)
                                                                                                                                                       tGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTACCACTTTAACAAGCCTTTGCTATAACGCACCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTAACTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTc                                                                                    >  1:423823/1‑251 (MQ=255)
                                                                                                                                                        gTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATTTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGt                                                                                     <  1:276986/249‑1 (MQ=255)
                                                                                                                                                         ttAATTACCACCGACGATCAACTCCCGCTCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACCCCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCCCCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCg                                                                                   <  1:53605/250‑1 (MQ=255)
                                                                                                                                                                      ccgaccAACTGCCGCGCTTTATCGATGTTCCCAATTTACCAGGCCTTTTCTATAACTCCCCGCTTACACCGCAGGGAA‑ATGCGCCGCTGCAACTTTCACACCCGCATCAAGCGGCTTATATCATCTTTACCTCGGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACg                                                                     <  2:69627/246‑1 (MQ=255)
                                                                                                                                                                      acgaTCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAAcc                                                                     >  2:356705/1‑250 (MQ=255)
                                                                                                                                                                              cTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGc                                                               >  2:394490/1‑249 (MQ=255)
                                                                                                                                                                               tGCCGCGCTCTCGCGCTGTTCCCCATCAAACAAGCCTTTGCTATCACCCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACCCGGCTTATATCATCTTTACCTCTGGCCCCACCGTCAGGCCGAAAGGGGTAATGGTCGGGCCGACGGCTACAGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAAATGACGTTGCCCAAAAAACGCCGTGCAGt                                                            >  2:443681/1‑251 (MQ=255)
                                                                                                                                                                                        ttAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACCCCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTCCCTCTGGCTCCCCCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGt                                                    <  2:87302/250‑1 (MQ=255)
                                                                                                                                                                                        ttAGCGATGTCCCCAATTTAACAAGCCTTTCCTATAACGCCCCGCTTCCCCCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTCCCTCTGGCTCCACCGGCAGCCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGt                                                    <  2:339174/250‑1 (MQ=255)
                                                                                                                                                                                          aGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGnctc                                                 >  2:408414/1‑251 (MQ=255)
                                                                                                                                                                                               tttccccaatagaaCAAGTCTTTGCTATACCGCCCCGCTTACACCGCAGG‑CACATGCGACGCTGCAACTTTCACAACCGCCTCACCCGCATTATATTATCTTTCCCCTTGGCTCTACCGGGCGGCCTAAAGGGGTAATGGTCGGGCAGACGGTTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGgtg                                             <  2:453953/238‑1 (MQ=255)
                                                                                                                                                                                               tGTTCCCAATTTAACAAGACTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGgtgg                                            >  1:351420/1‑250 (MQ=255)
                                                                                                                                                                                                gTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTCCACAACCGCATCACACGGCTTATATCATCTTTCCCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGgtgt                                            <  1:132450/250‑1 (MQ=255)
                                                                                                                                                                                                  tCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGCTGTATCGGTGTggg                                         >  1:213344/1‑251 (MQ=255)
                                                                                                                                                                                                  tCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTggg                                         >  1:85984/1‑251 (MQ=255)
                                                                                                                                                                                                                     tGCTATAACGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTGGGAGTTTTTCTGGCCGTTTAt                      >  2:378850/1‑251 (MQ=255)
                                                                                                                                                                                                                           aaCGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGtc                                                   >  2:234754/1‑216 (MQ=255)
                                                                                                                                                                                                                           aaCGCCCCGCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGtc                                                   <  1:234754/216‑1 (MQ=255)
                                                                                                                                                                                                                                   gCTTACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGATCCACCGGCAGGCCGCAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTGGGAGTTTTTCTGGCCGTTTATCGCAGGGGCAAAAc        >  1:58552/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ttACACCGCAGGGCA‑ATGCGCCGCTGCAACTTTCACAACCGCATCACAAGGCTTATATCATCTTTACCTCTGGCTCCACAGGCACGCCGAAAGGGGTAATGGTCGGGCAGATGGCTATCGTCAACCGCCTGCTTTGGCTGCAAAATCATTATCCGCTTACAGGCCAAGATGTCGTTGTCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTGGGAGCTTTTCTGGCCGTTTATCGCAGGGGCAACACTg      >  2:359555/1‑251 (MQ=255)
                                                                                                                                                                                                                                         cccGCAGGGCA‑TTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCGATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTGGGAGTTTTTCTGGCCGTTTATCGCAGGGGCAAAACTGGTGa  <  1:53818/250‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               
GGGGACAGCGTGGCGGTGGCACTACCGCGCTCGGTCTTTTTGACCCTGGCACTCCATGCGATAGTTGAAGCTGGAGCGGCCTGGCTACCGCTGGATACCGGCTATCCGGACGATCGCCTGAAAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCA‑GTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCCCAAAAAACGCCGTGCAGTTTTGATGTCTCGGTGTGGGAGTTTTTCTGGCCGTTTATCGCAGGGGCAAAACTGGTGA  >  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: ‑