Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,625,170 C→T H92H (CAC→CAT serB → phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,625,1700CT94.3% 150.5 / ‑2.0 53H92H (CAC→CATserBphosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (29/21);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTG  >  NC_000913/4624925‑4625399
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     
gAAGATGTCGCTTTATGGCCGGGTCTGCCTCTTTCATGAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAgc                                                                                                                                                                                                                                  <  2:247783/251‑1 (MQ=255)
 aaGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAgcg                                                                                                                                                                                                                                 >  1:377296/1‑251 (MQ=255)
   gATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAgcgc                                                                                                                                                                                                                                >  1:112013/1‑250 (MQ=255)
      gTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCaac                                                                                                                                                                                                                             <  2:112013/250‑3 (MQ=255)
        ctctTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCCGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCCGGTTCACTCACCGCACGGGCTACGCGTCTGGCCCATGAATCGCAACTgg                                                                                                                                                                                                                          >  1:224205/1‑251 (MQ=255)
            tgaTGGCCGGGTCTGCCTCTTTCAGTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCCGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGTCCCATGAAGCGCAACTGGATg                                                                                                                                                                                                                       <  2:439434/248‑1 (MQ=255)
            gtATGGCCTGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATTCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCTCATGAAGCGCAACTGGAt                                                                                                                                                                                                                        <  1:72363/248‑1 (MQ=255)
              aTGGCCGAGTCTGCCTCTTTCATTCAGTGGTGATGAAGTGCTGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGGGATTGTTGCCGCCTGGTGCGTGGCAGGTTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCCCGGGCTACGCGCCTGGCCCATGAAgccacactggctgtcc                                                                                                                                                                                                                    >  2:189948/1‑237 (MQ=255)
                             tctTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGAGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGGCCCGCTGGGGaaa                                                                                                                                                                                                      >  1:268837/1‑250 (MQ=255)
                                ttCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAAt                                                                                                                                                                                                    <  1:157606/249‑1 (MQ=255)
                                  cATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATccc                                                                                                                                                                                                 <  1:82185/250‑1 (MQ=255)
                                     tAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCAc                                                                                                                                                                                             >  2:441244/1‑251 (MQ=255)
                                              atgaAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACg                                                                                                                                                                                    >  2:166650/1‑251 (MQ=255)
                                                  aGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCgg                                                                                                                                                                                <  1:263466/251‑1 (MQ=255)
                                                   gTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCggg                                                                                                                                                                               >  1:240153/1‑251 (MQ=255)
                                                    tGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGt                                                                                                                                                                              >  1:3722/1‑251 (MQ=255)
                                                         ccACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGCTTGCt                                                                                                                                                                         >  2:96114/1‑251 (MQ=255)
                                                          cACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTg                                                                                                                                                                        <  1:166650/251‑1 (MQ=255)
                                                                  tACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATgg                                                                                                                                                                 >  1:59888/1‑250 (MQ=255)
                                                                     cACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTTATGGatat                                                                                                                                                             <  2:240153/251‑1 (MQ=255)
                                                                           ggTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGattc                                                                                                                                                       <  1:325873/251‑4 (MQ=255)
                                                                               gTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCGGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCAcc                                                                                                                                                   <  1:441244/251‑1 (MQ=255)
                                                                                                   ggTCGTGGGCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTAt                                                                                                                               >  1:352433/1‑251 (MQ=255)
                                                                                                    gTCGTTGGCTTGATAAGCATCGTCTGACACAATACTATAGCATACTGGGTGCGGCTATGGTGATTGTTGCCGCTTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTCCACGCACCGCACGGGCTACGTGCCTGGCCCACGAAGCGCAACTGTGTGTCGCCCCGCTGGGGAAAAGCTCGCACCTACGCACGCCGTGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTAtt                                                                                                                              <  2:7914/251‑1 (MQ=255)
                                                                                                         gggCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGAATGAATGTATtgatga                                                                                                                         <  1:374444/251‑1 (MQ=255)
                                                                                                          ggCTGGATAAACAGCGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGaa                                                                                                                        <  2:3722/251‑1 (MQ=255)
                                                                                                                ataaacagcgtgagaaccAATACAAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGAGCGTGGAAGATTATCAGGTGATTCGTATGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGAAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGATGGGGCTGGAAATGGATTCCACCGCAATGCAGATTGAATGGATTGAGGAAAtgtac                                                                                                                  >  2:319566/17‑247 (MQ=255)
                                                                                                                    acagcGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGCTTCCACCGCCATCCAGATTGCATGTATTGATGAAATTGCCAAAc                                                                                                              >  2:397743/5‑251 (MQ=255)
                                                                                                                    acagcGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGGGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCCCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGACTGTATTGATGAAATTGCCAAAc                                                                                                              >  1:236390/5‑251 (MQ=255)
                                                                                                                         gTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGcc                                                                                                         >  1:346208/1‑251 (MQ=255)
                                                                                                                                    tACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGGTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAAAGGGCgg                                                                                              >  2:446502/1‑250 (MQ=255)
                                                                                                                                       cAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCaact                                                                                                                                                                                                                            <  2:147840/122‑4 (MQ=255)
                                                                                                                                       cAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCaact                                                                                                                                                                                                                            >  1:147840/1‑119 (MQ=255)
                                                                                                                                        aaaacaaaaggggTGAGGCGATGTTGATTGGTGCCGCCTGGTGCGTGGAAGATTATCAGGTTATTCGGCTGGCAGGTTCACTCACCGCCTGGGCTACGCGCCGGGCCCATTAAGCGCAACTGTATGGCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGACGTGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGAtg                                                                                          <  2:262884/241‑1 (MQ=255)
                                                                                                                                           gCAAACTGGGTTCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTCCGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGAtggtg                                                                                       <  2:290144/251‑1 (MQ=255)
                                                                                                                                             aaaCTGAGTGCGGCGATGGTGAATGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCCCCTGCGCACGCCGGGGTTGCTGGTGCTGGACATGGATTCCACAGCCATCAAGCCTGCATGTATTGATGAAATTGCCCAACTGGCCGGAACGGGCGAGAtggtggg                                                                                     >  2:344464/1‑250 (MQ=255)
                                                                                                                                                          cgATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGCTTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGGGATGGTGGCGGAAGTAACCGa                                                                         >  2:254499/1‑250 (MQ=255)
                                                                                                                                                           aaTGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCAGCTTGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGGTTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGGAGTAACCGaa                                                                        >  2:269072/2‑250 (MQ=255)
                                                                                                                                                              ggTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACggg                                                                    >  1:4363/1‑251 (MQ=255)
                                                                                                                                                                 gATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCgg                                                                 >  2:162568/1‑250 (MQ=255)
                                                                                                                                                                  aTTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGAt                                                                >  1:247665/1‑251 (MQ=255)
                                                                                                                                                                  aTTGTTGCCGCCTGGTGCGTGGAAGATTAGCAGGTGTCTCGTCTGGCAGGTTCACTCACCGCACGGTCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGTACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGAt                                                                <  1:433473/251‑1 (MQ=255)
                                                                                                                                                                        gccgccgtgTGCATGGAAGATTATCAGGTGATTTGTCCGGCAGGTTCCCTCTCCGCCCGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGcgcg                                                           <  1:263503/242‑1 (MQ=255)
                                                                                                                                                                               ggggCGTGGAAGATTATCAGGTGATGCGTCGGGCAGGTTCGATCACCGCACGGGCTACGCGCCTGGTCCATGAAGAGCAACTGGATGTCGTCCCGCGGGGGAAAAGCCCGCAACTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGAGGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTc                                                   <  2:250018/248‑1 (MQ=255)
                                                                                                                                                                               ggTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTgca                                                                                                                                    >  2:109289/1‑168 (MQ=255)
                                                                                                                                                                               ggTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGaa                                                                                                                                    <  1:109289/170‑1 (MQ=255)
                                                                                                                                                                                           aTTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGCGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGcc                                       >  2:380378/1‑251 (MQ=255)
                                                                                                                                                                                                      aTTCGTCTGGCAGGTTCACTCACCGCACGGGCTACGCGCCTCGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCCGCCTgcgc                              >  2:391503/1‑249 (MQ=255)
                                                                                                                                                                                                                     tcactcacCGCACGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCt             >  1:390732/1‑251 (MQ=255)
                                                                                                                                                                                                                     accctcacCGCCCGGGGCACGCGCCTGGCCCATGAACCGCAATTGGATCTTCCCCCGTTGGGGAAAACCCCGCCCCGGTCCACGCCGGGTCTGCTGGTGAGTGACATGGATTCCGCCGCCATCCAGATGGAATGGATTGAAGAAATGGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGATCGGGGGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAgacgcttgtcgccgc              <  2:221253/247‑15 (MQ=255)
                                                                                                                                                                                                                                aCGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATtgatga                                                                                                                         >  1:128071/1‑132 (MQ=255)
                                                                                                                                                                                                                                aCGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATtgatga                                                                                                                         <  2:128071/132‑1 (MQ=255)
                                                                                                                                                                                                                                 cGGGCTACGCGCCTGGCCCATGAAGCGCAACTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGACATGGATTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGCGTGGCGACGCTGAAAGGCGCTg  >  1:198040/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     
GAAGATGTCTCTTTATGGCCGGGTCTGCCTCTTTCATTAAGTGGTGATGAAGTGATGCCACTGGATTACCACGCAGGTCGTAGCGGCTGGCTGCTGTATGGTCGTGGGCTGGATAAACAACGTCTGACCCAATACCAGAGCAAACTGGGTGCGGCGATGGTGATTGTTGCCGCCTGGTGCGTGGAAGATTATCAGGTGATTCGTCTGGCAGGTTCACTCACCGCACGGGCTACACGCCTGGCCCACGAAGCGCAGCTGGATGTCGCCCCGCTGGGGAAAATCCCGCACCTGCGCACGCCGGGTTTGCTGGTGATGGATATGGACTCCACCGCCATCCAGATTGAATGTATTGATGAAATTGCCAAACTGGCCGGAACGGGCGAGATGGTGGCGGAAGTAACCGAACGGGCGATGCGCGGCGAACTCGATTTTACCGCCAGCCTGCGCAGCCGTGTGGCGACGCTGAAAGGCGCTG  >  NC_000913/4624925‑4625399

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