Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,316,203 T→G Q47P (CAG→CCG)  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,2030TG86.0% 672.4 / 97.9 286Q47P (CAG→CCG) phnKATP‑binding component of phosphonate transport
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (121/125);  minor base A (21/19);  total (142/144)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 288 total aligned reads displayed.

CGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCTGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCC  >  NC_000913/4316113‑4316292
                                                                                          |                                                                                         
cggTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGc                                                                                         >  2:2795254/1‑93 (MQ=255)
  gTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGt                                                                                       >  1:394558/1‑93 (MQ=255)
   tCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTg                                                                                      >  2:5737307/1‑93 (MQ=255)
    cGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTgg                                                                                     >  2:2022394/1‑93 (MQ=255)
      gCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTc                                                                                   <  1:5523360/93‑1 (MQ=255)
       ccTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCt                                                                                  <  2:3808448/93‑1 (MQ=255)
        cTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCtt                                                                                 >  2:3116800/1‑93 (MQ=255)
         tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTc                                                                                >  1:725700/1‑93 (MQ=255)
         tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTc                                                                                <  2:407292/93‑1 (MQ=255)
         tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTc                                                                                <  2:4911279/93‑1 (MQ=255)
         tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTc                                                                                <  2:5453797/93‑1 (MQ=255)
          cGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTcc                                                                               >  2:519573/1‑93 (MQ=255)
           gCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTccc                                                                              <  2:6231210/93‑1 (MQ=255)
            cTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCg                                                                             >  1:4858846/1‑93 (MQ=255)
            cTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCg                                                                             <  1:875886/93‑1 (MQ=255)
             tCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCgg                                                                            <  2:1102453/93‑1 (MQ=255)
             tCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCgg                                                                            >  2:4408570/1‑93 (MQ=255)
             tCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCgg                                                                            >  1:2211301/1‑93 (MQ=255)
             tCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCgg                                                                            <  2:7376650/93‑1 (MQ=255)
              cATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGa                                                                           >  2:2369723/1‑93 (MQ=255)
              cATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGa                                                                           <  1:5181234/93‑1 (MQ=255)
                ttGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGc                                                                         >  1:6948126/1‑93 (MQ=255)
                ttGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGc                                                                         <  1:4974347/93‑1 (MQ=255)
                ttGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGc                                                                         <  1:5557009/93‑1 (MQ=255)
                 tGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGcc                                                                        >  1:252446/1‑93 (MQ=255)
                 tGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCCCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGcc                                                                        <  1:7507437/93‑1 (MQ=255)
                   cATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCgg                                                                      >  1:5027454/1‑93 (MQ=255)
                    aTACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGa                                                                     >  1:93888/1‑93 (MQ=255)
                    aTACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGa                                                                     >  2:3504798/1‑93 (MQ=255)
                     tACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGAt                                                                    >  2:2612275/1‑93 (MQ=255)
                      aCAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGAtt                                                                   >  1:1201067/1‑93 (MQ=255)
                       cAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTc                                                                  >  1:4083/1‑93 (MQ=255)
                        aGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTcc                                                                 >  2:3867691/1‑93 (MQ=255)
                           gAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCg                                                              >  1:3830535/1‑93 (MQ=255)
                            aaCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGa                                                             >  1:4882048/1‑93 (MQ=255)
                            aaCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGa                                                             >  2:251467/1‑93 (MQ=255)
                            aaCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGa                                                             >  2:7249144/1‑93 (MQ=255)
                            aaCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGa                                                             >  1:1035496/1‑93 (MQ=255)
                             aCGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGAc                                                            <  2:6807030/93‑1 (MQ=255)
                                 ttCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATg                                                        >  2:2228520/1‑93 (MQ=255)
                                  tctcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGc                                                       >  2:2732784/1‑93 (MQ=255)
                                   ctcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGcc                                                      <  2:5429747/93‑1 (MQ=255)
                                    tcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGccc                                                     <  2:7100783/93‑1 (MQ=255)
                                    tcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGccc                                                     <  1:1088135/93‑1 (MQ=255)
                                     cGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCa                                                    <  1:1841318/93‑1 (MQ=255)
                                      gTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCGGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAg                                                   <  1:5092855/93‑1 (MQ=255)
                                      gTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAg                                                   >  1:5613007/1‑93 (MQ=255)
                                        aGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCa                                                 <  2:703473/93‑1 (MQ=255)
                                        aGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCa                                                 <  1:5323642/93‑1 (MQ=255)
                                        aGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCa                                                 <  2:4141982/93‑1 (MQ=255)
                                         gTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCAc                                                <  2:6689725/93‑1 (MQ=255)
                                            aaTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGCCTTCAGCGGCGTGGTATTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTc                                             <  1:3324161/93‑1 (MQ=255)
                                             aTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTcc                                            <  1:1677718/93‑1 (MQ=255)
                                              tttCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTccc                                           >  1:538811/1‑93 (MQ=255)
                                               ttCCCCCTGCTGCGGCGTCAGGCGCGCGGAAATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTcccc                                          <  1:3186666/93‑1 (MQ=255)
                                                tCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTccccc                                         >  1:2903233/1‑93 (MQ=255)
                                                  ccccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCgg                                       >  2:1199445/1‑93 (MQ=255)
                                                  ccccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCgg                                       >  2:2471553/1‑93 (MQ=255)
                                                     ctgctgCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                                    >  1:6285896/1‑93 (MQ=255)
                                                     ctgctgCGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                                    >  1:6604460/1‑93 (MQ=255)
                                                       gctgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATa                                  >  2:1521599/1‑93 (MQ=255)
                                                       gctgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATa                                  <  2:5579051/93‑1 (MQ=255)
                                                         tgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATaaa                                <  1:2332100/93‑1 (MQ=255)
                                                         tgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATaaa                                <  1:452657/93‑1 (MQ=255)
                                                         tgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATaaa                                >  1:7179830/1‑93 (MQ=255)
                                                         tgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATaaa                                >  1:1495539/1‑93 (MQ=255)
                                                               gTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaaaa                          >  2:5238028/1‑93 (MQ=255)
                                                               gTCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaaaa                          <  1:3504798/93‑1 (MQ=255)
                                                                tCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAag                         <  2:1732884/93‑1 (MQ=255)
                                                                tCAGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAag                         >  1:4447948/1‑93 (MQ=255)
                                                                  aGGCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAagag                       >  1:6463109/1‑93 (MQ=255)
                                                                   ggCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAgaga                      <  1:4682152/93‑1 (MQ=255)
                                                                   ggCGCGCGGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAgaga                      <  1:2023578/93‑1 (MQ=255)
                                                                   ggCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAgaga                      <  1:3094324/93‑1 (MQ=255)
                                                                     cgcgcgGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACa                    <  1:2401994/93‑1 (MQ=255)
                                                                     cgcgcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACa                    <  2:326150/93‑1 (MQ=255)
                                                                      gcgcgGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACAt                   <  1:6513211/93‑1 (MQ=255)
                                                                      gcgcgGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACAt                   <  2:4144161/93‑1 (MQ=255)
                                                                      gcgcgGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACAt                   <  2:2745800/93‑1 (MQ=255)
                                                                      gcgcgGAGATCGACTCCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGGCAATGCCCAGCGCCTCCCCCGGCCATAAATCAAAATAGACAt                   <  1:6102771/93‑1 (MQ=255)
                                                                       cgcgGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATc                  >  1:438375/1‑93 (MQ=255)
                                                                       cgcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATc                  <  2:7502788/93‑1 (MQ=255)
                                                                        gcgGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCg                 >  2:4241499/1‑93 (MQ=255)
                                                                         cgGAGATCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGc                <  1:4227848/93‑1 (MQ=255)
                                                                            agaTCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTaa             <  1:5975342/93‑1 (MQ=255)
                                                                            agaTCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTaa             <  2:4982438/93‑1 (MQ=255)
                                                                             nntCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTaaa            <  1:3717046/91‑1 (MQ=255)
                                                                              aTCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAg           >  2:1366414/1‑93 (MQ=255)
                                                                              aTCGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAg           <  2:2944096/93‑1 (MQ=255)
                                                                                cGACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGcc         >  2:5977398/1‑93 (MQ=255)
                                                                                 gACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCt        <  2:1802992/93‑1 (MQ=255)
                                                                                 gACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCt        <  2:1909551/93‑1 (MQ=255)
                                                                                 gACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCt        <  2:4830552/93‑1 (MQ=255)
                                                                                 gACTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCt        <  2:6405401/93‑1 (MQ=255)
                                                                                   cTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCttt      >  2:861732/1‑93 (MQ=255)
                                                                                   cTTCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCttt      <  1:6726850/93‑1 (MQ=255)
                                                                                    ttCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTg     <  2:2153300/93‑1 (MQ=255)
                                                                                    ttCAGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTg     >  1:986792/1‑93 (MQ=255)
                                                                                    ttCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTg     >  2:2347235/1‑93 (MQ=255)
                                                                                       aGCGGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGccc  <  2:5796186/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCTGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCC  >  NC_000913/4316113‑4316292

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