Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,926,772 C→A G155G (GGC→GGA ppnN → nucleotide 5'‑monophosphate nucleosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,926,7720CA100.0% 137.6 / NA 40G155G (GGC→GGAppnNnucleotide 5'‑monophosphate nucleosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (21/19);  total (21/19)

AATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCAT  >  NC_000913/2926539‑2927008
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            
aaTCCCCCGGGAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACg                                                                                                                                                                                                                              >  1:176984/1‑250 (MQ=255)
   cccccGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGaaaa                                                                                                                                                                                                                          <  1:183366/251‑1 (MQ=255)
            gagGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTAttt                                                                                                                                                                                                                 <  1:86263/251‑1 (MQ=255)
             agGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTGCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGGTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGGGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTg                                                                                                                                                                                                                <  2:315110/251‑1 (MQ=255)
                    tGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGcc                                                                                                                                                                                                          <  1:154020/250‑1 (MQ=255)
                         aTGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAATCGAGTATTTGTATGCCcgtcgc                                                                                                                                                                                                    >  1:82092/1‑251 (MQ=255)
                             gCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCcgtcgcgtcg                                                                                                                                                                                                >  1:140086/1‑251 (MQ=255)
                                          gcTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGa                                                                                                                                                                                              >  1:125780/1‑240 (MQ=255)
                                          gcTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGa                                                                                                                                                                                              <  2:125780/240‑1 (MQ=255)
                                                         ttGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGCAACCAGCTGGGCATGCGTGAGCTGa                                                                                                                                                                     >  2:79132/1‑250 (MQ=255)
                                                                   tCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCAcc                                                                                                                                                          >  1:215959/1‑251 (MQ=255)
                                                                       gCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGCCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCt                                                                                                                                                      >  2:52993/1‑251 (MQ=255)
                                                                        cGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCtg                                                                                                                                                     <  1:67434/251‑1 (MQ=255)
                                                                           gACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTgg                                                                                                                                                  >  2:222253/1‑251 (MQ=255)
                                                                              aTTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTcc                                                                                                                                               <  1:141499/251‑1 (MQ=255)
                                                                                        tttACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATg                                                                                                                                     >  2:103746/1‑251 (MQ=255)
                                                                                              ggCAAATCCATAACACCGTTCGTGTTCCCAACCGGAATCTCGACAACTCCGTCCACATCACTAACCTGGGCGTTTACATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAAGAGCTGCACCGGCTGTGGTCCGGGAGCGAGGGAAgcg                                                                                                                               <  2:142601/251‑1 (MQ=255)
                                                                                               gCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAgcgc                                                                                                                              >  1:316483/1‑251 (MQ=255)
                                                                                                                            aCCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCCTGTGAGTGAAACGCCACATATGGTGGTCTGCTGGGGAGGTCCCTCAATTAACGAAAACGAGTCTTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGCGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGACGCGCCGATGCACGGGGCTGCGGTCGGACAcgcg                                                                                                 >  1:142645/1‑251 (MQ=255)
                                                                                                                             ccTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCGTCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACgcgc                                                                                                <  1:187017/251‑1 (MQ=255)
                                                                                                                                    ctcGACAACTCCGTCCATATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGAGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCg                                                                                         >  2:52549/1‑251 (MQ=255)
                                                                                                                                        ccAACTCCGTCCACAGCACTATCCTGGTCTGTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTAc                                                                                     <  2:240615/250‑1 (MQ=255)
                                                                                                                                              ccGTCCACATCTCTATCCTGGGCTTTGCCATGTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCAAAATATGGATGGTTGCTGGGGAGGTCACTCAATTAACGAAAACGATTATATGTATGCCCGCCGGGTCGGAAACCAGCTGGGCCTGCGTGAGCGGAATTTCTGCACCGGCTGTGGTCCGGGAGCGATGGGAGCGTTTATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTTCAAAGAc                                                                               <  1:322100/251‑1 (MQ=255)
                                                                                                                                              ccGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGAc                                                                               >  2:22795/1‑251 (MQ=255)
                                                                                                                                               cGTCCCCATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCTGTCGGACACGCGCAGCAGCGTTACAAAGACa                                                                              <  2:215959/251‑1 (MQ=255)
                                                                                                                                               cGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCATCGTTa                                                                                      <  2:164126/243‑1 (MQ=255)
                                                                                                                                               cGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTa                                                                                      >  1:164126/1‑243 (MQ=255)
                                                                                                                                                 tccacttgaCTAACATGGTCTTTTCAATCTTGCGTCACGCTCGCGCGCTGCATGTGGGAGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCCATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGACTATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCGGCAGCGTTACAACGACAGt                                                                            >  2:139284/9‑251 (MQ=255)
                                                                                                                                                      aTCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGttt                                                                       <  1:93542/251‑1 (MQ=255)
                                                                                                                                                                  gTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGAGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGCGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGCAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGa                                                            >  1:216551/1‑250 (MQ=255)
                                                                                                                                                                       ttCCCTCTTGCGTATCCCCCGCGCGATGCATGTGGGTGAAGCGCTATATATGGTGGTCTGCTGGGGATGTCACTCAATTAACGAAAACGAGTATTTGTATGCCTGTCGCCTCGGAAACCAGCTGTGTCTGCGTCCGCTGAAAATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGAGCCGAGGAAAGGTGCTGCGGTCGGAGACGCGCAGCAGCGTTACTAAGTCAGTCGTTCTATTGGTATGACagag                                                       <  1:191058/250‑1 (MQ=255)
                                                                                                                                                                         ccATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCg                                                    >  1:77595/1‑251 (MQ=255)
                                                                                                                                                                                       gCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATcgccgc                                      >  1:134835/1‑251 (MQ=255)
                                                                                                                                                                                          cgcgcgCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCt                                     <  1:222253/249‑1 (MQ=255)
                                                                                                                                                                                            cgcgcTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAAAGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTTTTGGTATGACAGAGCCGTCGATTATCGCCGCTGAAc                                 >  2:202737/1‑251 (MQ=255)
                                                                                                                                                                                                      gtgGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAAccc                       <  1:223662/251‑1 (MQ=255)
                                                                                                                                                                                                              aGCGCCAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAACCCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCa               <  2:79126/251‑1 (MQ=255)
                                                                                                                                                                                                                   cAAATATGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGa           >  2:231806/1‑250 (MQ=255)
                                                                                                                                                                                                                      atatGGTGGTCTGCTGGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTg       <  1:103746/251‑1 (MQ=255)
                                                                                                                                                                                                                           gtggtCTGCTAGGGAGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGACACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGatcat  >  2:211649/1‑251 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                            
AATCCCCCGGAAGAGGCTTTTGTCGATGGGCGAATTATTCGCGCTTTGCAGGCCAACTTGTTCGCGGTCCTGCGTGACATTCTCTTCGTTTACGGGCAAATCCATAACACCGTTCGTTTTCCCAACCTGAATCTCGACAACTCCGTCCACATCACTAACCTGGTCTTTTCCATCTTGCGTAACGCTCGCGCGCTGCATGTGGGTGAAGCGCCAAATATGGTGGTCTGCTGGGGCGGTCACTCAATTAACGAAAACGAGTATTTGTATGCCCGTCGCGTCGGAAACCAGCTGGGCCTGCGTGAGCTGAATATCTGCACCGGCTGTGGTCCGGGAGCGATGGAAGCGCCGATGAAAGGTGCTGCGGTCGGACACGCGCAGCAGCGTTACAAAGACAGTCGTTTTATTGGTATGACAGAGCCGTCGATTATCGCCGCTGAACCGCCTAACCCGCTGGTCAACGAATTGATCAT  >  NC_000913/2926539‑2927008

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: