Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,318,840 A→T 77.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,8400AT77.9% 317.2 / 91.4 158L256Q (CTG→CAG) phnJcarbon‑phosphorus lyase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (19/16);  new base T (58/65);  total (77/81)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.96e-01

Only 100 of 160 total aligned reads displayed.

GACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAATCGAGACT  >  NC_000913/4318741‑4318937
                                                                                                   |                                                                                                 
  ccGAAAGTAACGGTTGATTCATTGGTTTTTTTTCTTGTTCTTTTGGCGGCAATAATCGGTATCGTTGCAGACAAACATGCGGTTTCCGGCGTCATCCTGc                                                                                                 >  2:3664337/1‑100 (MQ=255)
  ccGAAAGTAACGGTTGATTCATTGGTTTTTTGTCTCTTTTTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGc                                                                                                 >  2:1790791/1‑100 (MQ=255)
  ccGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGc                                                                                                 <  2:562561/100‑1 (MQ=255)
   cGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGca                                                                                                <  2:4061057/100‑1 (MQ=255)
    gAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGcac                                                                                               <  2:1891973/100‑1 (MQ=255)
     aaaGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGcacc                                                                                              <  1:1560619/100‑1 (MQ=255)
     aaaGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGcacc                                                                                              <  2:3290034/100‑1 (MQ=255)
     aaaGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGcacc                                                                                              >  1:2196404/1‑100 (MQ=255)
      aaGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGcacca                                                                                             >  1:1633288/1‑100 (MQ=255)
       aGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGcaccac                                                                                            >  1:280107/1‑100 (MQ=255)
           aCGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTtca                                                                                        <  1:3773763/100‑1 (MQ=255)
            tggTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTtcat                                                                                       >  1:3552518/2‑100 (MQ=255)
            cGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTtcat                                                                                       <  1:3237432/100‑1 (MQ=255)
                 gATTCATTGGTTTTTTTTTTTTTTTTGTTGGCGGCAATAACCGGTATCAGAGCAGTCAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGa                                                                                  >  2:2653698/1‑100 (MQ=255)
                 gATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGa                                                                                  >  1:3461558/1‑100 (MQ=255)
                    tCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAg                                                                               <  2:1126307/100‑1 (MQ=255)
                     cATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGc                                                                              >  1:1969621/1‑100 (MQ=255)
                          gTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTgg                                                                         <  2:1334904/100‑1 (MQ=255)
                           ttttttGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTggg                                                                        <  2:2156154/100‑1 (MQ=255)
                             ttttGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTc                                                                      >  1:3021733/1‑100 (MQ=255)
                             ttttGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTc                                                                      <  1:1775445/100‑1 (MQ=255)
                             ttttGCCTCGCGCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTc                                                                      <  1:1683753/100‑1 (MQ=255)
                                tGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGAt                                                                   <  2:4137483/100‑1 (MQ=255)
                                 gCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATc                                                                  <  1:2341324/100‑1 (MQ=255)
                                  ccTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATcc                                                                 >  1:661482/1‑100 (MQ=255)
                                   cTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCCCTTCATCAAGATAGCTGTGGGTCGATCCg                                                                >  2:3269921/1‑100 (MQ=255)
                                    tCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGc                                                               >  1:3243621/1‑100 (MQ=255)
                                     cgcgCTGTTGGCGGCAATAAGCGGGATCGGAGCAGCGAAGCATGCGGTTGGCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCa                                                              <  2:3723195/100‑1 (MQ=255)
                                     cGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCa                                                              <  1:2445982/100‑1 (MQ=255)
                                       ctctGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGa                                                            >  1:1181978/1‑100 (MQ=255)
                                           gTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGgc                                                        >  2:3851621/1‑100 (MQ=255)
                                           gTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGgc                                                        >  1:2049582/1‑100 (MQ=255)
                                           gTTGGCGGCAATAATCGCTATCGGAGCACACAACCATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATgg                                                         <  2:1638736/99‑1 (MQ=255)
                                              ggcggcAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCa                                                     >  2:648109/1‑100 (MQ=255)
                                              ggcggcAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCa                                                     <  1:3851981/100‑1 (MQ=255)
                                              ggcggcAATAACCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCCTCAAGATAGCTGCGGATCGACCCACAgctatctt                                                      <  2:1446730/99‑8 (MQ=37)
                                               gcggcAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCAt                                                    >  1:3656030/1‑100 (MQ=255)
                                               gcggcAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCAGCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCAt                                                    <  2:3764724/100‑1 (MQ=255)
                                                cggcAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATg                                                   >  1:3832915/1‑100 (MQ=255)
                                                  gcAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGc                                                 <  2:2094391/100‑1 (MQ=255)
                                                   cAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCt                                                <  1:25214/100‑1 (MQ=255)
                                                   cAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGACGCATGGCt                                                >  1:1421316/1‑100 (MQ=255)
                                                   cAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATAAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCt                                                >  1:568672/1‑100 (MQ=255)
                                                    aataatCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTc                                               >  1:487601/1‑100 (MQ=255)
                                                     ataatCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCa                                              >  2:469552/1‑100 (MQ=255)
                                                      taatCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCAt                                             <  1:1175809/100‑1 (MQ=255)
                                                         tCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATccc                                          <  1:417814/100‑1 (MQ=255)
                                                         tCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGCCATCCTGCACCACTGCATCAAGATCGCTGTGGCACAATCCGCAGACGCCGCATGGCCCATcc                                           >  2:3376358/1‑99 (MQ=255)
                                                          cGGTATCGAAGCAAACAAACATGCGGCTTCCGGCGTCATCCTGCACCACTTCATCAAGATCGCTGTGGGTCCATCCGCAAATGGCGCATGGCTCAcccc                                          <  2:1647929/99‑1 (MQ=255)
                                                             tATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActg                                      >  1:3776798/1‑100 (MQ=255)
                                                              aTCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActgc                                     >  2:2886632/1‑100 (MQ=255)
                                                              aTCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActgc                                     <  2:2082090/100‑1 (MQ=255)
                                                               tCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActgct                                    <  2:239618/100‑1 (MQ=255)
                                                               tCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActgct                                    <  2:673576/100‑1 (MQ=255)
                                                               tCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActgct                                    >  1:3251772/1‑100 (MQ=255)
                                                               tCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCActgct                                    >  1:197272/1‑100 (MQ=255)
                                                                cGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGGTCCGCAGATGGCGCATGGCTCATCCCActgctg                                   <  2:3978089/100‑1 (MQ=255)
                                                                 ggAGCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtga                                  >  1:179551/1‑100 (MQ=255)
                                                                   aGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCCATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaac                                <  2:2223032/100‑1 (MQ=255)
                                                                    gCAGACAAACCTGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATCGCTGTGGGTCGATCCGCCGATCGCGCATGCGTCCTCCCACTGCtga                                  >  2:3647755/1‑97 (MQ=255)
                                                                    gCAGACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacg                               <  1:1855065/100‑1 (MQ=255)
                                                                     cAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGTtgaacgg                              <  1:2192634/100‑1 (MQ=255)
                                                                     cAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacgg                              <  2:1109797/100‑1 (MQ=255)
                                                                       gACAAACATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtg                            >  2:1860239/1‑100 (MQ=255)
                                                                       gACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtg                            <  2:1286525/100‑1 (MQ=255)
                                                                          aaaCATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtgaac                         <  1:880392/100‑1 (MQ=255)
                                                                           aaCATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtgaacg                        <  1:873216/100‑1 (MQ=255)
                                                                           aaCATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtgaacg                        >  1:797816/1‑100 (MQ=255)
                                                                            aCATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCtgaacggtgaacgg                       >  2:1545028/1‑100 (MQ=255)
                                                                             cATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACggg                      >  1:2140354/1‑100 (MQ=255)
                                                                             cATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACggg                      >  1:3329502/1‑100 (MQ=255)
                                                                             cATGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACggg                      <  1:458744/100‑1 (MQ=255)
                                                                              aTGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGt                     >  1:2903239/1‑100 (MQ=255)
                                                                               tGCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCTCTGCTGAACGGTGAACGGGTg                    <  1:992545/100‑1 (MQ=255)
                                                                               tGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTg                    >  1:2913991/1‑100 (MQ=255)
                                                                                gCGGTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTAGCTCCGCAGATGGCGCATGGCTCCTCCCACTGCtgaacggtgcacgg                       >  2:1542777/1‑96 (MQ=255)
                                                                                 cGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAt                  >  1:3449183/1‑100 (MQ=255)
                                                                                   gTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGACGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtcg                >  2:3328219/1‑100 (MQ=255)
                                                                                   gTTTCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGCCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtcg                >  1:736271/1‑100 (MQ=255)
                                                                                   gTTTCCGGCGTCATCCTGCACCACTTCAGCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtcg                <  1:4087960/100‑1 (MQ=255)
                                                                                      tCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCGTCCCACTGCTGAACGGTGAACGGGTGAtcgtcg             <  1:3660795/100‑1 (MQ=255)
                                                                                      tCCGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtcgtcg             <  2:1552022/100‑1 (MQ=255)
                                                                                       ccGGCGTCATCCTGCGCCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGa            <  2:3834388/100‑1 (MQ=255)
                                                                                       ccGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGa            >  2:895396/1‑100 (MQ=255)
                                                                                       ccGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGa            >  1:3738929/1‑100 (MQ=255)
                                                                                       ccGGCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGa            >  2:3129/1‑100 (MQ=255)
                                                                                       ccGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGa            >  2:232634/1‑100 (MQ=255)
                                                                                       ccGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGa            >  1:2060993/1‑100 (MQ=255)
                                                                                         ggCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGACGATCCGCAGATGGCGCATTGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGaaa          <  2:2274487/100‑1 (MQ=255)
                                                                                         ggCGTCATCCAGCACCACTTCATCAAGATAGGTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGAtcgtcgca           >  1:917478/1‑97 (MQ=255)
                                                                                          gCGTCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAAt         <  1:77023/100‑1 (MQ=255)
                                                                                             tCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAATCga      >  2:1203853/1‑100 (MQ=255)
                                                                                             tCATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAATCg       <  2:1425923/99‑1 (MQ=255)
                                                                                              cATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCTACTGCTGAACGGTGAACGGGTGATCGTCGAAATCgag     >  1:1346520/1‑100 (MQ=255)
                                                                                              cATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAATCgag     <  1:2718517/100‑1 (MQ=255)
                                                                                              cATCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAATCgag     <  2:2794337/100‑1 (MQ=255)
                                                                                               aTCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAATCgaga    <  1:915297/100‑1 (MQ=255)
                                                                                               aTCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGGGAACGGGGGATCGTCGAAACCgaga    >  2:3138332/1‑100 (MQ=255)
                                                                                                tCCTGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAATCGAGAc   >  1:847203/1‑100 (MQ=255)
                                                                                                 ccAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAATCGAGACt  >  1:722846/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
GACCGAAAGTAACGGTTGATTCATTGGTTTTTTGCCTCGCTCTGTTGGCGGCAATAATCGGTATCGGAGCAGACAAACATGCGGTTTCCGGCGTCATCCAGCACCACTTCATCAAGATAGCTGTGGGTCGATCCGCAGATGGCGCATGGCTCATCCCACTGCTGAACGGTGAACGGGTGATCGTCGAAATCGAGACT  >  NC_000913/4318741‑4318937

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: