Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,621,302 T→C D140D (GAT→GAC deoD → purine nucleoside phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,621,3020TC98.6% 227.5 / NA 71D140D (GAT→GACdeoDpurine nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (41/29);  total (41/30)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

ATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACT  >  NC_000913/4621063‑4621542
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                
aTGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGc                                                                                                                                                                                                                                       <  2:579985/251‑1 (MQ=255)
  gggTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCa                                                                                                                                                                                                                                      <  1:498800/250‑1 (MQ=255)
         ggTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGtgt                                                                                                                                                                                                                              >  1:344305/1‑249 (MQ=255)
            aTGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTg                                                                                                                                                                                                                            >  1:508958/1‑250 (MQ=255)
              gggTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGCAAATTATCCGCGTGGGTTCCTGTGGCGCCGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACCGGGTGTTggtg                                                                                                                                                                                                                         >  2:267503/1‑248 (MQ=255)
                     ccGTCCTGCTCCATCTACAACAAAGAACTGATCACCGATTGCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTAGCGGTATGGGTGCCTGCACCGATTCCAAAGTGAACCGTATCCGTGTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGt                                                                                                                                                                                                                  <  2:120955/251‑1 (MQ=255)
                      cGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTg                                                                                                                                                                                                                 >  1:261271/1‑251 (MQ=255)
                         ccTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTgggg                                                                                                                                                                                                              >  2:209549/1‑250 (MQ=255)
                            gCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTCTCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCCCGTAAAACTGCGCGACGTCGTTATCGGTCTGGGTGCCTGCACCGATTCCAAAGTTAACAGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGAGCGTAACGCGGTAGACGAGGCAAAAGCACTGGGTGTTGATGcacgggttgggaac                                                                                                                                                                                                           >  2:444011/1‑238 (MQ=255)
                            gCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGGGTTGATGCTCGCGTGGGTAAc                                                                                                                                                                                                           >  1:422671/1‑251 (MQ=255)
                              tCCATCTACACCAAAGAACTGATCACCGATTTCGGCGAGAAGAACATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAACGACCATGAATTTGCAGCTATAGCTGACTTTGACATGGTGCGTAACGCGGAAGACGCGGCTAAAGCACTGGGTGGTGATGCACGCGTGTGTAACCt                                                                                                                                                                                                         >  2:586317/1‑251 (MQ=255)
                              cccATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCt                                                                                                                                                                                                         <  1:325602/250‑1 (MQ=255)
                                cATCTACACCAAAGAACTGATCACCGATTTCGGCTTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGt                                                                                                                                                                                                       >  1:169865/1‑251 (MQ=255)
                                cATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGt                                                                                                                                                                                                       >  1:385104/1‑251 (MQ=255)
                                cATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGAAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCCCTAGGTGTTGATGATCGCGTGGGGAAccagt                                                                                                                                                                                                       >  2:498748/1‑248 (MQ=255)
                                             gAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGCGTTGATGCTCGCGTGGGTCACCTGTTCTCCGCTGCCCt                                                                                                                                                                                          >  2:645574/1‑251 (MQ=255)
                                               aCTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGt                                                                                                                                                                                        >  1:71038/1‑251 (MQ=255)
                                                   aTCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGGAACCTGTTCTCCGCTGACCTGTTCt                                                                                                                                                                                     >  2:250254/1‑250 (MQ=255)
                                                        cGATTTCGGCTTGAAGAAACTTCACCGAGTGGATTACTGTGGCGCAGTTCGGCCGCACGTAAAACTGCACAACGTCGTTATCGATATGGGTGCCTGCACCGATTCCAAAGTTAAACGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGAGGTAGACGCGGGCAACGCACTGGGGGGTGAAGCTCGCGGGGGGAAACTGTTCTCCACTGACCTGTTCTACtctc                                                                                                                                                                               >  2:448315/1‑251 (MQ=255)
                                                         gATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTcc                                                                                                                                                                              <  1:645574/251‑1 (MQ=255)
                                                           tttCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCgg                                                                                                                                                                            >  1:171848/1‑251 (MQ=255)
                                                              cGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATTCTCTCGTGTGTAACCTGTGCTCCGCTGTCCTGTTCTACTCTCCGGAc                                                                                                                                                                          <  1:250254/250‑1 (MQ=255)
                                                              cGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACg                                                                                                                                                                         >  1:149135/1‑251 (MQ=255)
                                                              cGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGGGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGGGGGTAAACTGTTCTCCGCTGACCTGTTCTACTCTCCGGACg                                                                                                                                                                         >  1:325725/1‑251 (MQ=255)
                                                                        aaacttATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAACACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTCTCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTACAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCAGCTGCCCTGTTCTTCTCTCCGGACGGCgcactgat                                                                                                                                                               >  2:273699/5‑244 (MQ=255)
                                                                           aTTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGCCCTGTTCTACTCTCCGGACGGCGAAATGTTCgc                                                                                                                                                            >  1:237976/1‑250 (MQ=255)
                                                                           aTTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCg                                                                                                                                                             >  1:302058/1‑250 (MQ=255)
                                                                              aTCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGt                                                                                                                                                         <  1:267503/251‑1 (MQ=255)
                                                                              aTCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGt                                                                                                                                                         >  1:496489/1‑251 (MQ=255)
                                                                                          tgcTGTGGCGCAGTTCTGCCGCACGGAAAAATGCGCGAAGTCGTTATCGGGATGGGTGCATGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAAGGTTCGACGTGATGGAAAAAt                                                                                                                                              <  2:659433/248‑1 (MQ=255)
                                                                                                ggCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTAGCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCTTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCa                                                                                                                                        <  2:385104/250‑1 (MQ=255)
                                                                                                  cgcAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACATGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGGAATGTTCGACGGGATGGAAAAATACGGcctcc                                                                                                                                     >  2:394583/1‑247 (MQ=255)
                                                                                                  cgcAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTAGCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCAtcc                                                                                                                                     <  1:391212/251‑3 (MQ=255)
                                                                                                   gcAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGCAACCTGTTATACGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTCCGACGTGATGGAACAATACGGCAtccc                                                                                                                                    >  2:342470/1‑248 (MQ=255)
                                                                                                     aGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCg                                                                                                                                  <  1:115184/251‑1 (MQ=255)
                                                                                                     aGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCg                                                                                                                                  >  1:450076/1‑251 (MQ=255)
                                                                                                       ttCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGc                                                                                                                                >  1:50365/1‑251 (MQ=255)
                                                                                                           gCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTAGCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCGCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTg                                                                                                                             <  1:435591/250‑1 (MQ=255)
                                                                                                                aCGTAAAACTGCGCGACGTCGTTAGCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATg                                                                                                                       <  1:49427/251‑1 (MQ=255)
                                                                                                                    aaaaCTGCTCGACGTCTTTATCGGTATGGGTTCCTGCTCCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAg                                                                                                                   <  2:388744/251‑1 (MQ=255)
                                                                                                                    aaaaCTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGaa                                                                                                                    <  1:335086/250‑1 (MQ=255)
                                                                                                                     aaaCTGCGCGACTGCGGTATCGGTATGATTGCCTGCACCGATTCCATAGTTAACCTTATCCGTTTTAAAGACCATGGCTTTGCCGCTATCGCTGACTTGGACATGGTGCGTAACGCGGTAGACGCGGCGAAAGCACTGGGTGTTGATGCGCGCGTGGGTATCCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGc                                                                                                                  <  2:450076/251‑1 (MQ=255)
                                                                                                                              gACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTAt                                                                                                         >  1:122353/1‑251 (MQ=255)
                                                                                                                                cgtcgtTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCt                                                                                                       <  1:209549/251‑1 (MQ=255)
                                                                                                                                         cGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGCTGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGCATCTACGGCGTCg                                                                                              >  1:533805/1‑251 (MQ=255)
                                                                                                                                         cGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCg                                                                                              <  1:80385/251‑1 (MQ=255)
                                                                                                                                         cGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGAGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGACTGGGGTTTACGGCGTcc                                                                                              >  1:535636/1‑250 (MQ=255)
                                                                                                                                                ggTGCCTTCACCGATTCCAAAGTTAACCGTCTCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGa                                                                                       <  2:508958/251‑1 (MQ=255)
                                                                                                                                                 gTGCCTGCACCGATTCCAAAGTTAACCGTTTCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCTGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGag                                                                                      <  2:171848/251‑2 (MQ=255)
                                                                                                                                                    ccTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGttt                                                                                   >  1:43502/1‑251 (MQ=255)
                                                                                                                                                          ccGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGCAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGcgcg                                                                             >  1:300355/1‑251 (MQ=255)
                                                                                                                                                            gATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGa                                                                                       >  1:44297/1‑239 (MQ=255)
                                                                                                                                                            gATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGGTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGa                                                                                       <  2:44297/239‑1 (MQ=255)
                                                                                                                                                                 cAAAGTTAACCGTATCCATTTTAAAGATCATGACGTTGTCGCTATCGTTGACTTTGACAGGGTGAGTAACGCTGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCTGCTGATCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATTGAAGCTGCTGGTATCTACGGCGTAGATGCAGAGTTTGGCGCGAAAGccc                                                                      <  2:472841/251‑1 (MQ=255)
                                                                                                                                                                      ttAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTAct                                                                                                                                                                                  <  1:393334/138‑1 (MQ=255)
                                                                                                                                                                      ttAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATAGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTTTAct                                                                                                                                                                                  >  2:393334/1‑138 (MQ=255)
                                                                                                                                                                          ccgtaTCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGCCGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGCCGTGATGGACAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCTTCGCTGCAGAGTTTGGCGCGAAAGCCCTGCCCATCt                                                             >  2:462255/5‑251 (MQ=255)
                                                                                                                                                                                    tttAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTCACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAACGGCTGGTATTTACGGCGGTGCTGCAGAGTTTGGCCCGGAAGCCCTGACCctctcaccctggtc                                                   >  2:529508/1‑241 (MQ=255)
                                                                                                                                                                                    tttAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGtgtc                                                   >  1:370306/1‑248 (MQ=255)
                                                                                                                                                                                       aaaGTCCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTg                                                 <  1:86812/250‑1 (MQ=255)
                                                                                                                                                                                       aaaGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGAAACGCGGTAGACGCGGATAAAGCACTGGGTGTTGATGCTCGCGTGGGAAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAACATACGGAATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTAAGGCGCCGCCGCAGAGTTTGGCGCGAAAGCCCTGTCCATCCGCACCgggtctgc                                                >  1:471707/1‑244 (MQ=255)
                                                                                                                                                                                          gACCATGACTTGGCCGCTTACGCCGCTTTCGACAAGATGCGTAACGCCGTCGACGCGGCTAAAGCACTGGGGGATGATGCTCGCGTGGGTAACCTGGTCACAGTTGTCCCGTTCTACTCTCCGGCCGGCGAAATGTTCGACGTGATGGAAAAATAGGGCATCCCCGGCGCGGAAATGGAAGCGGCTGGTACCCTCGGCGTCGCTTCAGAGTTTGGCTCGAAAGTCCTGACCATCTGCTCCGTTTCTGTCca                                             <  2:35702/251‑1 (MQ=255)
                                                                                                                                                                                             cATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATTCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGGGGAAATGGAAGCGGCTGGTATCTAAGGCGTCACTGCAGAGTTTGGGGCGAAAACCCTGAACATCTGCACCgagtcaggccacac                                          >  2:599847/1‑238 (MQ=255)
                                                                                                                                                                                                      gCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTc                                  <  1:127473/250‑1 (MQ=255)
                                                                                                                                                                                                                  gtctttgACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCGACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCa                      <  2:50365/244‑1 (MQ=255)
                                                                                                                                                                                                                                  gtaccgctgTAGACGCGGCTAAAGCACTGGTTGTTTATGCCCGCGTGGGTAACCTGTTCTCCGCTGACCTTTTTTCCTCCCCGTCCGGCGTAATTGTCGCCGTGATGGAAAAATCCGGCATCCCCGGCGTGGAAATGGAAGCGTCTGGTATCTCCGGTGTCGCTTCTGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCTCTTCCGCCCTCACGAGCAGACCACTTCCGCTGAGCGTCAg     <  2:422671/243‑1 (MQ=255)
                                                                                                                                                                                                                                  gTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGCGCGTc       >  1:229620/1‑249 (MQ=255)
                                                                                                                                                                                                                                   tAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAAt                                                                                                                                              >  1:430738/1‑113 (MQ=255)
                                                                                                                                                                                                                                   tAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTACTACTCTCCGGACGTCGAAATGTTCGACGTGATGGAAAAAt                                                                                                                                              <  2:430738/113‑1 (MQ=255)
                                                                                                                                                                                                                                    aaCGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGAc   >  1:587402/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aCGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCTGACCATCTGCACCGTGTCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGAcc  <  1:462255/251‑2 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                
ATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCTGACCATCTGCACCGTATCTGACCACATCCGCACTCACGAGCAGACCACTGCCGCTGAGCGTCAGACT  >  NC_000913/4621063‑4621542

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