Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,401,893 G→A Q455* (CAA→TAA)  abgB ← p‑aminobenzoyl‑glutamate hydrolase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,401,8930GA100.0% 127.8 / NA 37Q455* (CAA→TAA) abgBp‑aminobenzoyl‑glutamate hydrolase subunit B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/15);  total (22/15)

ATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATAT  >  NC_000913/1401661‑1402131
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
aTAATCAAATAGATAAAGAGCAGAAAAGGATGGGGACCCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGccc                                                                                                                                                                                                                              <  1:99337/251‑1 (MQ=255)
aTAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCATGCAATAGccc                                                                                                                                                                                                                              >  2:7037/1‑251 (MQ=255)
      aaaTAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATc                                                                                                                                                                                                                        <  2:75024/251‑1 (MQ=255)
         tAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGGGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCa                                                                                                                                                                                                                       >  2:93036/1‑249 (MQ=255)
            aTAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGa                                                                                                                                                                                                                   >  1:384106/1‑250 (MQ=255)
                agagCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAgaga                                                                                                                                                                                                              <  1:93036/251‑1 (MQ=255)
                agagCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTGCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAgaga                                                                                                                                                                                                              <  1:291113/251‑1 (MQ=255)
                             aTGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGAAGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTGTCGGGATAGGGCAGTGATACTGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATATCCCTGAATCAAGGAAGAAATTGACTGTGGTtg                                                                                                                                                                                                 >  2:35728/1‑251 (MQ=255)
                                    aCCTTGTTACCAATTCTCTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTACTCGGTTGCTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGCGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCAATCTTGCAATCGACCTGAATCAATGAAGAGATTGACTGTGGTTGATGCCAt                                                                                                                                                                                          >  1:326892/1‑251 (MQ=255)
                                                aTTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTcgccgc                                                                                                                                                                              >  1:66419/1‑251 (MQ=255)
                                                          cccAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCAtt                                                                                                                                                                    >  1:361313/1‑251 (MQ=255)
                                                          cccAGCCATAGAGCTTCCCGGATTGGGAGGAAGACGGTATGGATGACATGCTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCAtt                                                                                                                                                                    <  2:66419/251‑1 (MQ=255)
                                                                   agagCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCCTAGTTTTCGCCGCCAGCAGCATTCCTTTATg                                                                                                                                                            >  2:15495/1‑250 (MQ=255)
                                                                                        acgccgGTATTGATTACATACTCATGGGCATTCCTCGGTTGTTGTTGTTGCGTTGTTTGTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCTTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATCGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTGCGCCCCTgg                                                                                                                                      <  2:17400/247‑1 (MQ=255)
                                                                                         cgacgGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGc                                                                                                                                     >  1:315381/1‑251 (MQ=255)
                                                                                                      gACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGcc                                                                                                                         >  1:87497/1‑250 (MQ=255)
                                                                                                           gctgaTGGGCAGTCCTCGGTTGTTGTGTTTGCGGTGGTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCTTGTCCGTTACTTGCTGATGCTCTTTTTAGCATTCTTGCAATAGCCCTGAATCAAGTAAGAGATTGACTGTGGGTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTGCGCCCCTGGCTAACCAGTTGCCACGTa                                                                                                                    <  2:451920/246‑1 (MQ=255)
                                                                                                           cataaggggCTTTCCTCGGTTGTTTTGTTTGCGTGGTTGTTGTTATTGTAAAGGTGACGGTGTCGCGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTTCTTGCTGATGCGCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAATAGATTGACTGTGTTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGGTGTTCGCCCCTGGCTAACCAGTTGCCACggat                                                                                                                   <  2:431972/245‑4 (MQ=255)
                                                                                                           aCTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTAt                                                                                                                   >  1:416212/1‑251 (MQ=255)
                                                                                                                           gggTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCGTGTTAGCATTCGTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGGTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGAc                                                                                                   <  2:222423/250‑1 (MQ=255)
                                                                                                                                                 gttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCCGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAacac                                                                             >  1:281419/1‑251 (MQ=255)
                                                                                                                                                 gttgttATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAaca                                                                              <  1:364202/250‑1 (MQ=255)
                                                                                                                                                      tATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTggg                                                                         <  1:391678/250‑1 (MQ=255)
                                                                                                                                                       aTTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCa                                                                       <  1:150964/251‑1 (MQ=255)
                                                                                                                                                            aaaGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAgg                                                                  >  1:251510/1‑251 (MQ=255)
                                                                                                                                                                gTGACGGTGCCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGCATCTTGGAAGAGATTGTCTGTGGTTGCTGCCATCGTTTTCGCCGCTAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGCCTAACCAGTTGCCACTTATGTAGCGGTGTACCGACCGCACAACTGGGCCTGAAACACTTGGCAACAGGCAGt                                                              >  2:196718/1‑251 (MQ=255)
                                                                                                                                                                         tCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTg                                                     >  1:136320/1‑251 (MQ=255)
                                                                                                                                                                                       aTAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATc                                       >  1:165260/1‑251 (MQ=255)
                                                                                                                                                                                          gggCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGCAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGt                                    >  2:357474/1‑251 (MQ=255)
                                                                                                                                                                                              aGTGATACGCCTGCGCGTCCGTTACTTGCTGATGCTCATGTTAGCATTCTTGCATTAGCCCTGAATCAACGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTACCACTGACGTCGCCGACATCAGTCGAt                                >  2:67642/1‑251 (MQ=255)
                                                                                                                                                                                                 gATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATGCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGcc                             <  1:16942/251‑1 (MQ=255)
                                                                                                                                                                                                  aTACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCg                            >  1:155046/1‑251 (MQ=255)
                                                                                                                                                                                                             tgtCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGACGACATCAGTCGATGCCGCAAGCACGtt                  >  1:277549/1‑250 (MQ=255)
                                                                                                                                                                                                              gtCCGTTACTTGCTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGGTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTa                 <  1:453160/250‑1 (MQ=255)
                                                                                                                                                                                                                          cTGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATTGGTGGCGGCat    >  1:112249/1‑251 (MQ=255)
                                                                                                                                                                                                                           tGATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCat    >  2:84816/1‑250 (MQ=255)
                                                                                                                                                                                                                            gATGCTCTTGTTAGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCatat  <  1:341943/251‑1 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
ATAATCAAATAGATAAAGAGCAGAAAAGGATGGGGAACCTTGTTACCAATTCTTTCGACCCAGCCATAGAGCTTCCCGGATTGGGAGGACGACGGTATGGATGACATACTCATGGGCATTCCTCGGTTGTTGTGTTTGCGTTGTTGTTGTTATTTTAAAGGTGACGGTGTCACGTTTTTCGGGATAGGGCAGTGATACGGTTGCGTGTCCGTTACTTGCTGATGCTCTTGTTGGCATTCTTGCAATAGCCCTGAATCAAGGAAGAGATTGACTGTGGTTGCTGCCATAGTTTTCGCCGCCAGCAGCATTCCTTTATGAGCAATAGATGTTCGCCCCTGGCTAACCAGTTGCCACGTATGTAGCGGTGTACCGACGGCAAAACAGGGGCTGAAACACTGGGCAACAGGCAGTTTCCAACTGACGTCGCCGACATCAGTCGATGCCGCAAGCACGTTATCGGTGGCGGCATAT  >  NC_000913/1401661‑1402131

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: