Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,318,840 A→T 48.9% L256Q (CTG→CAG)  phnJ ← carbon‑phosphorus lyase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,318,8400AT48.9% 18.9 / 268.3 182L256Q (CTG→CAG) phnJcarbon‑phosphorus lyase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (46/47);  new base T (51/38);  total (97/85)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

GACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAATCGAGAC  >  NC_000913/4318741‑4318936
                                                                                                   |                                                                                                
gACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCa                                                                                                  <  2:2096192/100‑1 (MQ=255)
gACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCa                                                                                                  >  1:2342453/1‑100 (MQ=255)
  ccGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGc                                                                                                >  2:2857061/1‑100 (MQ=255)
   cGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGca                                                                                               >  2:1478149/1‑100 (MQ=255)
   cGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGca                                                                                               <  1:1998850/100‑1 (MQ=255)
    gAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGcac                                                                                              >  2:3675507/1‑100 (MQ=255)
     aaaGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGcacc                                                                                             <  1:3030621/100‑1 (MQ=255)
           aCGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTtca                                                                                       <  2:3726750/100‑1 (MQ=255)
           aCGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGGATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTtca                                                                                       <  1:3536920/100‑1 (MQ=255)
             ngTTGATTCATTGGTTTTTTTCTTTTTTTTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTtcatc                                                                                     >  1:3948198/2‑100 (MQ=255)
             ggTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGTAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATTCAGCACCACTtcatc                                                                                     >  2:4204042/1‑100 (MQ=255)
              gTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTtcatca                                                                                    <  2:3934431/100‑1 (MQ=255)
               ttGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCaa                                                                                   <  2:2674737/100‑1 (MQ=255)
                tGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAg                                                                                  <  2:4015636/100‑1 (MQ=255)
                 gATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGa                                                                                 >  2:2463595/1‑100 (MQ=255)
                    tCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACCAACCTGCGGCCTCCCGCGTCATCCTGCCCCACTTCATCAAGCTa                                                                               >  1:406664/1‑99 (MQ=255)
                     cATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGc                                                                             >  1:1761936/1‑100 (MQ=255)
                     cATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGc                                                                             >  2:655558/1‑100 (MQ=255)
                      aTTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCt                                                                            >  1:710848/1‑100 (MQ=255)
                      aTTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCt                                                                            >  2:4166206/1‑100 (MQ=255)
                       ttGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGTCGTCATCCTGCACCACTTCATCACGATAGc                                                                             >  1:2582261/1‑98 (MQ=255)
                         ggTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCtgtg                                                                         >  2:3638967/1‑100 (MQ=255)
                         ggTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCtgtg                                                                         <  2:2346827/100‑1 (MQ=255)
                         ggTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCtgtg                                                                         <  1:1528647/100‑1 (MQ=255)
                         ggTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCtgtg                                                                         <  1:986725/100‑1 (MQ=255)
                         gattttttGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCtgtg                                                                         >  2:2884628/3‑100 (MQ=255)
                          gTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTgg                                                                        >  2:4023410/1‑100 (MQ=255)
                          gTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTgg                                                                        <  1:3677350/100‑1 (MQ=255)
                           ttttttGCCTTGCTCTGTTGGCGGCAAAAATCGGTATAGGACCAGGCAACCATGCGGTTTCCGCCGTCATCCTGCACCACTTCATCAAGATAGCTGAggg                                                                       >  1:3227135/1‑100 (MQ=255)
                           ttttttGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTggg                                                                       >  2:202089/1‑100 (MQ=255)
                           ttttttGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTggg                                                                       <  2:3370190/100‑1 (MQ=255)
                            tttttGCCTTGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGt                                                                      >  2:3848737/1‑100 (MQ=255)
                            tttttGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGt                                                                      >  1:1907415/1‑100 (MQ=255)
                               ttGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGa                                                                   <  2:1294530/100‑1 (MQ=255)
                               ttGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGa                                                                   <  2:2956160/100‑1 (MQ=255)
                                tGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGAt                                                                  >  1:598554/1‑100 (MQ=255)
                                  ccTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGAccc                                                                >  1:3322916/1‑100 (MQ=255)
                                  ccTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATcc                                                                >  1:2012257/1‑100 (MQ=255)
                                  ccTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATcc                                                                <  1:2778796/100‑1 (MQ=255)
                                     cGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCa                                                             <  1:4024982/100‑1 (MQ=255)
                                      gCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAg                                                            <  2:2719513/100‑1 (MQ=255)
                                       ctctGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGa                                                           >  1:2447347/1‑100 (MQ=255)
                                         ctGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATg                                                         <  1:3788081/100‑1 (MQ=255)
                                          tGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATgg                                                        <  2:696219/100‑1 (MQ=255)
                                             tGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATACAGATAGCTGTGGCTCGATCCGCAgccg                                                         >  1:1537847/1‑93 (MQ=255)
                                              ggcggcAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCa                                                    >  2:2350348/1‑100 (MQ=255)
                                               gcggcAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCAt                                                   >  1:2392972/1‑100 (MQ=255)
                                               gcggcAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCa                                                    >  1:845546/1‑99 (MQ=255)
                                               gcggcAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCAt                                                   >  2:1896052/1‑100 (MQ=255)
                                                cggcAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATg                                                  >  1:653019/1‑100 (MQ=255)
                                                cggcAATAATCGGTATCGGAGCAGACAAACATGCGGTGTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATg                                                  >  1:2288583/1‑100 (MQ=255)
                                                  gcAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGc                                                <  2:1161061/100‑1 (MQ=255)
                                                   cAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCt                                               <  2:3535130/100‑1 (MQ=255)
                                                   cAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCt                                               <  2:1724446/100‑1 (MQ=255)
                                                   cAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCt                                               <  2:3071341/100‑1 (MQ=255)
                                                       aatCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATc                                           >  1:3799594/1‑100 (MQ=255)
                                                           ggTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCAc                                       <  2:1329537/100‑1 (MQ=255)
                                                           ggTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCAc                                       >  2:2641232/1‑100 (MQ=255)
                                                           ggTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGACTCATCCCAc                                       >  1:280528/1‑100 (MQ=255)
                                                            gTATCGGAGCAGACAAACATGCGGTTTCCGTCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCAc                                       <  2:3902747/99‑1 (MQ=255)
                                                              aTCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActgc                                    <  1:3477997/100‑1 (MQ=255)
                                                              aTCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActgc                                    <  1:799083/100‑1 (MQ=255)
                                                              aTCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActgc                                    >  1:3872287/1‑100 (MQ=255)
                                                               tCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActgct                                   <  2:1622423/100‑1 (MQ=255)
                                                                cGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActgc                                    >  2:3381918/1‑98 (MQ=255)
                                                                 ggAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCCCCACTTCAGCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtga                                 <  1:291559/100‑1 (MQ=255)
                                                                 ggAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtga                                 <  2:1251364/100‑1 (MQ=255)
                                                                  gAGTAGACGAACCGGCGGGTTCCGGCGTCATGCTCCACGACTTCATCAAGATAGCTGGGGGTCGATCCGCAGATGGCGCATGGCTCATccc                                         <  1:1690450/91‑1 (MQ=255)
                                                                  gAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActgctg                                  >  1:4259429/1‑98 (MQ=255)
                                                                    gCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacg                              >  2:1664107/1‑100 (MQ=255)
                                                                     cAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacgg                             <  1:1986058/100‑1 (MQ=255)
                                                                     cAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacgg                             <  2:2263929/100‑1 (MQ=255)
                                                                     cAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacgg                             <  2:772499/100‑1 (MQ=255)
                                                                      aGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggt                            >  2:831350/1‑100 (MQ=255)
                                                                      aGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggt                            >  1:1343459/1‑100 (MQ=255)
                                                                      aGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggt                            <  1:220645/100‑1 (MQ=255)
                                                                      aGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCACATGGCGCATGGCTCATCCCACTGCtgaacggt                            <  2:4053113/100‑1 (MQ=255)
                                                                       gACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtg                           >  2:3803283/1‑100 (MQ=255)
                                                                       gACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtg                           >  2:3170128/1‑100 (MQ=255)
                                                                        aCAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtga                          <  1:1816702/100‑1 (MQ=255)
                                                                         cAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtgaa                         >  1:3265760/1‑100 (MQ=255)
                                                                          aaaCATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtgaac                        <  2:4128943/100‑1 (MQ=255)
                                                                           aaCATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtgaacg                       <  2:2984349/100‑1 (MQ=255)
                                                                           aaCATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtgaacg                       <  2:1372091/100‑1 (MQ=255)
                                                                            aCATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtgaacgg                      <  1:4247571/100‑1 (MQ=255)
                                                                               tGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtgaacgg                      >  1:1296611/1‑97 (MQ=255)
                                                                                gCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGa                  <  2:3900252/100‑1 (MQ=255)
                                                                                 cGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAt                 >  1:4169384/1‑100 (MQ=255)
                                                                                   gTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtcg               >  1:3173869/1‑100 (MQ=255)
                                                                                   gTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtcg               >  1:756092/1‑100 (MQ=255)
                                                                                     ttCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtcgtc             >  2:3993093/1‑100 (MQ=255)
                                                                                     ttCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtcgtc             >  2:3500187/1‑100 (MQ=255)
                                                                                     ttCCGGCGTCATCCTGCACCACTTCAGCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtcgtc             <  1:464924/100‑1 (MQ=255)
                                                                                       ccGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGCTGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGa           >  1:1384067/1‑100 (MQ=255)
                                                                                        cGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGa           >  1:1163201/1‑99 (MQ=255)
                                                                                        cGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGaa          >  2:252405/1‑100 (MQ=255)
                                                                                         ggCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGaaa         >  1:101529/1‑100 (MQ=255)
                                                                                          gCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGGGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGa           <  1:1002063/97‑1 (MQ=255)
                                                                                              cATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAATCgag    >  2:2400382/1‑100 (MQ=255)
                                                                                                tCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAATCGAGAc  >  1:1216461/1‑100 (MQ=255)
                                                                                                   |                                                                                                
GACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAATCGAGAC  >  NC_000913/4318741‑4318936

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