Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,618,604 T→C G51G (GGT→GGC deoD → purine‑nucleoside phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,618,6040TC86.0% 551.3 / 75.3 228G51G (GGT→GGCdeoDpurine‑nucleoside phosphorylase
Reads supporting (aligned to +/- strand):  ref base T (11/21);  new base C (84/112);  total (95/133)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 229 total aligned reads displayed.

GGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCG  >  NC_000913/4618512‑4618695
                                                                                            |                                                                                           
ggCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGt                                                                                             >  2:3776097/1‑93 (MQ=255)
  cGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCAc                                                                                           <  2:1357024/93‑1 (MQ=255)
   gACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACt                                                                                          >  2:2809836/1‑93 (MQ=255)
    aCCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACtt                                                                                         <  2:4644978/93‑1 (MQ=255)
      ccGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTAc                                                                                       >  1:1844050/1‑93 (MQ=255)
        gCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACaa                                                                                     >  1:4367999/1‑93 (MQ=255)
         cTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACaaa                                                                                    <  1:6514855/93‑1 (MQ=255)
         cTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACaaa                                                                                    <  1:4950821/93‑1 (MQ=255)
          tgcgtgcgAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAg                                                                                   <  2:1047130/93‑1 (MQ=255)
            cgtgcgAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGc                                                                                 <  1:1717634/93‑1 (MQ=255)
            cgtgcgAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGc                                                                                 <  2:6259391/93‑1 (MQ=255)
            cgtgcgAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGc                                                                                 <  1:2654748/93‑1 (MQ=255)
                  aaGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCaaa                                                                           >  2:5348233/1‑93 (MQ=255)
                  aaGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCaaa                                                                           >  2:2603371/1‑93 (MQ=255)
                   aGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCaaaa                                                                          <  2:5809045/93‑1 (MQ=255)
                   aGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCaaaa                                                                          >  2:7442557/1‑93 (MQ=255)
                    gTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAAt                                                                         >  1:785941/1‑93 (MQ=255)
                      atatTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAAttt                                                                       >  1:5812604/1‑93 (MQ=255)
                      atatTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAAttt                                                                       >  2:5486741/1‑93 (MQ=255)
                      atatTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAAttt                                                                       <  1:1150599/93‑1 (MQ=255)
                       tatTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTc                                                                      <  2:4574874/93‑1 (MQ=255)
                         ttGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCg                                                                    >  1:2172920/1‑93 (MQ=255)
                          tGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGt                                                                   <  2:3997207/93‑1 (MQ=255)
                           gCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTa                                                                  >  1:4091745/1‑93 (MQ=255)
                           gCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTa                                                                  <  2:3146825/93‑1 (MQ=255)
                           gCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTa                                                                  >  1:4673704/1‑93 (MQ=255)
                            cTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTaa                                                                 >  2:1287504/1‑93 (MQ=255)
                            cTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTaa                                                                 <  1:1934019/93‑1 (MQ=255)
                            cTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTaa                                                                 <  2:6879380/93‑1 (MQ=255)
                             tGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAAt                                                                <  2:7183060/93‑1 (MQ=255)
                             tGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAAt                                                                <  1:3365158/93‑1 (MQ=255)
                             tGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAAt                                                                >  1:6399849/1‑93 (MQ=255)
                             tGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAAt                                                                <  2:5937757/93‑1 (MQ=255)
                             tGAAACATTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAAt                                                                >  2:4771393/1‑93 (MQ=255)
                              gAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATg                                                               >  1:6050968/1‑93 (MQ=255)
                              gAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATg                                                               >  1:4041334/1‑93 (MQ=255)
                              gAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATg                                                               >  1:5412479/1‑93 (MQ=255)
                              gAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATg                                                               <  1:709434/93‑1 (MQ=255)
                               aaaCTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATgg                                                              <  2:1024860/93‑1 (MQ=255)
                                aaCTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATggg                                                             <  2:3949883/93‑1 (MQ=255)
                                aaCTTTCCTTGAAGACGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATggg                                                             >  2:1621202/1‑93 (MQ=255)
                                 aCTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGt                                                            <  1:7131150/93‑1 (MQ=255)
                                   tttCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCa                                                          <  2:3852292/93‑1 (MQ=255)
                                   tttCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCa                                                          >  2:3343035/1‑93 (MQ=255)
                                   tttCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCa                                                          <  1:1997533/93‑1 (MQ=255)
                                     tCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACg                                                        <  1:529272/93‑1 (MQ=255)
                                      ccTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACgg                                                       >  2:7343436/1‑93 (MQ=255)
                                        ttGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTa                                                     <  1:3694848/93‑1 (MQ=255)
                                        ttGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTa                                                     >  2:60458/1‑93 (MQ=255)
                                          gAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATg                                                   <  1:1096368/93‑1 (MQ=255)
                                          gAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATg                                                   <  2:4264624/93‑1 (MQ=255)
                                           aaGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATgg                                                  >  2:1460176/1‑93 (MQ=255)
                                           aaGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATgg                                                  <  1:1533390/93‑1 (MQ=255)
                                            aGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATggg                                                 <  2:5205305/93‑1 (MQ=255)
                                            aGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATggg                                                 >  2:2880635/1‑93 (MQ=255)
                                             gATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGt                                                <  1:3540021/93‑1 (MQ=255)
                                              aTGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTa                                               <  2:4086603/93‑1 (MQ=255)
                                               tGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTAt                                              >  2:5108013/1‑93 (MQ=255)
                                                 cccGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATcc                                            >  2:6494562/1‑93 (MQ=255)
                                                 cccGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATcc                                            <  1:6472626/93‑1 (MQ=255)
                                                  ccGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATccc                                           >  2:1969018/1‑93 (MQ=255)
                                                   cGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCg                                          >  2:3653657/1‑93 (MQ=255)
                                                     tgaagtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTc                                        <  1:3856953/93‑1 (MQ=255)
                                                      gaagtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTcc                                       <  1:4316828/93‑1 (MQ=255)
                                                       aagtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCt                                      <  1:3028034/93‑1 (MQ=255)
                                                       aagtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCt                                      <  1:563305/93‑1 (MQ=255)
                                                         gtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGc                                    >  1:6233589/1‑93 (MQ=255)
                                                         gtgaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTGACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGc                                    >  1:7509277/1‑93 (MQ=255)
                                                           gaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTc                                  >  2:6848767/1‑93 (MQ=255)
                                                           gaaCAACGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTc                                  >  2:4431845/1‑93 (MQ=255)
                                                             acaacGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCa                                <  2:1754795/93‑1 (MQ=255)
                                                               aacGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATc                              <  1:2600086/93‑1 (MQ=255)
                                                               aacGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATc                              <  1:6023703/93‑1 (MQ=255)
                                                                acGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCt                             <  2:3054652/93‑1 (MQ=255)
                                                                acGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCt                             >  2:6587139/1‑93 (MQ=255)
                                                                acGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCt                             >  1:1470374/1‑93 (MQ=255)
                                                                acGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCt                             <  2:2125192/93‑1 (MQ=255)
                                                                acGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCt                             >  1:4490631/1‑93 (MQ=255)
                                                                 cGTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTa                            <  2:7230594/93‑1 (MQ=255)
                                                                  gTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTac                           <  2:5660287/93‑1 (MQ=255)
                                                                  gTTCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTac                           >  2:4902621/1‑93 (MQ=255)
                                                                   ttCGCGGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTaca                          <  2:6613615/93‑1 (MQ=255)
                                                                       cgGTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCaa                      <  1:7390994/93‑1 (MQ=255)
                                                                        ggTATGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCaaa                     >  2:6413267/1‑93 (MQ=255)
                                                                         gTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAg                    >  1:4144294/1‑93 (MQ=255)
                                                                            tGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAAc                 <  2:2631334/93‑1 (MQ=255)
                                                                            tGCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAAc                 <  1:5961847/93‑1 (MQ=255)
                                                                             gCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACt                <  1:1701597/93‑1 (MQ=255)
                                                                             gCTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACt                >  1:398045/1‑93 (MQ=255)
                                                                              cTGGGCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTg               <  2:1499006/93‑1 (MQ=255)
                                                                                gggCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGAt             <  2:3755554/93‑1 (MQ=255)
                                                                                gggCTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGAt             <  1:2452405/93‑1 (MQ=255)
                                                                                   cTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCAc          <  2:4476723/93‑1 (MQ=255)
                                                                                   cTTCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCAc          <  2:6279091/93‑1 (MQ=255)
                                                                                    ttCACCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCAcc         >  2:1157854/1‑93 (MQ=255)
                                                                                       aCCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGAt      <  2:2401131/93‑1 (MQ=255)
                                                                                       aCCGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGAt      >  1:630355/1‑93 (MQ=255)
                                                                                        ccGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGAtt     <  1:891545/93‑1 (MQ=255)
                                                                                         cGGCACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGAttt    >  1:97209/1‑93 (MQ=255)
                                                                                           gTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCg  <  1:5191021/93‑1 (MQ=255)
                                                                                            |                                                                                           
GGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCG  >  NC_000913/4618512‑4618695

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: