Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 267,354 A→G W218R (TGG→CGG)  yafZ ← CP4‑6 prophage; conserved protein YafZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913267,3540AG100.0% 1194.2 / NA 288W218R (TGG→CGG) yafZCP4‑6 prophage; conserved protein YafZ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (133/155);  total (133/155)

Only 100 of 289 total aligned reads displayed.

ACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCT  >  NC_000913/267256‑267451
                                                                                                  |                                                                                                 
 ccGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAg                                                                                                <  1:455881/101‑1 (MQ=255)
 ccGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAg                                                                                                <  1:456015/101‑1 (MQ=255)
   gcgAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGt                                                                                              <  1:456142/101‑1 (MQ=255)
    cgAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTc                                                                                             >  1:456375/1‑101 (MQ=255)
    cgAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTc                                                                                             <  2:456037/101‑1 (MQ=255)
     gAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCa                                                                                            <  1:455776/101‑1 (MQ=255)
      aaCGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCAt                                                                                           >  1:456379/1‑101 (MQ=255)
        cGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTg                                                                                         >  1:456383/1‑101 (MQ=255)
         ggCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGc                                                                                        <  1:455850/101‑1 (MQ=255)
         ggCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGc                                                                                        <  2:456002/101‑1 (MQ=255)
         ggCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGc                                                                                        <  1:456060/101‑1 (MQ=255)
           cACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGctc                                                                                      >  2:456384/1‑101 (MQ=255)
               ggTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCAt                                                                                  <  2:456131/101‑1 (MQ=255)
                gTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATc                                                                                 >  2:456389/1‑101 (MQ=255)
                   tgtgACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTg                                                                              <  2:455968/101‑1 (MQ=255)
                   tgtgACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTg                                                                              <  2:456237/101‑1 (MQ=255)
                   tgtgACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTg                                                                              >  1:456393/1‑101 (MQ=255)
                    gtgACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGc                                                                             <  1:455736/101‑1 (MQ=255)
                      gACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCa                                                                           >  1:456396/1‑101 (MQ=255)
                       aCCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAg                                                                          <  1:455889/101‑1 (MQ=255)
                         cGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCAt                                                                                  >  2:456402/1‑91 (MQ=255)
                         cGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAgcg                                                                        <  1:456044/101‑1 (MQ=255)
                         cGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAgcg                                                                        >  2:456399/1‑101 (MQ=255)
                         cGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCAGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAg                                                                          <  1:455929/99‑1 (MQ=255)
                            cctcctTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAgcggcg                                                                     >  1:456404/1‑101 (MQ=255)
                            cctcctTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAgcggcg                                                                     >  2:456407/1‑101 (MQ=255)
                               cctTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAgg                                                                  <  1:456127/101‑1 (MQ=255)
                               cctTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAgg                                                                  >  1:456410/1‑101 (MQ=255)
                                ctTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAggg                                                                 <  1:456319/101‑1 (MQ=255)
                                ctTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAggg                                                                 >  2:456411/1‑101 (MQ=255)
                                   tGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGgag                                                              >  2:456416/1‑101 (MQ=255)
                                    ggCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGgaga                                                             >  1:456417/1‑101 (MQ=255)
                                      cATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAgg                                                           >  2:456421/1‑101 (MQ=255)
                                      cATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTACTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAgg                                                           <  1:455934/101‑1 (MQ=255)
                                       aTTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGa                                                          >  2:456425/1‑101 (MQ=255)
                                        ttACGGCCACTTAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGAt                                                         >  2:456428/1‑101 (MQ=255)
                                           cGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTg                                                      <  2:455792/101‑1 (MQ=255)
                                            ggCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGc                                                     >  1:456431/1‑101 (MQ=255)
                                            ggCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGc                                                     <  2:456256/101‑1 (MQ=255)
                                            ggCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGc                                                     <  2:456043/101‑1 (MQ=255)
                                            ggCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGc                                                     >  2:456432/1‑101 (MQ=255)
                                              ccACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATc                                                                                 <  1:456071/71‑1 (MQ=255)
                                              ccACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGa                                                   <  2:456046/101‑1 (MQ=255)
                                               cACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGAt                                                  <  1:455981/101‑1 (MQ=255)
                                               cACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGAt                                                  >  1:456437/1‑101 (MQ=255)
                                                aCTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAgg                                                           <  2:456144/91‑1 (MQ=255)
                                                aCTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGAtt                                                 <  1:456088/101‑1 (MQ=255)
                                                 cTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTc                                                <  1:455948/101‑1 (MQ=255)
                                                  tGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCCTCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCa                                               <  1:456005/101‑1 (MQ=255)
                                                  tGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCa                                               >  1:456440/1‑101 (MQ=255)
                                                    aGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGt                                             >  2:456441/1‑101 (MQ=255)
                                                    aGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGt                                             <  1:456075/101‑1 (MQ=255)
                                                     gcccgcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTc                                            <  1:456156/101‑1 (MQ=255)
                                                       ccgcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGTAGAGGATCTGCGATTCAGTCAc                                          >  1:456445/1‑101 (MQ=255)
                                                       ccgcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCAc                                          <  2:455957/101‑1 (MQ=255)
                                                       ccgcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCAc                                          >  2:456446/1‑101 (MQ=255)
                                                         gcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCg                                        <  2:456013/101‑1 (MQ=255)
                                                         gcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCg                                        <  2:456108/101‑1 (MQ=255)
                                                         gcccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCg                                        <  2:455989/101‑1 (MQ=255)
                                                          cccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCgg                                       >  2:456450/1‑101 (MQ=255)
                                                          cccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCgg                                       <  2:455912/101‑1 (MQ=255)
                                                          cccTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCgg                                       <  2:455764/101‑1 (MQ=255)
                                                            cTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCt                                     >  1:456453/1‑101 (MQ=255)
                                                             ttAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAgcggc                                                                      <  2:455941/67‑1 (MQ=255)
                                                             ttAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCtg                                    >  2:456456/1‑101 (MQ=255)
                                                                aTCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCtggtg                                 <  2:456033/101‑1 (MQ=255)
                                                                aTCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCtggtg                                 >  2:456459/1‑101 (MQ=255)
                                                                 tCAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCtgg                                   >  2:456462/1‑98 (MQ=255)
                                                                  cAGGTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTggtggt                               >  1:456463/1‑101 (MQ=255)
                                                                     gTTCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCtt                            >  2:456468/1‑101 (MQ=255)
                                                                      ttCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTc                           <  1:455905/101‑1 (MQ=255)
                                                                      ttCTCCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTc                           >  1:456469/1‑101 (MQ=255)
                                                                        ctcCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCAc                         <  2:456059/101‑1 (MQ=255)
                                                                        ctcCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCAc                         <  1:455827/101‑1 (MQ=255)
                                                                         tcCTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCAcc                        >  1:456472/1‑101 (MQ=255)
                                                                          ccTGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCa                       <  2:456095/101‑1 (MQ=255)
                                                                            tGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCaaa                     <  1:456194/101‑1 (MQ=255)
                                                                            tGAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCaaa                     <  2:455856/101‑1 (MQ=255)
                                                                             gAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCg                                                    >  2:456479/1‑69 (MQ=255)
                                                                             gAATACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAg                    >  1:456477/1‑101 (MQ=255)
                                                                              aaTACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGc                   <  2:456378/101‑1 (MQ=255)
                                                                               aTACGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCg                  <  2:456177/101‑1 (MQ=255)
                                                                                  cGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCtggtgg                                >  1:456485/1‑84 (MQ=255)
                                                                                  cGCTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTa               >  2:456483/1‑101 (MQ=255)
                                                                                    cTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTAtg             >  1:456489/1‑101 (MQ=255)
                                                                                    cTGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTAtg             <  2:456308/101‑1 (MQ=255)
                                                                                     tGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCg                  >  2:456490/1‑95 (MQ=255)
                                                                                     tGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTAtgt            <  2:455891/101‑1 (MQ=255)
                                                                                     tGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTAtgt            <  1:456096/101‑1 (MQ=255)
                                                                                     tGGTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTAtgt            <  1:455970/101‑1 (MQ=255)
                                                                                      ggTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTAtgtg           <  2:456078/101‑1 (MQ=255)
                                                                                      ggTACGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTAtgtg           <  2:456051/101‑1 (MQ=255)
                                                                                         aCGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAgg        >  2:456494/1‑101 (MQ=255)
                                                                                         aCGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAgg        <  2:455831/101‑1 (MQ=255)
                                                                                          cGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAg                    >  2:456501/1‑88 (MQ=255)
                                                                                          cGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGa          >  1:456500/1‑98 (MQ=255)
                                                                                          cGTGGTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAggg       >  1:456498/1‑101 (MQ=255)
                                                                                             ggTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcgg    <  2:456021/101‑1 (MQ=255)
                                                                                             ggTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcgg    >  1:456505/1‑101 (MQ=255)
                                                                                              gTCCGCAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGggcggc   <  1:456226/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
ACCGCGAACGGCACGGGTATGTGACCGTCCTCCTTTGGCATTACGGCCACTGAGCCCGCCCTTAATCAGGTTCTCCTGAATACGCTGGTACGTGGTCCACAGGTCATTGCTCTCATCCTGCCAGCGGCGAGGGGAGAGGATCTGCGATTCAGTCACCGGCTGGTGGTCTTCACCAAAGCGGTATGTGAGGGCGGCT  >  NC_000913/267256‑267451

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