Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,417,189 A→G A66A (GCT→GCC recE ← Rac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,417,1890AG100.0% 237.5 / NA 65A66A (GCT→GCCrecERac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/30);  total (35/30)

CTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTC  >  NC_000913/1417042‑1417329
                                                                                                                                                   |                                                                                                                                            
ctcAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCg                                                                                                                                           >  2:110850/1‑151 (MQ=255)
 tcAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTcggcgg                                                                                                                                          >  1:391939/1‑151 (MQ=255)
  cAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCggg                                                                                                                                         >  1:370814/1‑151 (MQ=255)
   aaTCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCggg                                                                                                                                         >  2:659491/1‑150 (MQ=255)
   aaTCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGc                                                                                                                                        >  1:167237/1‑151 (MQ=255)
   aaTCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGc                                                                                                                                        >  2:334990/1‑151 (MQ=255)
        cAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTc                                                                                                                                   <  1:640195/151‑1 (MQ=255)
            cATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCAtt                                                                                                                                <  2:78272/150‑1 (MQ=255)
             aTGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTa                                                                                                                               >  1:462795/1‑150 (MQ=255)
             aTGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTa                                                                                                                               >  1:740426/1‑150 (MQ=255)
             aTGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTa                                                                                                                               >  1:606341/1‑150 (MQ=255)
               gTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACg                                                                                                                            <  1:20653/151‑1 (MQ=255)
                tCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACg                                                                                                                            >  1:672544/1‑150 (MQ=255)
                    cGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAgg                                                                                                                       >  2:387108/1‑151 (MQ=255)
                       cGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGaa                                                                                                                     >  2:248298/1‑150 (MQ=255)
                           gACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTc                                                                                                                <  1:160878/151‑1 (MQ=255)
                            aCGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCTCCGTCGGCGGGCAGGTGATTAACGACAGGAAAAtt                                                                                                                 <  2:356225/149‑1 (MQ=255)
                            aCGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCg                                                                                                               >  1:432771/1‑151 (MQ=255)
                             cGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGt                                                                                                              >  1:458979/1‑151 (MQ=255)
                               ttGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTg                                                                                                             <  1:262954/150‑1 (MQ=255)
                                    aTTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAAc                                                                                                        <  2:617525/150‑1 (MQ=255)
                                    aTTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACa                                                                                                       <  1:409678/151‑1 (MQ=255)
                                    aTTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACa                                                                                                       <  1:875394/151‑1 (MQ=255)
                                    aTTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACa                                                                                                       <  1:713681/151‑1 (MQ=255)
                                       gCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCTCCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGc                                                                                                     <  1:313419/150‑1 (MQ=255)
                                        ccTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCtt                                                                                                   >  1:497396/1‑151 (MQ=255)
                                        ccTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCtt                                                                                                   <  1:584852/151‑1 (MQ=255)
                                        ccTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCtt                                                                                                   >  2:660750/1‑151 (MQ=255)
                                               aGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAAt                                                                                             <  1:794200/150‑1 (MQ=255)
                                               aGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATa                                                                                            <  1:414011/151‑1 (MQ=255)
                                                gTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAAt                                                                                             >  2:357768/1‑149 (MQ=255)
                                                 tGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAgc                                                                                          >  1:267774/1‑151 (MQ=255)
                                                  gAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAgc                                                                                          <  1:119852/150‑1 (MQ=255)
                                                       ttGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTc                                                                                     >  2:66697/1‑150 (MQ=255)
                                                          tCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGtt                                                                                 >  2:257841/1‑151 (MQ=255)
                                                             ggTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGtt                                                                                 <  1:679017/148‑1 (MQ=255)
                                                                gctgctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTAc                                                                           >  1:782167/1‑151 (MQ=255)
                                                                 ctgctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTAcc                                                                          >  1:758546/1‑151 (MQ=255)
                                                                  tgctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCt                                                                         >  1:347656/1‑151 (MQ=255)
                                                                  tgctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCt                                                                         >  1:622805/1‑151 (MQ=255)
                                                                   gctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTAcc                                                                          <  1:787819/149‑1 (MQ=255)
                                                                     tCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGAc                                                                      >  2:618526/1‑151 (MQ=255)
                                                                      cccGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGAc                                                                      <  1:162237/150‑1 (MQ=255)
                                                                         ggTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTa                                                                   >  2:2502/1‑150 (MQ=255)
                                                                          gTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAAc                                                                 <  1:510621/151‑1 (MQ=255)
                                                                              tAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGAt                                                              <  1:38776/150‑1 (MQ=255)
                                                                               aGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCa                                                            <  2:669126/151‑1 (MQ=255)
                                                                                gTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATc                                                             >  2:535947/1‑149 (MQ=255)
                                                                                             ggAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGt                                               <  2:583403/150‑1 (MQ=255)
                                                                                                gTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGca                                           <  1:248324/151‑1 (MQ=255)
                                                                                                   tttGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGcac                                          >  1:585054/1‑149 (MQ=255)
                                                                                                         gAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGAt                                  >  2:139718/1‑151 (MQ=255)
                                                                                                            ttGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTc                                >  1:772305/1‑150 (MQ=255)
                                                                                                                  gCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCg                         <  2:819749/151‑1 (MQ=255)
                                                                                                                   cGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGt                        <  1:357315/151‑1 (MQ=255)
                                                                                                                   cGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGt                        >  2:319008/1‑151 (MQ=255)
                                                                                                                    gTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGtt                       >  2:198689/1‑151 (MQ=255)
                                                                                                                      tCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGttg                      <  1:330314/150‑1 (MQ=255)
                                                                                                                           ccAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGcc                >  2:399783/1‑151 (MQ=255)
                                                                                                                             aGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGccc               >  2:85248/1‑150 (MQ=255)
                                                                                                                              ggTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCcac             <  2:782116/151‑1 (MQ=255)
                                                                                                                               gTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCcac             >  2:235669/1‑150 (MQ=255)
                                                                                                                                    aTCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGgacg       <  1:215549/151‑1 (MQ=255)
                                                                                                                                        aTCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGt   >  2:735321/1‑151 (MQ=255)
                                                                                                                                            cACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTc  <  1:659540/148‑1 (MQ=255)
                                                                                                                                            cACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTAGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTTCCCACAGgacgac     <  2:318/145‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
CTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTC  >  NC_000913/1417042‑1417329

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