Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,316,203 T→A Q47L (CAG→CTG)  phnK ← ATP‑binding component of phosphonate transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,316,2030TA100.0% 742.3 / NA 197Q47L (CAG→CTG) phnKATP‑binding component of phosphonate transport
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (94/103);  total (94/103)

Only 100 of 207 total aligned reads displayed.

CGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCTGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCC  >  NC_000913/4316113‑4316292
                                                                                          |                                                                                         
cggTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCagcagc                                                                                         <  2:5717109/93‑1 (MQ=255)
 ggTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCg                                                                                        <  1:4108545/93‑1 (MQ=255)
  gTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGt                                                                                       >  2:5521566/1‑93 (MQ=255)
   tCGGCCTCGCTCATTGCCTACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTg                                                                                      <  2:5470850/93‑1 (MQ=255)
   tCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTg                                                                                      >  1:1418265/1‑93 (MQ=255)
     ggCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGt                                                                                    >  1:2451277/1‑93 (MQ=255)
      gCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTc                                                                                   <  2:199762/93‑1 (MQ=255)
      gCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTc                                                                                   >  1:324914/1‑93 (MQ=255)
      gCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTc                                                                                   <  1:3226156/93‑1 (MQ=255)
      gCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTc                                                                                   >  1:5272254/1‑93 (MQ=255)
      gCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTc                                                                                   >  1:2881274/1‑93 (MQ=255)
        cTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCtt                                                                                 >  2:3562734/1‑93 (MQ=255)
        cTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCtt                                                                                 >  1:1774168/1‑93 (MQ=255)
         tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTc                                                                                >  1:4049161/1‑93 (MQ=255)
         tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTc                                                                                <  1:5783369/93‑1 (MQ=255)
         tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTc                                                                                >  2:1495216/1‑93 (MQ=255)
           gCTCATTGCATATAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTccc                                                                              <  2:2844111/93‑1 (MQ=255)
           gCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTccc                                                                              <  1:4829976/93‑1 (MQ=255)
            cTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCg                                                                             >  1:5891825/1‑93 (MQ=255)
            cTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCg                                                                             <  1:1021279/93‑1 (MQ=255)
             tCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCATCGTGGTCTTcccag                                                                            >  2:3325626/1‑91 (MQ=255)
             tCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCgg                                                                            >  2:4064357/1‑93 (MQ=255)
             gcATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCgg                                                                            >  1:3853750/2‑93 (MQ=255)
              cATTGCATACGGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGAATTCAGCAGCGTGGTCTTCCCGGa                                                                           >  1:3497014/1‑93 (MQ=255)
              cATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGa                                                                           >  2:406977/1‑93 (MQ=255)
               aTTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAg                                                                          >  2:2794488/1‑93 (MQ=255)
                ttGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGc                                                                         >  1:1678346/1‑93 (MQ=255)
                    aTACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGa                                                                     <  2:3286715/93‑1 (MQ=255)
                    aTACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGa                                                                     <  2:1426381/93‑1 (MQ=255)
                        aGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTcc                                                                 <  2:198232/93‑1 (MQ=255)
                         gCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTccc                                                                <  1:943077/93‑1 (MQ=255)
                         gCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTccc                                                                <  1:5724076/93‑1 (MQ=255)
                         gCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTccc                                                                >  1:4621025/1‑93 (MQ=255)
                            aaCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGa                                                             >  1:500816/1‑93 (MQ=255)
                             aCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGAc                                                            >  1:2027088/1‑93 (MQ=255)
                               ggTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACaa                                                          <  2:3450669/93‑1 (MQ=255)
                                gTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAAt                                                         >  2:1341667/1‑93 (MQ=255)
                                  tctcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGc                                                       <  2:5594847/93‑1 (MQ=255)
                                   ctcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGcc                                                      >  1:2979166/1‑93 (MQ=255)
                                    tcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGccc                                                     >  1:510467/1‑93 (MQ=255)
                                         gTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGTAGCCGGATTCCCCGACAATGCCCAGCAc                                                <  1:5027831/93‑1 (MQ=255)
                                         gTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCACCGACAATGCCCAGCAc                                                >  1:3596833/1‑93 (MQ=255)
                                          tGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACt                                               >  1:5247391/1‑93 (MQ=255)
                                          tGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACt                                               >  2:3018773/1‑93 (MQ=255)
                                            aaTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTc                                             <  2:2115558/93‑1 (MQ=255)
                                            aaTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTc                                             <  1:6096999/93‑1 (MQ=255)
                                             aTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGACGGATTCCCCGACAATGCCCAGCACTTcc                                            >  2:311957/1‑93 (MQ=255)
                                              tttCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTccc                                           >  2:4602603/1‑93 (MQ=255)
                                              tttCCCACTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACCATGCCCAGCACTTccc                                           >  1:6051109/1‑93 (MQ=255)
                                               tgcccccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTcccc                                          <  2:4900720/91‑1 (MQ=255)
                                                tCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTccccc                                         <  1:5535989/93‑1 (MQ=255)
                                                 cccccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCg                                        >  2:2188660/1‑93 (MQ=255)
                                                  ccccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCgg                                       >  2:3817391/1‑93 (MQ=255)
                                                   cccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGc                                      <  2:2052241/93‑1 (MQ=255)
                                                   cccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGc                                      >  2:2961756/1‑93 (MQ=255)
                                                    ccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGcc                                     <  2:3523238/93‑1 (MQ=255)
                                                    ccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGcc                                     <  2:5353605/93‑1 (MQ=255)
                                                     ctgctgCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                                    <  2:5044765/93‑1 (MQ=255)
                                                     ctgctgCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                                    >  2:3375329/1‑93 (MQ=255)
                                                      tgctgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCAt                                   >  2:3984523/1‑93 (MQ=255)
                                                         tgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATaaa                                >  1:271576/1‑93 (MQ=255)
                                                         tgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATaaa                                <  2:4813878/93‑1 (MQ=255)
                                                          gcggcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAAt                               <  1:1640499/93‑1 (MQ=255)
                                                            ggcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCa                             <  1:198599/93‑1 (MQ=255)
                                                              cgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaaa                           >  1:838855/1‑93 (MQ=255)
                                                               gTCAGGCGCGAGTAGATCGACTTCAGCAGCGTGGTCTTCCGGGAGCCGGATCCCCCGGCAATGCCCAGCACTTCCCCCGGCCATAAATCaaaa                          <  2:2269610/93‑1 (MQ=255)
                                                                tCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAag                         >  2:4260384/1‑93 (MQ=255)
                                                                 taGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAaga                        <  2:4161299/92‑1 (MQ=255)
                                                                 cAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAaga                        >  1:2310778/1‑93 (MQ=255)
                                                                   ggCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAgaga                      <  1:5158209/93‑1 (MQ=255)
                                                                   ggCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGAAAATGCCCAGCACTTCCACCGGCCATAAATCAAAAgaga                      >  1:3957860/1‑93 (MQ=255)
                                                                      gcgcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACAt                   >  1:4464646/1‑93 (MQ=255)
                                                                      gcgcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACAt                   >  2:6003989/1‑93 (MQ=255)
                                                                       cgcgGAGATCGACTTCAGCAGGGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATc                  <  2:1735238/93‑1 (MQ=255)
                                                                        gcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCg                 >  1:3469368/1‑93 (MQ=255)
                                                                         cgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGc                <  1:4177976/93‑1 (MQ=255)
                                                                         cgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGc                <  2:2463546/93‑1 (MQ=255)
                                                                         cgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGc                <  1:4071557/93‑1 (MQ=255)
                                                                          ggAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCt               <  1:4341159/93‑1 (MQ=255)
                                                                           gagaTCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTa              >  2:1866061/1‑93 (MQ=255)
                                                                            agaTCGGCTTCAGCAGAGTGGTCTTCCCGGAGCGGGATTCACCGACAATGCCGCGCCGTTCCCCCGGCCATAAATCAAAAGAGACATCGCTaa             <  1:653462/93‑1 (MQ=255)
                                                                            agaTCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTaa             <  2:5976030/93‑1 (MQ=255)
                                                                             gaTCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCTACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTaaa            <  1:4140377/93‑1 (MQ=255)
                                                                              aTCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAg           <  2:2336129/93‑1 (MQ=255)
                                                                               tCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGc          >  1:3694943/1‑93 (MQ=255)
                                                                               tCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGc          >  1:3155693/1‑93 (MQ=255)
                                                                                cGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGcc         <  2:2572122/93‑1 (MQ=255)
                                                                                 gACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCt        <  2:4570299/93‑1 (MQ=255)
                                                                                 gACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCt        >  2:1655769/1‑93 (MQ=255)
                                                                                 gACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCt        >  2:6025263/1‑93 (MQ=255)
                                                                                 gACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCt        >  2:1213950/1‑93 (MQ=255)
                                                                                  aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       <  2:613603/93‑1 (MQ=255)
                                                                                  aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       >  2:661026/1‑93 (MQ=255)
                                                                                  aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       >  2:4465838/1‑93 (MQ=255)
                                                                                  aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       <  1:160407/93‑1 (MQ=255)
                                                                                  aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       <  1:1266391/93‑1 (MQ=255)
                                                                                   ctacagcagCGTGGTCTTCCCGGAGCCGGATGCCCCGACAATGCCCAGCACTTCCCCCGGCCAGAAATCAAAATAGACATCGCTAAAGCCttt      <  1:6142367/90‑1 (MQ=255)
                                                                                   cTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCttt      <  2:4849872/93‑1 (MQ=255)
                                                                                    ttCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTg     <  1:2184012/93‑1 (MQ=255)
                                                                                       agcagcGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGccc  >  2:3930383/1‑93 (MQ=255)
                                                                                          |                                                                                         
CGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCTGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCC  >  NC_000913/4316113‑4316292

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: