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,2030TA99.6% 798.8 / NA 223Q47L (CAG→CTG) phnKATP‑binding component of phosphonate transport
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (107/115);  minor base C (1/0);  total (108/115)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 246 total aligned reads displayed.

GGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCTGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCC  >  NC_000913/4316114‑4316292
                                                                                         |                                                                                         
ggTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCg                                                                                        <  1:2388858/93‑1 (MQ=255)
 gTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGt                                                                                       <  2:2106832/93‑1 (MQ=255)
 gTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGt                                                                                       <  1:4457726/93‑1 (MQ=255)
  tCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGGTCGACTTCGGCAGCGTg                                                                                      >  1:1997768/1‑93 (MQ=255)
  tCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTg                                                                                      <  1:2766617/93‑1 (MQ=255)
    ggCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGAGGt                                                                                    >  2:2818583/1‑93 (MQ=255)
     gCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTc                                                                                   <  2:2750890/93‑1 (MQ=255)
      ccTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCt                                                                                  >  2:1519090/1‑93 (MQ=255)
      ccTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCt                                                                                  >  1:1314194/1‑93 (MQ=255)
       cTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCtt                                                                                 >  1:3490331/1‑93 (MQ=255)
       cTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCtt                                                                                 <  2:3413977/93‑1 (MQ=255)
        tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTc                                                                                >  2:205141/1‑93 (MQ=255)
        tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTc                                                                                >  1:594654/1‑93 (MQ=255)
        tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTc                                                                                >  2:220068/1‑93 (MQ=255)
         cGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTcc                                                                               <  2:74242/93‑1 (MQ=255)
         cGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTcc                                                                               <  1:2494371/93‑1 (MQ=255)
          gCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTccc                                                                              >  2:4255158/1‑93 (MQ=255)
          gCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTccc                                                                              >  2:2879085/1‑93 (MQ=255)
          gCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTccc                                                                              >  1:3270566/1‑93 (MQ=255)
          ccTCATTGCATACAGCGAACTGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTccc                                                                              <  2:1374318/92‑1 (MQ=255)
           cTCATTACATACAGCGAACGTTTCTCGTAGTAAATTCCCCCCTGCGGCGGCGTCGGGCGCGGGGAGATCGACTTCAGCAGCGTGGTCTCCCCg                                                                             <  2:2399519/93‑1 (MQ=255)
            tCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCgg                                                                            >  2:356973/1‑93 (MQ=255)
             cATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGa                                                                           >  1:4313853/1‑93 (MQ=255)
             cATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGa                                                                           <  2:3750359/93‑1 (MQ=255)
             cATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGGCTTCCCGGa                                                                           >  2:1921406/1‑93 (MQ=255)
             cATTGCACACAGCGAACGGTTCTCGTAGTGAATTGCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGa                                                                           <  2:4386814/93‑1 (MQ=255)
              aTTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAg                                                                          >  1:3505116/1‑93 (MQ=255)
               ttGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGc                                                                         >  1:4615420/1‑93 (MQ=255)
                   aTACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGa                                                                     >  1:4197502/1‑93 (MQ=255)
                      cAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTc                                                                  <  2:3906361/93‑1 (MQ=255)
                       aGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTcc                                                                 >  1:1770171/1‑93 (MQ=255)
                       aGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTcc                                                                 <  2:4408840/93‑1 (MQ=255)
                       aGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTcc                                                                 >  1:3751337/1‑93 (MQ=255)
                        gCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTccc                                                                >  1:2039803/1‑93 (MQ=255)
                         cGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTcccc                                                               <  2:2716108/93‑1 (MQ=255)
                            aCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGAc                                                            >  1:2260023/1‑93 (MQ=255)
                              ggTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACaa                                                          <  1:2531852/93‑1 (MQ=255)
                              ggTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACaa                                                          <  2:2495238/93‑1 (MQ=255)
                              ggTTCTCGTAGTGAATTCCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACaa                                                          <  2:3871408/93‑1 (MQ=255)
                               gTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAAt                                                         >  1:3361754/1‑93 (MQ=255)
                                ttCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATg                                                        >  2:2655707/1‑93 (MQ=255)
                                   tcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGccc                                                     >  1:164640/1‑93 (MQ=255)
                                   tcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGccc                                                     >  1:830357/1‑93 (MQ=255)
                                   tcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGccc                                                     >  2:1812919/1‑93 (MQ=255)
                                     gTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAg                                                   <  1:995294/93‑1 (MQ=255)
                                      tAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGc                                                  >  1:3144114/1‑93 (MQ=255)
                                        gTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCAc                                                <  1:4079345/93‑1 (MQ=255)
                                        gTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCAc                                                <  2:2677362/93‑1 (MQ=255)
                                         tGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACt                                               <  1:3574425/93‑1 (MQ=255)
                                         tGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACt                                               <  2:3806271/93‑1 (MQ=255)
                                          gAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACtt                                              <  1:2733835/93‑1 (MQ=255)
                                           aaTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTc                                             <  1:7062/93‑1 (MQ=255)
                                           aaTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTc                                             <  2:1350676/93‑1 (MQ=255)
                                           aaTGTCCCCCTGCTGCGGCGTCAGGCGCGCGTAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCTACAATGCCCAGCACTTc                                             <  2:4073021/93‑1 (MQ=255)
                                            aTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTcc                                            >  1:2278816/1‑93 (MQ=255)
                                            aTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTcc                                            <  2:4494906/93‑1 (MQ=255)
                                                  cccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCCGc                                      >  1:921674/1‑93 (MQ=255)
                                                   ccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGcc                                     >  2:834155/1‑93 (MQ=255)
                                                    ctgctgCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                                    >  2:2344020/1‑93 (MQ=255)
                                                    ctgctgCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                                    >  1:2851509/1‑93 (MQ=255)
                                                    ctgctgCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                                    >  1:4419166/1‑93 (MQ=255)
                                                     tgctgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCAt                                   >  1:2181331/1‑93 (MQ=255)
                                                        tgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATaaa                                >  1:1722664/1‑93 (MQ=255)
                                                        tgcGGCGTCAGGCGCGCGGAGATCGACTACAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATaaa                                <  2:4168238/93‑1 (MQ=255)
                                                         gcggcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAAt                               <  2:4467446/93‑1 (MQ=255)
                                                          cggcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATc                              >  1:1360102/1‑93 (MQ=255)
                                                           ggcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCa                             <  2:1620405/93‑1 (MQ=255)
                                                           ggcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCa                             <  2:1702888/93‑1 (MQ=255)
                                                            gcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaa                            >  1:1274285/1‑93 (MQ=255)
                                                              gTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTCCCCGGAGCCGGATGCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaaaa                          <  2:3619235/93‑1 (MQ=255)
                                                               tCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAag                         >  1:2666530/1‑93 (MQ=255)
                                                                cAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAaga                        <  2:107536/93‑1 (MQ=255)
                                                                cAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAaga                        >  1:1594034/1‑93 (MQ=255)
                                                                cAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAaga                        >  2:1313625/1‑93 (MQ=255)
                                                                cAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAaga                        <  1:3657957/93‑1 (MQ=255)
                                                                 aGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAagag                       <  1:2214911/93‑1 (MQ=255)
                                                                  ggggcgcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAgaga                      <  1:1033305/90‑1 (MQ=255)
                                                                  ggCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAgaga                      <  2:3708534/93‑1 (MQ=255)
                                                                   gcgcgcGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGAc                     <  1:4662625/93‑1 (MQ=255)
                                                                    gacgcgGAGATCGCCATCAGCAGCGTGGTCTTCACGGAGCCGGATTCCCCGACAAGGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACa                    <  2:3316666/91‑1 (MQ=255)
                                                                    cgcgcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACa                    <  2:1254084/93‑1 (MQ=255)
                                                                     gcgcgGAGATCGACTTCGGCAGCGTGGTCTTCCCTGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACAt                   <  1:4169196/93‑1 (MQ=255)
                                                                        cgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGc                <  1:629925/93‑1 (MQ=255)
                                                                         ggAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCt               <  2:2241666/93‑1 (MQ=255)
                                                                         ggAGATCGACTTCAGCAGCGTGGTCTGCCCAGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCt               <  2:3091205/93‑1 (MQ=255)
                                                                           agaTCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTaa             <  1:2607000/93‑1 (MQ=255)
                                                                            gaTCGACTTCAGCAGCGTGGCCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTaaa            <  1:339557/93‑1 (MQ=255)
                                                                             aTCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAg           >  1:1394222/1‑93 (MQ=255)
                                                                                gACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCt        <  2:599509/93‑1 (MQ=255)
                                                                                 aCTTCAGCAGCGTGTTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       <  1:3001557/93‑1 (MQ=255)
                                                                                 aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       <  1:3951903/93‑1 (MQ=255)
                                                                                 aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       <  1:3824879/93‑1 (MQ=255)
                                                                                 aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       <  2:499474/93‑1 (MQ=255)
                                                                                 aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       >  2:566304/1‑93 (MQ=255)
                                                                                  cTTCAGCAGCGTGTTGCTCCCGGAGCCCGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCttt      <  2:350837/93‑1 (MQ=255)
                                                                                  cTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCttt      >  1:1696649/1‑93 (MQ=255)
                                                                                  cTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCttt      >  2:1999605/1‑93 (MQ=255)
                                                                                  cTTCAGCAGCGTGGTCTTCCCGGAGCCGGATCCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCttt      <  1:2866122/93‑1 (MQ=255)
                                                                                      agcagcGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGccc  <  1:1818787/93‑1 (MQ=255)
                                                                                      agcagcGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGccc  >  2:4654096/1‑93 (MQ=255)
                                                                                         |                                                                                         
GGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCTGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCC  >  NC_000913/4316114‑4316292

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