Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,530,797 C→T Q285Q (CAG→CAA sgcX ← putative endoglucanase with Zn‑dependent exopeptidase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,530,7970CT98.7% 248.2 / NA 78Q285Q (CAG→CAAsgcXputative endoglucanase with Zn‑dependent exopeptidase domain
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (39/38);  minor base A (0/1);  total (39/39)
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.92e-01

TGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAAC  >  NC_000913/4530560‑4531042
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
tGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCCACCTCCGCCGGAGAATGGGTATAGCGACACGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATCTAAACGGGTTCTGTTATCACGCCGGGCGCCGCTTCTCGCGCCAACGGAATGTTGAGTACAAGAGCCGTTTGTTACACCATcc                                                                                                                                                                                                                                          >  2:394728/1‑251 (MQ=255)
 gAATCAGGCTCAACGGGGAAATGTGCTGCTGACATAACTGCCAGAGAGGTCACTCGGCGGATGCACTCCGTCAAATAACGAAGGCTGGAAAACTCCGCCGGAGAATGGGTATAACGACAAGGAATAGAGGGACTGGCGCCGGGAATACAATCCTGCTCAACCCGGATAAAAACGGGTTATGTTAGCACGCCGGGAGCAACTTCTAGTTGAAACGGAATGTTGTGTTCAAGAGACGTTTGTTcacaaaaccg                                                                                                                                                                                                                                         >  2:532900/1‑242 (MQ=255)
  aaTCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGTTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGt                                                                                                                                                                                                                                        <  1:405411/251‑1 (MQ=255)
   aTCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGGTATAACCGGTTTCTGTTATCACGCCGGGAGACACTTCTCGCTGCACCGGAATGTTGTGTTCAAGTGCCGTTTGTTCCAACATCCGt                                                                                                                                                                                                                                        >  2:482912/1‑250 (MQ=255)
           cAACGGGGAAATTGGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCg                                                                                                                                                                                                                                <  2:476591/250‑1 (MQ=255)
              cGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAACCGTTTGTTCCAACATCCGTATCAGGCgagg                                                                                                                                                                                                                             >  2:85250/1‑250 (MQ=255)
                  gAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATTTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCg                                                                                                                                                                                                                        <  1:161650/251‑1 (MQ=255)
                     aTGTGCTTCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGa                                                                                                                                                                                                                     <  1:55471/251‑1 (MQ=255)
                            gctgACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATcc                                                                                                                                                                                                              <  1:288694/251‑1 (MQ=255)
                               gACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAAt                                                                                                                                                                                                                <  1:559995/246‑1 (MQ=255)
                               gACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAAt                                                                                                                                                                                                                >  2:559995/1‑246 (MQ=255)
                                          cAGAGCGGTCAATAGGCGGATGTAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCCCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTcccc                                                                                                                                                                                                <  1:199616/251‑1 (MQ=255)
                                                  tCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCg                                                                                                                                                                                                                                >  2:83951/1‑211 (MQ=255)
                                                  tCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCg                                                                                                                                                                                                                                <  1:83951/211‑1 (MQ=255)
                                                        ggCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTc                                                                                                                                                                                  >  1:443293/1‑251 (MQ=255)
                                                           ggATGCACTCAGTCAAATCACGCAGGCTGGCCACCTCCGCCGGAGAATGGGTATAGCGAAAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGCTTAATCCGGCCAGCGTTCCCCGCACTTGGTAGTTcagg                                                                                                                                                                               >  2:279929/1‑251 (MQ=255)
                                                            gagtcagtcattgaaaTCACGCAGGCGGGCATCCTCCGCCGTATAATGGGTAGAGCGACAAGTAATATAGAGGCTGGCTCAGTGAATACCTTCCTGCTCAACCTTGATATAACCGGTTTCTGTTATCACGTCTGGCGCCTCTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTTTTCCAACTTCCTTATCAGGCGAGGAGGCGTGATTAATCCGTCCAGCTTTCCCCGACCATTGTAGTTcagtc                                                                                                                                                                              <  2:313941/238‑3 (MQ=255)
                                                                  aTCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGa                                                                                                                                                                        <  2:132454/251‑1 (MQ=255)
                                                                      gTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATccc                                                                                                                                                                    <  1:160835/251‑1 (MQ=255)
                                                                       tCAAATCACGCAAGATAGAAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGCGAGACTGGAGCAGGGAATACCATCCTGCTCAACCAGGCTATAACCGGTTTCTGTTATCACGCAGGGAGCCACTTCTCGCTGCACAGAAATGTTGTATTCACGAGCCGCTTGGTCCAACATCCGTATCCGTCCATGAGGAGTGCTTAATCAGGCAATCGTACCCCTACCACGGTATTGCAGACAGGTGATccc                                                                                                                                                                    >  2:471560/1‑250 (MQ=255)
                                                                        cAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCa                                                                                                                                                                   >  1:331421/1‑250 (MQ=255)
                                                                          aaTCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACg                                                                                                                                                                >  1:547770/1‑251 (MQ=255)
                                                                           atgaCGCAGACTGGCAACCTCCGCCGGAGAATGGGCAAAGCGACCAGGACTAGAGAGACTGGCGCAGGGAATACCATCCTCCTCAACCTGGATATAACCGGTTTCTGTTATCACGACGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCAGTTTGTTCCAACATCCGTATCAGGCGAGGCGGCGGGATTAATCCGGCCAGAGTTCCCCGACCAGGGTAGTTAAGGCCGGTGCTCCCaccac                                                                                                                                                               >  2:539982/4‑247 (MQ=255)
                                                                                gCAGGCTGGCAACCGCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGa                                                                                                                                                          <  1:19162/251‑1 (MQ=255)
                                                                                  aGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGAtt                                                                                                                                                        <  1:405204/251‑1 (MQ=255)
                                                                                       ggCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATc                                                                                                                                                    <  1:389860/250‑1 (MQ=255)
                                                                                             cTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCa                                                                                                                                                                                                                                                <  2:174713/152‑1 (MQ=255)
                                                                                             cTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCa                                                                                                                                                                                                                                                >  1:174713/1‑152 (MQ=255)
                                                                                               ccgccgGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCg                                                                                                                                           >  1:409649/1‑251 (MQ=255)
                                                                                                  cgggAGAATGGGTTGTGCGACAAGGAATTGAGAGACGGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTACAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGaa                                                                                                                                        <  2:443293/249‑1 (MQ=255)
                                                                                                           gggTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCgtg                                                                                                                                <  1:398403/250‑1 (MQ=255)
                                                                                                           gggTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTGCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCgtg                                                                                                                                <  1:352385/250‑1 (MQ=255)
                                                                                                               aaaGCGCCAAGGAATAGAGAGACCGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCCCTGCTCGCTGCACCGGAAGGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGTCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTTTAAAt                                                                                                                           <  2:192665/249‑1 (MQ=255)
                                                                                                                   agACAAGGAAAAGATAGACTGGCGCAGGTAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTAACTCGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAgg                                                                                                                       <  2:370309/250‑1 (MQ=255)
                                                                                                                    gACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCg                                                                                                                                                                                               >  1:220216/1‑178 (MQ=255)
                                                                                                                    gACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCg                                                                                                                                                                                               <  2:220216/178‑1 (MQ=255)
                                                                                                                          gAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGTTGTc                                                                                                                 >  1:477145/1‑250 (MQ=255)
                                                                                                                              agagagACTGGCGCAGGGACTACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGaatac                                                                                                                                     >  2:48037/1‑225 (MQ=255)
                                                                                                                              agagagACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGaataa                                                                                                                                     <  1:48037/226‑1 (MQ=255)
                                                                                                                                agagaCTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCa                                                                                                                                                                                           >  1:360514/1‑170 (MQ=255)
                                                                                                                                agagaCTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTGCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTTTTCAAGAGCCGTTTGTTCCAACAGCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCa                                                                                                                                                                                           <  2:360514/170‑1 (MQ=255)
                                                                                                                                      tGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTg                                                                                                    >  1:365887/1‑251 (MQ=255)
                                                                                                                                            aGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAg                                                                                                                                                                                                                                   >  1:267911/1‑118 (MQ=255)
                                                                                                                                            aGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGGATCAg                                                                                                                                                                                                                                   <  2:267911/118‑1 (MQ=255)
                                                                                                                                               gAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCCACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGGGatctccat                                                                                            >  2:387631/1‑244 (MQ=255)
                                                                                                                                               gAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATa                                                                                           <  1:319758/251‑1 (MQ=255)
                                                                                                                                                aaTACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGCGTGATATCAATAc                                                                                          >  1:270166/1‑251 (MQ=255)
                                                                                                                                                 aTACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTCGTTCAGGCCGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGAGGGAGTGATATCAATAcc                                                                                         >  2:136041/1‑251 (MQ=255)
                                                                                                                                                   aCCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCa                                                                                        >  1:284566/1‑250 (MQ=255)
                                                                                                                                                      aTCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATAcc                                                                                         <  2:223363/246‑1 (MQ=255)
                                                                                                                                                      aTCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATAcc                                                                                         >  1:223363/1‑246 (MQ=255)
                                                                                                                                                          tGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGcc                                                                                  >  1:97357/1‑249 (MQ=255)
                                                                                                                                                           gCTCATCCTGTAGATAACCGGTTTCTGTTATCATGCCGGGGGCCACTACGCGCGGCACCGGAATGTTGTGTTCAAGAGCCGTTTTTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCATGCAGGTGATCCCAACGCCCTTATTAATACTGACCTCGGAATAATCGTGGAAATCAGGGTTGTCACACGATGGAGTGATATCAATACCATTCGCcaga                                                                               <  2:263647/251‑2 (MQ=255)
                                                                                                                                                           gCTCAACCTGGATGTCACCGGTTTCTGTTATCACGCCGGGAGCCCCTTCTCGCCGCACCGTAATGTTGTGTTCAAGAGCCGTTTGTTCCAACACCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATTCCATCGCCCTGATTAATCCTGACCTCGGAAGAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCca                                                                                 <  1:382050/249‑1 (MQ=255)
                                                                                                                                                              cAACCGGGATATAACCGGTTTCTGTTACCACGCCGGGAGCCACTTCGCGCTGCACCGGAATGTGGTTTTCAAGAGCCGTTTGGTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCATTCGCcagg                                                                               <  2:88849/248‑1 (MQ=255)
                                                                                                                                                                   tGGATATAACCGGATTCTGTTATCACGCCGGGAGCCCCGTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCCTGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAAGGGTGTCACCAGATGTAGTGCTATACATACCAATACCCCGGTCAGGttt                                                                       >  2:308120/1‑251 (MQ=255)
                                                                                                                                                                               gTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTa                                                            >  1:169533/1‑250 (MQ=255)
                                                                                                                                                                                     gTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACgg                                                      <  1:279929/250‑1 (MQ=255)
                                                                                                                                                                                      ttATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGGGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAGTAAGGGGTGTCACACGATGGAGAGATATCATAAACAATCGCCAGGTCCGGTTTTCCACGGCGTAATACgggg                                                    >  2:142420/1‑250 (MQ=255)
                                                                                                                                                                                              ccGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGcc                                            <  1:448032/251‑1 (MQ=255)
                                                                                                                                                                                               cGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCcg                                           <  1:85250/251‑1 (MQ=255)
                                                                                                                                                                                                    gCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGTGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTCATACGGGAACAATg                                              >  2:9666/1‑243 (MQ=255)
                                                                                                                                                                                                    gCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATg                                              <  1:9666/243‑1 (MQ=255)
                                                                                                                                                                                                               cTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAAttctt                           >  1:58885/1‑251 (MQ=255)
                                                                                                                                                                                                               cTGCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATACGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGACCTGATTAATCCTGACCTCGGAATAATCGGGTAAATCAGGGGTGTCACACGATGGAGTGATAACAATACCAATCGCCAGGTCAGGGTTTACACGGCGTAATACGGGAAAAATGCAGCGGATAGTAAAttcat                           >  2:263472/1‑249 (MQ=255)
                                                                                                                                                                                                                 gcaccggtatgtgtggtgcAAGAGTCGTATGGTCCAACATCCGTATCAGGCGAGGAGGCGTGATTATACCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGTAATAATCGTGTAAATCAGGGTTGTCACACGATGGAGTGATATCAATACCAATCGACAGGTCAGGTTTTACATGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTcc                         <  2:188185/233‑1 (MQ=255)
                                                                                                                                                                                                                 gCACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCAAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACAATGGAAGTTCAGGCAGGTGATCCCAACGCCCGGATTAATCATGACGTCGGAATAATCGTGTAAATCAGGGATGTTACAAGATGAAGGGATATCAATACAAATCGCCAGGAAAGGTTTTACACGGCGTACTACGGGAACAATGCAGTGGATATTAAATTCTTcc                         >  1:230070/1‑251 (MQ=255)
                                                                                                                                                                                                                  cACCGGAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCt                        >  1:6264/1‑251 (MQ=255)
                                                                                                                                                                                                                  cACCGGAATGTTGTGTTCAAGAGCCGTTTGTGCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCt                        <  1:520217/251‑1 (MQ=255)
                                                                                                                                                                                                                      ggAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTAc                    <  1:538629/251‑1 (MQ=255)
                                                                                                                                                                                                                       gAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGcc                                                                                  <  2:109455/188‑1 (MQ=255)
                                                                                                                                                                                                                       gAATGTTGTGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGcc                                                                                  >  1:109455/1‑188 (MQ=255)
                                                                                                                                                                                                                            ttgtGTTCAAGAGCCGTTTGTTCCAACATCCGTATCAGGAGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGCCCATGGTAGTTCAGGCAGGGGATCCCAACGCCCTGATTAATCCTGACCTCGGCATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGc              >  2:79416/1‑251 (MQ=255)
                                                                                                                                                                                                                            ttgtGATCAAGAGCCGGTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGAACATGGTAGTTCAGGCAGGTGATACCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCTCACGATGGAGTGATATAACTACCAATCGCCAGGTCATGTTTTACACGGCGTAATACGGGAACAATGAATCGGATATTAAATTCTTACTGGACCGatgc              >  2:365930/1‑248 (MQ=255)
                                                                                                                                                                                                                                 ttCAAGAGCCGTTTGGTCCAACATCCGTATCAGGCGCGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGGTGAGGCAGGTGTTACCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGGTAGCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATACCGCGGATAATAAATTCTTCCTGTCACGAAGcacac          >  2:391820/1‑246 (MQ=255)
                                                                                                                                                                                                                                     agagCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGAt      >  1:266309/1‑250 (MQ=255)
                                                                                                                                                                                                                                       agCCGTTTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATaaa   <  1:58577/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gCCGTTTGTGCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAAc  <  1:394728/251‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                     
TGAATCAGGCTCAACGGGGAAATGTGCTGCTGACATACCTGCCAGAGCGGTCAATAGGCGGATGCAATCAGTCAAATCACGCAGGCTGGCAACCTCCGCCGGAGAATGGGTATAGCGACAAGGAATAGAGAGACTGGCGCAGGGAATACCATCCTGCTCAACCTGGATATAACCGGTTTCTGTTATCACGCCGGGAGCCACTTCTCGCTGCACCGGAATGTTGTGTTCAAGAGCCGTCTGTTCCAACATCCGTATCAGGCGAGGAGGCGTGATTAATCCGGCCAGCGTTCCCCGACCATGGTAGTTCAGGCAGGTGATCCCAACGCCCTGATTAATCCTGACCTCGGAATAATCGTGTAAATCAGGGGTGTCACACGATGGAGTGATATCAATACCAATCGCCAGGTCAGGTTTTACACGGCGTAATACGGGAACAATGCCGCGGATATTAAATTCTTCCTGTACCGAAGCCACCAGATAAAC  >  NC_000913/4530560‑4531042

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: