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% 924.5 / NA 255Q47L (CAG→CTG) phnKATP‑binding component of phosphonate transport
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (145/109);  minor base C (0/1);  total (145/110)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 276 total aligned reads displayed.

GGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCTGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCC  >  NC_000913/4316114‑4316292
                                                                                          |                                                                                         
 gTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGTGCGGAGATCGACTTCAGCAGCGt                                                                                       >  2:5612346/1‑93 (MQ=255)
 gTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGt                                                                                       >  1:4729259/1‑93 (MQ=255)
 gTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCCACTTCAGCAGCGt                                                                                       >  1:2749017/1‑93 (MQ=255)
  tCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTg                                                                                      <  2:1618958/93‑1 (MQ=255)
  tCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTg                                                                                      <  1:5196509/93‑1 (MQ=255)
   cGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTgg                                                                                     >  2:9103786/1‑93 (MQ=255)
   cGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTgg                                                                                     <  1:3681155/93‑1 (MQ=255)
     gCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCAACTTCAGCAGCGTGGTc                                                                                   <  1:4414554/93‑1 (MQ=255)
       cTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTAGTCt                                                                                  >  2:9980161/1‑93 (MQ=255)
       cTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCtt                                                                                 >  2:2514613/1‑93 (MQ=255)
        tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTc                                                                                >  1:8796848/1‑93 (MQ=255)
        tCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTc                                                                                <  1:5963742/93‑1 (MQ=255)
         cGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTcc                                                                               >  1:1998566/1‑93 (MQ=255)
          gCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTccc                                                                              >  1:9876123/1‑93 (MQ=255)
           cTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCg                                                                             >  2:8188113/1‑93 (MQ=255)
            tCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCgg                                                                            >  2:7064104/1‑93 (MQ=255)
            tCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCgg                                                                            >  2:1968466/1‑93 (MQ=255)
               ttGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGc                                                                         >  1:676098/1‑93 (MQ=255)
               ttGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGc                                                                         >  1:1669389/1‑93 (MQ=255)
                 tcATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCg                                                                       <  1:5692677/92‑1 (MQ=255)
                   aTACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGa                                                                     >  2:9403640/1‑93 (MQ=255)
                   aTACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGa                                                                     <  1:7391080/93‑1 (MQ=255)
                   aTACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGa                                                                     <  1:3844987/93‑1 (MQ=255)
                   aTACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGa                                                                     >  1:325299/1‑93 (MQ=255)
                    tACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGAt                                                                    <  2:6808465/93‑1 (MQ=255)
                    tACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGAt                                                                    >  1:1256068/1‑93 (MQ=255)
                     aCAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGAtt                                                                   >  1:8377997/1‑93 (MQ=255)
                         cGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTcccc                                                               <  2:2956718/93‑1 (MQ=255)
                          gAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCg                                                              >  1:8733147/1‑93 (MQ=255)
                          gAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCg                                                              <  1:3108188/93‑1 (MQ=255)
                             cGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACa                                                           >  2:2417360/1‑93 (MQ=255)
                              ggTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACaa                                                          >  2:7455162/1‑93 (MQ=255)
                              ggTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACaa                                                          <  2:7879738/93‑1 (MQ=255)
                               gTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGTATTCCCCGACAAt                                                         >  1:9324104/1‑93 (MQ=255)
                               gTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACaag                                                         >  1:7726513/1‑92 (MQ=255)
                               gTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAAt                                                         >  2:9915744/1‑93 (MQ=255)
                               gTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAAt                                                         >  2:124113/1‑93 (MQ=255)
                                ttCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATg                                                        >  1:6371670/1‑93 (MQ=255)
                                   tcGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGccc                                                     >  1:226703/1‑93 (MQ=255)
                                     gTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAg                                                   >  2:665141/1‑93 (MQ=255)
                                     gTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAg                                                   <  1:5437776/93‑1 (MQ=255)
                                      tAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGc                                                  <  2:4437085/93‑1 (MQ=255)
                                       aGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCa                                                 <  2:3176541/93‑1 (MQ=255)
                                        gTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCAc                                                >  2:6137695/1‑93 (MQ=255)
                                        gTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCAc                                                >  1:1905943/1‑93 (MQ=255)
                                         tGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACt                                               <  2:6625700/93‑1 (MQ=255)
                                         tGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACt                                               >  2:415058/1‑93 (MQ=255)
                                           aaTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTc                                             <  2:4251409/93‑1 (MQ=255)
                                           aaTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTc                                             >  2:257933/1‑93 (MQ=255)
                                           aaTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTc                                             >  1:2405478/1‑93 (MQ=255)
                                            aTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTcc                                            >  1:5459695/1‑93 (MQ=255)
                                            aTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTcc                                            <  2:1707268/93‑1 (MQ=255)
                                            aTTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTcc                                            <  2:2810550/93‑1 (MQ=255)
                                             tttCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTccc                                           <  1:8225537/93‑1 (MQ=255)
                                              ttCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTcccc                                          >  1:5794947/1‑93 (MQ=255)
                                               tCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACACTGCCCAGCCCTTccccc                                         >  1:9050559/1‑93 (MQ=255)
                                               tCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACACTGCCCAGCACTTccccc                                         >  1:7251768/1‑93 (MQ=255)
                                                  ccccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCgc                                       >  1:8586710/1‑92 (MQ=255)
                                                    ccTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGcc                                     >  2:4493456/1‑93 (MQ=255)
                                                     ctgctgCGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCa                                    >  1:9629211/1‑93 (MQ=255)
                                                      tgctgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCAt                                   >  1:6836396/1‑93 (MQ=255)
                                                       gctgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATa                                  >  1:7696240/1‑93 (MQ=255)
                                                         tgcGGCGTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATaaa                                >  2:9528548/1‑93 (MQ=255)
                                                           cggcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATTCCCAGCACTTCCCCCGGCCATAAATc                              <  1:550046/93‑1 (MQ=255)
                                                            ggcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCa                             >  2:2672693/1‑93 (MQ=255)
                                                            ggcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCa                             >  2:1108327/1‑93 (MQ=255)
                                                            ggcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCa                             >  2:3790922/1‑93 (MQ=255)
                                                            ggcgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCa                             >  2:8020842/1‑93 (MQ=255)
                                                              cgTCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCaaa                           <  2:3885015/93‑1 (MQ=255)
                                                                tCAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAag                         <  2:5107072/93‑1 (MQ=255)
                                                                 cAGGCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAaga                        <  2:3692067/93‑1 (MQ=255)
                                                                   ggCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAgaga                      <  2:178683/93‑1 (MQ=255)
                                                                   ggCGCGCGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAgaga                      <  2:7417650/93‑1 (MQ=255)
                                                                    gcgcgcGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGAc                     <  1:1087669/93‑1 (MQ=255)
                                                                    gcgcgcGGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGAc                     <  1:4772648/93‑1 (MQ=255)
                                                                     cgcgcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACa                    <  2:5399263/93‑1 (MQ=255)
                                                                      gcgcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACAt                   >  2:5246526/1‑93 (MQ=255)
                                                                        gcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCg                 >  1:183828/1‑93 (MQ=255)
                                                                        gcgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCg                 >  1:3549129/1‑93 (MQ=255)
                                                                         cgGAGATCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGc                <  1:156589/93‑1 (MQ=255)
                                                                             gaTCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATTCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTaaa            >  2:577428/1‑93 (MQ=255)
                                                                             gaTCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTaaa            >  2:4357780/1‑93 (MQ=255)
                                                                              aTCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAg           <  1:4577318/93‑1 (MQ=255)
                                                                               tCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGc          <  1:1432640/93‑1 (MQ=255)
                                                                               tCGACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGc          >  2:4312461/1‑93 (MQ=255)
                                                                                 gACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCt        <  2:48471/93‑1 (MQ=255)
                                                                                 gACTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCt        <  1:1214317/93‑1 (MQ=255)
                                                                                  aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       >  2:2096687/1‑93 (MQ=255)
                                                                                  aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       >  1:10016271/1‑93 (MQ=255)
                                                                                  aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       >  1:9403210/1‑93 (MQ=255)
                                                                                  aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       <  1:7931549/93‑1 (MQ=255)
                                                                                  aCTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCtt       >  2:9041939/1‑93 (MQ=255)
                                                                                   cTTCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAATGAGACATCGCTAAAGCCttt      <  2:9682717/93‑1 (MQ=255)
                                                                                    ttCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTg     <  1:6514845/93‑1 (MQ=255)
                                                                                     tCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGc    <  2:5867246/93‑1 (MQ=255)
                                                                                     tCAGCAGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGc    <  2:852721/93‑1 (MQ=255)
                                                                                      cagcagCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGcc   >  1:9067532/1‑93 (MQ=255)
                                                                                      cAGCGGCGTGTTCTTCCCGGTCGCGGATACCCCGACAATGCCAAGCATTGCCACCGGCCATAAATCAAAAGAGACATCGAGAAAGCCTTTGcc   <  2:6276324/93‑1 (MQ=255)
                                                                                       agcagcGTGGTCTTCCCGGAGCCGGATCCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGccc  <  1:1365450/93‑1 (MQ=255)
                                                                                       agcagcGTGGTCTTCCCGGAGCCGGATCCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGccc  <  1:29955/93‑1 (MQ=255)
                                                                                          |                                                                                         
GGTCGGCCTCGCTCATTGCATACAGCGAACGGTTCTCGTAGTGAATTTCCCCCTGCTGCGGCGTCAGGCGCGCGGAGATCGACTTCAGCTGCGTGGTCTTCCCGGAGCCGGATTCCCCGACAATGCCCAGCACTTCCCCCGGCCATAAATCAAAAGAGACATCGCTAAAGCCTTTGCCC  >  NC_000913/4316114‑4316292

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