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% 211.8 / NA 60A66A (GCT→GCC) recERac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/29);  total (31/29)

CCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTC  >  NC_000913/1417041‑1417335
                                                                                                                                                    |                                                                                                                                                  
cctcAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGc                                                                                                                                                  >  2:680765/1‑151 (MQ=255)
   cAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCggg                                                                                                                                               <  1:662232/151‑1 (MQ=255)
    aaTCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGc                                                                                                                                              >  1:422243/1‑151 (MQ=255)
      tctcAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGc                                                                                                                                              >  2:799299/1‑149 (MQ=255)
         cAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTc                                                                                                                                         >  2:735280/1‑151 (MQ=255)
             cATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTa                                                                                                                                     <  2:455162/151‑1 (MQ=255)
                     cGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAgg                                                                                                                             >  1:367754/1‑151 (MQ=255)
                        cGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGaaa                                                                                                                          <  1:782233/151‑1 (MQ=255)
                         gTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGaaa                                                                                                                          >  2:564646/1‑150 (MQ=255)
                            gACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTc                                                                                                                      <  1:618845/151‑1 (MQ=255)
                                  ggTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGc                                                                                                                 >  2:796829/1‑150 (MQ=255)
                                       tGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAg                                                                                                            <  2:345644/150‑1 (MQ=255)
                                        gCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGc                                                                                                           <  1:409692/150‑1 (MQ=255)
                                        gCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCt                                                                                                          <  2:712699/151‑1 (MQ=255)
                                          cTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCttt                                                                                                        >  2:273975/1‑151 (MQ=255)
                                           ttGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCttt                                                                                                        >  1:460538/1‑150 (MQ=255)
                                            tGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTa                                                                                                       <  1:123146/150‑1 (MQ=255)
                                             gATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTaa                                                                                                      <  1:850742/150‑1 (MQ=255)
                                                aGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAAt                                                                                                   <  2:590249/150‑1 (MQ=255)
                                                      cGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAgctgct                                                                                            <  2:743663/151‑1 (MQ=255)
                                                       gTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTc                                                                                           >  1:309519/1‑151 (MQ=255)
                                                              ggTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGttttt                                                                                    >  1:399014/1‑151 (MQ=255)
                                                                 gctgctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCTGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTAc                                                                                 <  1:187487/151‑1 (MQ=255)
                                                                  ctgctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTAcc                                                                                >  2:163146/1‑151 (MQ=255)
                                                                    gctCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTg                                                                              <  1:670852/151‑1 (MQ=255)
                                                                          ggTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTa                                                                         >  2:433641/1‑150 (MQ=255)
                                                                           gTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAAc                                                                       >  1:13264/1‑151 (MQ=255)
                                                                           gTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAAc                                                                       >  2:259680/1‑151 (MQ=255)
                                                                            ttttAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGATTTTTACCTGACTTAACg                                                                      >  1:706677/1‑151 (MQ=255)
                                                                             tttAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACg                                                                      >  2:723906/1‑150 (MQ=255)
                                                                              ttAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGAt                                                                    >  2:477175/1‑151 (MQ=255)
                                                                               tAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGAt                                                                    <  1:83154/150‑1 (MQ=255)
                                                                                    cccATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAAGAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTa                                                              <  1:322314/151‑1 (MQ=255)
                                                                                     ccATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAgg                                                                >  2:297138/1‑148 (MQ=255)
                                                                                         gTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTcc                                                          <  1:405175/150‑1 (MQ=255)
                                                                                          tCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTcc                                                          >  1:225962/1‑149 (MQ=255)
                                                                                           cATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCaga                                                       <  2:590452/151‑1 (MQ=255)
                                                                                           cATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCag                                                        <  2:309543/150‑1 (MQ=255)
                                                                                              ggAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGt                                                     <  2:369910/150‑1 (MQ=255)
                                                                                              ggAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTg                                                    >  2:370620/1‑151 (MQ=255)
                                                                                                     ttGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAgg                                             >  1:112888/1‑151 (MQ=255)
                                                                                                      tGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAgg                                             <  2:433841/150‑1 (MQ=255)
                                                                                                               gATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCaa                                    >  2:448386/1‑150 (MQ=255)
                                                                                                                 tAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAAt                                 >  2:858260/1‑151 (MQ=255)
                                                                                                                   gCGTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCg                               <  1:368306/151‑1 (MQ=255)
                                                                                                                     gTTCACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGtt                             <  2:558071/151‑1 (MQ=255)
                                                                                                                        cACTCCAGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGttgct                          >  1:193990/1‑151 (MQ=255)
                                                                                                                             cAGGTAAAATCGATCTCACCTTCGGGGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGccc                     >  2:521985/1‑151 (MQ=255)
                                                                                                                              aGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGccc                     >  1:377459/1‑150 (MQ=255)
                                                                                                                              aGGTAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTCGCAACAGCTTTAAAAGAGCTGCTCAGTTTTTTACCTGACTTACCGATCAGGTAGTACAGAGTTGCACAGGTCGATTCAAAATCGTTGCTTGCCca                    <  1:619798/151‑1 (MQ=255)
                                                                                                                                 tAAAATCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCcac                   >  1:522027/1‑149 (MQ=255)
                                                                                                                                  aaaaTCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAg                 <  2:395530/150‑1 (MQ=255)
                                                                                                                                    aaTCGATCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGgac              >  2:378052/1‑151 (MQ=255)
                                                                                                                                         aTCTCACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGt         >  2:560346/1‑151 (MQ=255)
                                                                                                                                             cACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGgacgac           <  1:485710/145‑1 (MQ=255)
                                                                                                                                             cACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGt         <  2:577916/147‑1 (MQ=255)
                                                                                                                                             cACCTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAg      <  1:440654/150‑1 (MQ=255)
                                                                                                                                               ccTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGtt   <  2:78530/151‑1 (MQ=255)
                                                                                                                                               ccTTCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGt    <  1:723863/150‑1 (MQ=255)
                                                                                                                                                 ttCGGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTc  >  2:853870/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCTCAATCTCAGTCATGTCTTCGCCGTTGACGTTGGTATTGCCTTGATAGTGAGCGTTGTCTGGTGCTGCTCCCGGTTTTAGTTCCCATGTCATGGAGTCTTTGCTGAGTTGATAGCGTTCACTCCAGGTAAAATCGATCTCACCTTCAGCGGGCAGGTCATTAACGACAGGAAAATTCGTGGCAACAGCTTTAAAATAGCTGCTCAGTTTTTTACCTGACTTAACGATCAGGTAGTCCAGAGTGGCACAGGTCGATTCAAAATCGTTGCTTGCCCACAGGACGACGTCAGGTTC  >  NC_000913/1417041‑1417335

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