Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,621,236 C→T R118R (CGC→CGT deoD → purine nucleoside phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,621,2360CT93.6% 133.9 / ‑1.8 47R118R (CGC→CGTdeoDpurine nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (24/20);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

GTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCT  >  NC_000913/4620995‑4621475
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
gtgaagtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGtttt                                                                                                                                                                                                                                        >  1:155509/1‑251 (MQ=255)
  gaagtgaaCAACTTTCGCGGTAGGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTaa                                                                                                                                                                                                                                      <  2:277862/251‑1 (MQ=255)
  gaagtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGACGAACATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTa                                                                                                                                                                                                                                       >  2:3327/1‑250 (MQ=255)
                    ggTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAGAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCCGCACCGATTCCAAAGTTACCCGTATCCGTTTTAGAGACCATGGCTTTGcagc                                                                                                                                                                                                                    >  2:201000/1‑248 (MQ=255)
                       aTGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAACATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATAGGGGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTAt                                                                                                                                                                                                                 >  2:64720/1‑251 (MQ=255)
                         gCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCg                                                                                                                                                                                                               >  1:467820/1‑251 (MQ=255)
                                ttCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACt                                                                                                                                                                                                         <  1:128526/250‑1 (MQ=255)
                                  cACCGGGACTTACAAAGGACGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGCCtttg                                                                                                                                                                                                      >  2:196876/1‑249 (MQ=255)
                                          cTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCACGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAACTTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGAGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCCAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGAAATGGTg                                                                                                                                                                                              >  2:467104/1‑251 (MQ=255)
                                              cAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGt                                                                                                                                                                                           >  1:347927/1‑250 (MQ=255)
                                                      gCAAAATTTCCGTAATGGGTCACGGTATGGGTCTCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTAGCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTa                                                                                                                                                                                  <  1:259204/251‑1 (MQ=255)
                                                         aaaTTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGac                                                                                                                                                                               <  1:278511/251‑2 (MQ=255)
                                                                cGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGCCCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGacgcggcc                                                                                                                                                                         >  2:121784/1‑243 (MQ=255)
                                                                      gggTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCa                                                                                                                                                                   <  1:175207/250‑1 (MQ=255)
                                                                                  gggTATCCCGTCCTGCTCCATCTACACCAAAGAACAGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCTCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCTCCGATTCCAAAGTTAACCGTATCCGTTTTAAATACCATGACGTTGCCGCTATCGCTGATTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATg                                                                                                                                                      <  1:230529/251‑1 (MQ=255)
                                                                                                 cTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGGGTTGATGCTCGCGTGGGTAAcc                                                                                                                                       >  1:147850/1‑251 (MQ=255)
                                                                                                 cTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGGGTTGATGCTCGCGTGGGTAAcc                                                                                                                                       >  1:304555/1‑251 (MQ=255)
                                                                                                    cATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTTTTGATGCTCGCGTGGGTAACCTGt                                                                                                                                    <  1:347029/251‑1 (MQ=255)
                                                                                                                         gacgagatcgaccatcaaatcggcaaTCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTAct                                                                                                               <  1:355799/226‑1 (MQ=255)
                                                                                                                         cACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTAct                                                                                                               >  2:420575/1‑251 (MQ=255)
                                                                                                                              aTTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCg                                                                                                          >  1:381747/1‑251 (MQ=255)
                                                                                                                                   ggCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGGTCTCCGCTGGCCTGTTCTACTCTCCGGACgg                                                                                                     >  2:143762/1‑251 (MQ=255)
                                                                                                                                              aaTTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCg                                                                                          >  1:380352/1‑251 (MQ=255)
                                                                                                                                              aaTTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCg                                                                                          <  1:420575/251‑1 (MQ=255)
                                                                                                                                                 tATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGAc                                                                                        <  2:381747/250‑1 (MQ=255)
                                                                                                                                                     cgcgTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGAt                                                                                   <  1:143762/251‑1 (MQ=255)
                                                                                                                                                            gTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGaaaaa                                                                            <  1:161775/251‑1 (MQ=255)
                                                                                                                                                               ccTGGGGCGCAGTTCTGCCGCACGAAAAACGGCGCGACGTCGTTATCGGTATGGGGGCATGCACCGATTCCAAAGTTAACCGCAGCCGTGTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACAGGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGGGGGTAACCTGTTCTCCGCTGAACTGTTCTAATCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATa                                                                          <  2:155509/250‑1 (MQ=255)
                                                                                                                                                                  gtgGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGc                                                                      >  1:23434/1‑251 (MQ=255)
                                                                                                                                                                                ccGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTgg                                                         >  1:164849/1‑250 (MQ=255)
                                                                                                                                                                                     cGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATgg                                                   <  1:201000/251‑1 (MQ=255)
                                                                                                                                                                                        aaaaCTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGCTTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGGAATGTTCGCCGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGagg                                                >  2:470667/1‑249 (MQ=255)
                                                                                                                                                                                         aaaCTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCCAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGCCGGCGAAATGTTCGACGTGATGGAAAAATACGGCATACTCGGCGTGGAAATTGaa                                                 >  2:95056/1‑249 (MQ=255)
                                                                                                                                                                                         aaaCTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGc                                               >  1:228044/1‑251 (MQ=255)
                                                                                                                                                                                         aaaCGGCGCGACGTCGTTAGCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGc                                               <  2:228044/251‑1 (MQ=255)
                                                                                                                                                                                            cggcgcGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTAGCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGc                                            <  2:36583/249‑1 (MQ=255)
                                                                                                                                                                                              gcgcGACGTCTTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGTGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGc                                            <  2:304555/249‑1 (MQ=255)
                                                                                                                                                                                               cgctaCGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTgg                                         <  2:23434/247‑1 (MQ=255)
                                                                                                                                                                                                    cgtcgtTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCt                                    >  1:186352/1‑251 (MQ=255)
                                                                                                                                                                                                          tATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGc                               <  2:186352/250‑1 (MQ=255)
                                                                                                                                                                                                           aTCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCg                              <  1:405699/250‑1 (MQ=255)
                                                                                                                                                                                                               gTATGGGTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCCTCCTCGGCGGGGAAATGGAAGCGGCTGGGATATACGGCGTCGCt                         >  2:307685/1‑251 (MQ=255)
                                                                                                                                                                                                                   gggTGCCTGCACCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCa                      <  1:76856/250‑1 (MQ=255)
                                                                                                                                                                                                                           gCACCGATTCCAAAGTTAACCTTATCCGTTTTAAAGTCCATGACTTTGTCGCTATCGCTGACTTTGTCATGGTGCGTAACGCGGGAGTCGCGGCTAAAGCACTGGGTGTTTATGCTCGCGTGGGTAACCTGGTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTTATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGgc             <  1:229477/251‑1 (MQ=255)
                                                                                                                                                                                                                             aCCGATTCCAAAGTTAACCGTATCCGTTTTAACGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGCAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGGATCTACGGCGGCGCGGCCGAGTTTGgcgc           >  2:270502/1‑251 (MQ=255)
                                                                                                                                                                                                                             aCCGATTCCAAAGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGgcgc           >  1:429131/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aaaGTTAACCGTATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTTGACATGGTGCGTAACGCGGTAGACGCGGCTAAAGCACTGGGTGTTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATCCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAGTTTGGCGCGAAAGCCCt  >  2:201175/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               
GTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCGGCGTGAAGAAAATTATCCGCGTGGGTTCCTGTGGCGCAGTTCTGCCGCACGTAAAACTGCGCGACGTCGTTATCGGTATGGGTGCCTGCACCGATTCCAAAGTTAACCGCATCCGTTTTAAAGACCATGACTTTGCCGCTATCGCTGACTTCGACATGGTGCGTAACGCAGTAGATGCAGCTAAAGCACTGGGTATTGATGCTCGCGTGGGTAACCTGTTCTCCGCTGACCTGTTCTACTCTCCGGACGGCGAAATGTTCGACGTGATGGAAAAATACGGCATTCTCGGCGTGGAAATGGAAGCGGCTGGTATCTACGGCGTCGCTGCAGAATTTGGCGCGAAAGCCCT  >  NC_000913/4620995‑4621475

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