Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,625,179 G→A Q95Q (CAG→CAA serB → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,625,1790GA98.6% 205.8 / NA 72Q95Q (CAG→CAAserBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (35/36);  minor base T (0/1);  total (35/37)
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.05e-01

TTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATT  >  NC_000913/4624936‑4625410
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       
tttATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGGGCGTGGAAGATTATCAGGGGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATg                                                                                                                                                                                                                                  >  2:183725/1‑251 (MQ=255)
tttATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGATGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGCTTCGTCTGGCAGGTCCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATg                                                                                                                                                                                                                                  >  2:234863/1‑251 (MQ=255)
  tATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGGGATTGTTGCCGCCTGGGGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCCCCGCCCGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTc                                                                                                                                                                                                                                >  2:32862/1‑251 (MQ=255)
          ggTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCt                                                                                                                                                                                                                        <  1:98957/251‑1 (MQ=255)
             cTGCCTCTTTCATTAAGGGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGAGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCCTGAAGCGCAACTCGATGTCGCCCCGCTggg                                                                                                                                                                                                                     >  2:448516/1‑251 (MQ=255)
               gCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTAGCAGGTGATTCGTCTGGCAGGTTCACGCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTgg                                                                                                                                                                                                                      <  2:592008/248‑1 (MQ=255)
                ccTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGa                                                                                                                                                                                                                   >  1:246453/1‑250 (MQ=255)
                     ttCATTAAGTGGTGATGAAGTGAGGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATc                                                                                                                                                                                                              <  2:239405/250‑1 (MQ=255)
                             gtggtgATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATACACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCAGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTg                                                                                                                                                                                                     >  2:416062/1‑251 (MQ=255)
                             gtggtgATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGTGGAAAATCCCGCACCTg                                                                                                                                                                                                     <  1:198811/251‑1 (MQ=255)
                              tggtgATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCCTGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCACACCTgg                                                                                                                                                                                                    >  1:428646/1‑250 (MQ=255)
                               ggtgATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTgc                                                                                                                                                                                                    <  1:167061/250‑1 (MQ=255)
                                gtgATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTgcgc                                                                                                                                                                                                  >  2:481940/1‑251 (MQ=255)
                                 tgatgaCGTGATGCCACTGGATTACCACGCAGGCCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAACCAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCCCTGGGGAAAATCCCGCACCTgcgc                                                                                                                                                                                                  >  2:580433/1‑250 (MQ=255)
                                  gatgaAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCAc                                                                                                                                                                                                >  1:136957/1‑251 (MQ=255)
                                         tGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCgggg                                                                                                                                                                                         >  2:193096/1‑250 (MQ=255)
                                          gATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGtt                                                                                                                                                                                        <  2:197836/251‑1 (MQ=255)
                                           aTGCCACTGGATTACCAAGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCCATACCAGCGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCCGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGttt                                                                                                                                                                                       >  2:619497/1‑251 (MQ=255)
                                                    gATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGAt                                                                                                                                                                              <  1:319995/251‑1 (MQ=255)
                                                        aCCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGac                                                                                                                                                                          <  2:150455/251‑2 (MQ=255)
                                                          catgCAGGTCTGAGCGGCTGGCTGCGGTATGTGCGTGGGCTGGATAAACATCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGTTGCGTTGAAGATTATCAGGGGATTCGTCTGGCAGGGTCACTCACCGCACGGGCTACGCGCCTGGCCCGTTAAGCGCATCTGGAGGTCGCCTCGCTGGTGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGacat                                                                                                                                                                        <  1:37631/248‑4 (MQ=255)
                                                              cAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTCTGCTGGTGATGGACATGGat                                                                                                                                                                    >  1:410005/1‑250 (MQ=255)
                                                                    gTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACCTGGATTCCAcc                                                                                                                                                              >  2:204602/1‑251 (MQ=255)
                                                                      aGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGc                                                                                                                                                            <  1:498341/251‑1 (MQ=255)
                                                                        cGGCTGGCTGCGGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCa                                                                                                                                                          <  1:481940/251‑1 (MQ=255)
                                                                             ggctgCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGGCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAg                                                                                                                                                     <  1:314307/251‑1 (MQ=255)
                                                                                tgctgTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGAtt                                                                                                                                                  >  1:275042/1‑251 (MQ=255)
                                                                                tgctgTATGGACGACGGCTGGATAAACAGCGTCTAACCCAATACCAGAGCAAACTGGGTGCGGAGATGGTGCTTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGATCACTCACCGCACGGGATACGAGACTGGCCCATGAAGCGAAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACATGCGCACGCCGGGTTTGCTGGTGATGGACATGGCTTCAACCGCCATCCCGAtt                                                                                                                                                  >  2:97577/1‑251 (MQ=255)
                                                                                 gctgTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTGCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCACCCAGAt                                                                                                                                                   <  2:187735/249‑1 (MQ=255)
                                                                                   tgTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGa                                                                                                                                                >  2:212889/1‑250 (MQ=255)
                                                                                    gTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGaa                                                                                                                                               <  1:438046/250‑1 (MQ=255)
                                                                                     tATGGACGTGGTCTGGATAAACAGCGTCTGTACCAAGTCCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAAt                                                                                                                                              <  2:427274/250‑1 (MQ=255)
                                                                                              gggCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATtgatg                                                                                                                                     >  1:335424/1‑250 (MQ=255)
                                                                                                gCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGaaa                                                                                                                                  <  2:189215/251‑1 (MQ=255)
                                                                                                   ggATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTg                                                                                                                               <  2:497007/251‑1 (MQ=255)
                                                                                                   ggATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATT‑‑ACCGCCATCCAGATTGAACGTATTGATGAAATTg                                                                                                                               <  1:580433/249‑1 (MQ=255)
                                                                                                   ggATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCAGGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACTCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAAt                                                                                                                                 <  2:90236/249‑1 (MQ=255)
                                                                                                    gATAAACAGCGTCTGACCCAATACCAGAGCAAACTTGGTGTGGCGATGGTGATTGTTGCCGCCTGGTGCGTCGAAGACTATCAGGCGCTTCGTCTCGCAGGTTCACTCACCGCACGAGCTACGCGCCTGGCCCCTGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGCTTGCTGGAGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGCAAtgac                                                                                                                              >  1:329398/1‑248 (MQ=255)
                                                                                                       aaaCAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGAGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCCCGGGCTACGCGCCTGGCCCATGAAGCGCGACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCa                                                                                                                            >  1:567620/1‑250 (MQ=255)
                                                                                                           agcGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTg                                                                                                                       <  1:399209/249‑1 (MQ=255)
                                                                                                           agcGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGAGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTAGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATGCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTg                                                                                                                       >  2:174512/3‑251 (MQ=255)
                                                                                                             cGTCTGACCCAATAACAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTACGTGGAAGATTATCAGGCGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGc                                                                                                                     >  2:235296/1‑251 (MQ=255)
                                                                                                              gTCTGACCCAATACCAGAGCAAACGGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCGGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAAGCCCGCACCTGCGCTCGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGcc                                                                                                                    <  2:291020/251‑1 (MQ=255)
                                                                                                                cTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGCTTGAATGTATTGATGAAATTGCCAAACTGGCCgg                                                                                                                  >  1:39861/1‑251 (MQ=255)
                                                                                                                   aCCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAAc                                                                                                               <  1:378152/251‑1 (MQ=255)
                                                                                                                   aCCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAAc                                                                                                               <  1:193096/251‑1 (MQ=255)
                                                                                                                        aTACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCg                                                                                                          >  1:352433/1‑251 (MQ=255)
                                                                                                                                        ggTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGt                                                                                          <  2:39861/251‑1 (MQ=255)
                                                                                                                                             ggcgATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCGCGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCCTCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCg                                                                                     >  2:261785/1‑251 (MQ=255)
                                                                                                                                                gATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAAc                                                                                  <  1:67601/251‑1 (MQ=255)
                                                                                                                                                gATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAAc                                                                                  <  1:488448/251‑1 (MQ=255)
                                                                                                                                                      gATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGTGCGAGATGGTGGCGGAAGTAACCGAACGGGCGa                                                                            <  1:183725/251‑1 (MQ=255)
                                                                                                                                                      gATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCg                                                                             >  2:208473/1‑250 (MQ=255)
                                                                                                                                                          gttgCCGCCTGGTGCGTGGAAGATTCTCAGGTGATTCATATGGCAGGTTCACTGACCGCACGGGCTACGCGCCTGGCCCATGACGCGCAACTGGATGTCACACCGCTGGGGAAAATCCCGCCCATGCGCACGCAGGGTTTGCTGGTGATGGACATGGATCCCACCGCCATCCAGATTAAATGAATTGATGAAATTGCCAAACTGGCCGGAACGGGCGCGCTGGTGGCGGAGGTAACCGCACGGGCGATgcg                                                                        >  2:573654/1‑251 (MQ=255)
                                                                                                                                                             gccgccTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGcgcg                                                                      <  1:17115/250‑1 (MQ=255)
                                                                                                                                                             gccgccTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCgcgg                                                                     <  2:396077/251‑1 (MQ=255)
                                                                                                                                                                  cTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCACACTGGCCGGAACGGGCGAGATGGAGGCGGAAGTAACCGAACGGGCGAGGCGCGGCGAAc                                                                >  2:298683/1‑251 (MQ=255)
                                                                                                                                                                        cgaggAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGAtt                                                          >  2:410468/4‑251 (MQ=255)
                                                                                                                                                                             aaGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTAcc                                                     <  2:546368/251‑1 (MQ=255)
                                                                                                                                                                               gATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGc                                                   >  1:390110/1‑251 (MQ=255)
                                                                                                                                                                                     cAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCAGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCt                                             >  1:54629/1‑251 (MQ=255)
                                                                                                                                                                                      gggttaTTCGTCTGGCAGTTTCACGCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAACCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGTATTCCACCGCCATCCAGATTGAATGTATTGAGGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTg                                            <  2:399398/246‑1 (MQ=255)
                                                                                                                                                                                      aGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCGCGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTg                                            <  1:261785/251‑1 (MQ=255)
                                                                                                                                                                                                          tcactcacCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCt                        >  1:69006/1‑251 (MQ=255)
                                                                                                                                                                                                           cactcacCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACCTGGATTCCCCCGCCATCCAGACTGACTGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAAATCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTg                       >  2:592176/1‑251 (MQ=255)
                                                                                                                                                                                                                 ccGCACGGGCTTCGCGCCTGCCCCATGAGGCGCAACTGGATGTCGCCCCGCTGGGGAAATTCTCTCACCTTCGCACACCGGTTCTGCTGGTGAGGGACATTGAGTCCACCGCCATCCAGATTGAATGTATTGATTAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGTAAGTATCCTAACGGGCGATGCGCGGCGAACTCGATTGTACCGCCAGCCTGCGCAGCCGCTTGGCGACGCTGAAAGgc                 <  1:81099/251‑1 (MQ=255)
                                                                                                                                                                                                                    cACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAACCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCAGCCAGAGTGAATGTATTGATGAAAGTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGAGGCGACGCTGAAAGGCGCt              <  1:448516/251‑1 (MQ=255)
                                                                                                                                                                                                                           tacgcGCCTGGCCCAAGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGCTTGCATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGCAACCGAACGGGCGATGCGCGGCGACCTCCATTTTACCGCCAGCCTGCGCAGCCCCCTGGCGCCGCTGACACGCGCTgcccccc       >  2:240806/5‑245 (MQ=255)
                                                                                                                                                                                                                            acgcGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCac      >  1:487265/4‑250 (MQ=255)
                                                                                                                                                                                                                            acgcGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACAGGGATTCCACCGCCATCCAGATTGAATGGATTGATGAAATGGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGAGGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGcc        <  1:235296/246‑1 (MQ=255)
                                                                                                                                                                                                                              gcGCCGGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCCCCTGCGCCCGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAata    <  2:182630/250‑1 (MQ=255)
                                                                                                                                                                                                                                gCCTGGCCCATGAAGCGCAACGGGATGTCGCCCCGCTGGGGAAAATCCCGCCCCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTGACGCCAATAtt  <  2:136957/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       
TTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTGACGCCAATATT  >  NC_000913/4624936‑4625410

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