Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,096,754 G→A G25D (GGT→GAT)  rdgB → dITP/XTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,096,7540GA100.0% 184.8 / NA 63G25D (GGT→GAT) rdgBdITP/XTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (40/23);  total (40/23)

GGCCAGGCCAACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGA  >  NC_000913/3096515‑3096991
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
ggccaggccaACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCCACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCg                                                                                                                                                                                                                                     >  2:627052/1‑251 (MQ=255)
         aaCAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGAGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAACTCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAACGACTTCGATCTTGATATCGTGGCCCaa                                                                                                                                                                                                                             >  2:119797/1‑250 (MQ=255)
         aaCAGTCATCTGGTAAAGTTTCCCGTAAAGAAATGCCGGGTTGCCAAAAGCCAGGTGTGGATTGAAAAAGGCAAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGGCGCTGCTTATCGACTTCGATCTTGATATCGTGGCCCaaa                                                                                                                                                                                                                            <  1:141923/251‑1 (MQ=255)
                  cTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGAATCCTCTGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCAGGCGTCGCTGCTTAGAGACTTCGATCTTGATATCGTGGCCCCAACAGACCt                                                                                                                                                                                                                     >  2:557513/1‑249 (MQ=255)
                    ggTGAAGTTTCCAGGTACTCAATGCCGGGTTCCCAAAAGCCAGGGGGTGATTGAAAAAGGCGAACGTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGc                                                                                                                                                                                                                 <  2:330133/251‑1 (MQ=255)
                     gTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCg                                                                                                                                                                                                                >  1:280395/1‑251 (MQ=255)
                                              gggTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGATCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCt                                                                                                                                                                                       >  1:74403/1‑251 (MQ=255)
                                               ggTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTg                                                                                                                                                                                      >  1:20284/1‑251 (MQ=255)
                                                   gCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCt                                                                                                                                                                                  >  1:664378/1‑251 (MQ=255)
                                                            cAGGAGGTGATTGACAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGATGTCCTCGCAACCAGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATATTGATATCTTGGCCCAAACAGACCGCGGCGTTGATTCCGCTGCAGAACCCGGCCTGACCTTTATCGaa                                                                                                                                                                          >  2:532233/1‑250 (MQ=255)
                                                             aGGTGGAGATTGAAAAAGGCGAAATTGGCCGCCACAAACACATTAAAATCATTAATCCGCAACAAATCCCACAAGAAATCACGGCGGTAATTAATTAGGTATCCTATGCAACAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTTAGCTGGCATCGCTGCTTAGCGACTTCGATCTTGATCTCGTGGCACAAACAGACCTCGGCGTTGATTCCGCTGAAGAACCCGGCCGGACCTTTATCGaaa                                                                                                                                                                         >  2:402624/1‑250 (MQ=255)
                                                               gtggggATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATACCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCTCTGCTTAGCGACTTCGATATTGATATCGTGGCCCTAACAGACATCGGCGTTGATTCCGCTGAAGAAAACGTCCTTACCTTTATCGAAAAcg                                                                                                                                                                      >  2:488069/1‑251 (MQ=255)
                                                                 ggtgATTGAAAAAGGCGACCTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAAcgcg                                                                                                                                                                    <  1:582594/251‑1 (MQ=255)
                                                                   tgATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGAt                                                                                                                                                                  >  1:313551/1‑251 (MQ=255)
                                                                       ggAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTg                                                                                                                                                              <  1:16452/250‑1 (MQ=255)
                                                                          aaaaGGCAAACTTGTCCGCCACAAACAAATTAAAAGCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGaaa                                                                                                                                                           <  2:264156/251‑1 (MQ=255)
                                                                            agggggagctgggCCGGCACAAACAAATGAAAATCAGTAACCCGGAACAAATCTCGCCAGAAAGCGGGGCGTTAAATAAGTAGGTATGCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGGCGCTGCGTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGGTGATTCCGCGGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAgc                                                                                                                                                         <  2:66018/240‑1 (MQ=255)
                                                                                   aCTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTgaaga                                                                                                                                                                                                >  1:597503/1‑205 (MQ=255)
                                                                                   aCTTGGCCGCCAAAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTgaaga                                                                                                                                                                                                <  2:597503/205‑1 (MQ=255)
                                                                                    cTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAAAAAATCCCGACAGAAATCGAGGCGTTAATTACTTAGATATCCTATGAAAAACGTTGTACTAGCAACCGGCAATGTCGGTAAAATGCGTGAAATGGCATAGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAATACACGGCGTTAATTCCGCTGACGAAAACGCACTGACCTTTCTCGAACACGCGGTTCTGaaacctcgcc                                                                                                                                                    >  2:635399/1‑241 (MQ=255)
                                                                                     cgggCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCTCGCCAt                                                                                                                                                  <  2:617287/247‑1 (MQ=255)
                                                                                         ccgccACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑Gcc                                                                                                                                           >  2:168409/1‑250 (MQ=255)
                                                                                                aaacaaaTTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTg                                                                                                                                     >  1:102743/1‑250 (MQ=255)
                                                                                                   caaaTTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCg                                                                                                                                 >  1:412341/1‑251 (MQ=255)
                                                                                                         aaaaTCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGGTTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTAc                                                                                                                           >  1:252660/1‑251 (MQ=255)
                                                                                                                          acaaATCCCGCCCGAAATCGCGGCGTTAATTAATTAGGTATCATATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGTGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCgccgac                                                                                                          >  1:516720/1‑251 (MQ=255)
                                                                                                                          acaaATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGATGTACTCGCAACCGGAAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTCCCGGCACTTGCCgacgac                                                                                                          >  2:85340/1‑251 (MQ=255)
                                                                                                                             aaTCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGActct                                                                                                       >  1:293581/1‑251 (MQ=255)
                                                                                                                                 ccGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGCCCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCAGACGACtctggtc                                                                                                   >  1:668385/1‑251 (MQ=255)
                                                                                                                                       gAAATAGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAACAGGAGATTCTGAAAGCGCGCAATGCG‑GCAAAAGTGAACGCTTTACCGGCATATGCCGCAGACCCTGGTCAGGCgg                                                                                             >  2:306574/1‑251 (MQ=255)
                                                                                                                                        aaaTCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGAACTCGGCGTTGCTTCCGCTGAGGAAACCGGCCTGACCTTTATCGAAGACGCGATTCTGAAAGCCCGCAATGCG‑CCAAAAGGGACCGCTTTACCGGCAATTGCCGCCGACTCTGGTCTGGCggc                                                                                            >  1:380311/1‑250 (MQ=255)
                                                                                                                                              cggcgTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGAtgt                                                                                      <  1:627052/251‑1 (MQ=255)
                                                                                                                                                 ggTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCGCGCAACCGGAAATGTCGGTAAAGTGCGTGAGCTGGCGGCGCGGCTAAGCGACTGCGATCTTGATATCGTGGCCCGAACAGACCCCGGCGTTTATTCCGCTGAAGGAGCGGGCCTTGCCTTTATCGAAAACGGGAGTCTGGACGCGCGTCATGCG‑GCAAGAGTGGGCGCTTTACCGTCAATGGCCGCTGACTCGGGTCTGGCGGTAGGTGTGCt                                                                                   <  2:644168/250‑1 (MQ=255)
                                                                                                                                                     aattaattAGGTATCCTATGCAAAAAGTCGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGTCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTggc                                                                               <  1:402624/251‑1 (MQ=255)
                                                                                                                                                     aattaattAGGAATCCTATGCAAAAAGTTGTACTCGCAACAGGCCATGTCGGTCAAGTGCGTGACCTGGAGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCACACAGACATCGGCGTTGATTCCGATGAAGAAACCGGCATGACATTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGTCCGCTTTCCCAGCAATTGCCGACGACTCTGGTCTAGCGGGAGATGCGCtttgg                                                                               >  2:670190/1‑248 (MQ=255)
                                                                                                                                                       ttaattaGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTggcgg                                                                             >  1:536112/1‑251 (MQ=255)
                                                                                                                                                              ggTATCCTATGCAAAAAGTTGTCCTCGCAACCGGAAATGTCGGTAAAGTGCGTGAGCTGGCGTCACTGCTTAACGCCTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGCAACCGGCCTGCCCTATATCGAACACCCTATTCTGAAAGCGCGACATGCGCGCAAA‑GAGACCGCTTTCACGACACTTACCGCCGACTCTGGTCTAGCgtggaagtgtctttgagagcagacg                                                                      >  2:404685/1‑227 (MQ=255)
                                                                                                                                                               gTATCCTATGCAAAAAGTTGTACTCGCAACCGGCCATGTCGGTAAAGCGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCACTTGCCGACGCCTCTGGTCTGGAGGTAGATGTGCTTGGCAGCGCGCCgg                                                                     >  2:151322/1‑251 (MQ=255)
                                                                                                                                                                        gCAACAAGTTGACCTCGCAACCGGCACTGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGTCCTTTATAGCAAACGATATTCTGAAAGCGCGCCATGCG‑GAAAAAGTGACCGCCTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGCAGTGCTTGGCCGGGCGCCCGGCGtttga                                                             >  2:650654/1‑248 (MQ=255)
                                                                                                                                                                        gCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGGCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCTTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACt                                                            <  2:371517/251‑1 (MQ=255)
                                                                                                                                                                         cAAAAAGTTGTACTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAACCGCGCTTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGCCTCTGGTCTGGCGGGTGATGTGCTTGGCGGCGCGCCGGGCCTTTACTc                                                           >  2:447835/1‑251 (MQ=255)
                                                                                                                                                                            aaaGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTcc                                                          >  1:238196/1‑249 (MQ=255)
                                                                                                                                                                            aaaGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCcgc                                                        >  1:188045/1‑251 (MQ=255)
                                                                                                                                                                            aaaGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAAATGTGCTTTGCGGCGCGCCGGGCATTTACTCcgc                                                        >  2:537552/1‑251 (MQ=255)
                                                                                                                                                                            aaaGTTGTCCTCGCAACCGGCAATGTCGGTAAAGGGCGTGAGCTGGCGTCGCTGCTTAGCGACTACGATCTTGAGATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTcc                                                          <  2:238196/249‑1 (MQ=255)
                                                                                                                                                                                       cGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTcc                                              <  2:357207/250‑1 (MQ=255)
                                                                                                                                                                                         cAACCGGCAATGGCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGt                                           <  2:293581/251‑1 (MQ=255)
                                                                                                                                                                                            ccGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCgggga                                         >  1:185111/1‑247 (MQ=255)
                                                                                                                                                                                               gCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAc                                     >  1:378097/1‑251 (MQ=255)
                                                                                                                                                                                               gCAATGTCGGTAAAGTGCGTGAGCTGGCGGCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAc                                     <  1:29771/251‑1 (MQ=255)
                                                                                                                                                                                                 aaTGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCGTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAcgc                                   <  1:168409/251‑1 (MQ=255)
                                                                                                                                                                                                  aTGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAcgcg                                  >  1:619941/1‑251 (MQ=255)
                                                                                                                                                                                                    gTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGAc                                >  1:415704/1‑251 (MQ=255)
                                                                                                                                                                                                      cGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACGCCGCGCGTTATTCCGGTGAAGACGCGAcc                               <  1:205292/250‑1 (MQ=255)
                                                                                                                                                                                                      cGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACCTCGATCTTGAGATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGAcc                               <  2:536112/250‑1 (MQ=255)
                                                                                                                                                                                                        gTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGTCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGAt                            <  2:102743/251‑1 (MQ=255)
                                                                                                                                                                                                          aaaGTGCGTGAGCTGGCGGCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTAGTCCGGTGAAGACGCGACCGATCa                          <  2:280395/251‑1 (MQ=255)
                                                                                                                                                                                                                cGTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGCCGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCACGCGTTATTCCGGTGAAGACGCGTCCGATCCAAAgc                     >  2:401347/1‑249 (MQ=255)
                                                                                                                                                                                                                 gTGAGCTGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGAAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGGTAAGACGCGACCGATCAAAAGAAt                   >  2:155933/1‑251 (MQ=255)
                                                                                                                                                                                                                       tGGCGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCa              >  1:525344/1‑250 (MQ=255)
                                                                                                                                                                                                                          cGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCAt                                                                 >  1:136446/1‑196 (MQ=255)
                                                                                                                                                                                                                          cGTCGCTGCTTAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCAt                                                                 <  2:136446/196‑1 (MQ=255)
                                                                                                                                                                                                                                  cgtAGCGACTTCGATCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGa  <  2:537145/249‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
GGCCAGGCCAACAGTCATCTGGTGAAGTTTCTCGGTAAGCAATTCCGGGTTGCCAAAAGCCAGGTGGTGATTGAAAAAGGCGAACTTGGCCGCCACAAACAAATTAAAATCATTAATCCGCAACAAATCCCGCCAGAAATCGCGGCGTTAATTAATTAGGTATCCTATGCAAAAAGTTGTCCTCGCAACCGGCAATGTCGGTAAAGTGCGTGAGCTGGCGTCGCTGCTTAGCGACTTCGGTCTTGATATCGTGGCCCAAACAGACCTCGGCGTTGATTCCGCTGAAGAAACCGGCCTGACCTTTATCGAAAACGCGATTCTGAAAGCGCGCCATGCG‑GCAAAAGTGACCGCTTTACCGGCAATTGCCGACGACTCTGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGA  >  NC_000913/3096515‑3096991

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: ‑