Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 454,856 G→A 17.4% intergenic (+67/‑277) bolA → / → tig stationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC/peptidyl‑prolyl cis/trans isomerase (trigger factor)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913454,8560GA17.4% 213.2 / 33.4 92intergenic (+67/‑277)bolA/tigstationary‑phase morphogene, transcriptional repressor for mreB; also regulator for dacA, dacC, and ampC/peptidyl‑prolyl cis/trans isomerase (trigger factor)
Reads supporting (aligned to +/- strand):  ref base G (42/34);  new base A (4/12);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 5.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

CCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCG  >  NC_000913/454757‑454955
                                                                                                   |                                                                                                   
cctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCg                                                                                                     >  1:1991462/1‑100 (MQ=255)
 ctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCgg                                                                                                    >  2:1718985/1‑100 (MQ=255)
  tccCTGTGGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCggg                                                                                                   >  1:1371761/1‑100 (MQ=255)
  tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCggg                                                                                                   <  2:661653/100‑1 (MQ=255)
   cccTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGc                                                                                                  <  2:1773834/100‑1 (MQ=255)
     cTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCg                                                                                                <  1:1689182/100‑1 (MQ=255)
     cTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCg                                                                                                >  1:1618715/1‑100 (MQ=255)
        tCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgttt                                                                                             <  2:1982805/100‑1 (MQ=255)
          gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCTgttttg                                                                                           >  1:1820034/1‑100 (MQ=255)
          gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttg                                                                                           >  1:2027864/1‑100 (MQ=255)
            ggAGCAGGAAGCATCTCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                         <  1:2561694/100‑1 (MQ=255)
            ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTATAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                         <  2:759620/100‑1 (MQ=255)
            ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                         >  1:1030197/1‑100 (MQ=255)
            ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                         >  1:360593/1‑100 (MQ=255)
            ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                         <  1:2232534/100‑4 (MQ=255)
            ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                         <  2:2057021/100‑4 (MQ=255)
            ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttc                                                                                          <  1:1659901/99‑3 (MQ=255)
             gggCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTGATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                        <  1:1244108/98‑5 (MQ=38)
             gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                        <  1:2362178/100‑5 (MQ=39)
             gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                        <  1:1408183/100‑5 (MQ=39)
              aGCGGGAAGCATCGCGTAAAAACGTATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcata                                                                                       <  2:946519/100‑6 (MQ=38)
              aGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                       <  1:2122000/100‑1 (MQ=255)
              aGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcata                                                                                       <  1:1425491/100‑6 (MQ=39)
               gCAGTAAGCATCTCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGTCCgtttttcat                                                                                        <  2:2244221/98‑5 (MQ=255)
               gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                       >  2:1673358/1‑99 (MQ=255)
               gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                      >  1:2085251/1‑100 (MQ=255)
               gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                      <  2:1732902/100‑1 (MQ=255)
                cAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcataaa                                                                                     <  2:1341835/100‑8 (MQ=38)
                  ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTGCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                                   <  2:1671823/100‑1 (MQ=255)
                     aGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGa                                                                                <  1:646171/100‑1 (MQ=255)
                       cATCGCGTAAAAACGCATTTTCAACTGTCGGCGCTTTTCCAGTAGGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGCtt                                                                              >  2:1589011/1‑100 (MQ=255)
                         tctcgTAATATCGCATTAGCGACTTTCGGCGCTTCTCCACTATGTGCCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                         <  2:2546414/84‑4 (MQ=16)
                         tCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAtttt                                                                            >  2:1286979/1‑100 (MQ=255)
                         tCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAtttt                                                                            <  2:2027853/100‑1 (MQ=255)
                          cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                           >  2:1671456/1‑100 (MQ=255)
                             gTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgcgc                                                                        >  2:148744/1‑100 (MQ=255)
                              tAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCt                                                                       >  1:859540/1‑100 (MQ=255)
                                aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttt                                                                     >  1:495523/1‑100 (MQ=255)
                                aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttt                                                                     >  2:2587791/1‑100 (MQ=255)
                                 aaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCtttt                                                                    <  1:363594/100‑1 (MQ=255)
                                  aaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttttt                                                                   >  2:1866784/1‑100 (MQ=255)
                                   acctaTTTGCTACTGTCGTCGCTTCCCCAGTATGTTGCTAAATATTTTATGAAACACGGCCTGCGGGCCGTTTTGTTGTGTCTGGATTTTGCGCttttcg                                                                  <  2:1912833/96‑3 (MQ=255)
                                    cGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTATGGATTTTGCGCTTTTTGc                                                                 >  2:694805/1‑100 (MQ=255)
                                      cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                               >  2:383109/1‑100 (MQ=255)
                                      cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                               <  2:1913454/100‑1 (MQ=255)
                                          tGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGc                                                            >  2:927034/1‑99 (MQ=255)
                                           gCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAt                                                          <  2:2574747/100‑1 (MQ=255)
                                            cAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAtt                                                         >  2:243213/1‑100 (MQ=255)
                                            cAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAtt                                                         >  1:1821025/1‑100 (MQ=255)
                                              aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                       <  2:833511/100‑1 (MQ=255)
                                              aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                       <  1:1532575/100‑1 (MQ=255)
                                               cTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAg                                                      <  2:218767/100‑1 (MQ=255)
                                                 gTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGAc                                                    <  1:850784/100‑1 (MQ=255)
                                                      cgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaaaa                                               <  1:23813/100‑1 (MQ=255)
                                                      cgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaaaa                                               >  2:2042912/1‑100 (MQ=255)
                                                          tttCCAGTATGTTGCTAAAGATGTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTg                                            <  2:2140841/99‑1 (MQ=255)
                                                            tCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGccc                                         >  1:691993/1‑100 (MQ=255)
                                                            tCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGccc                                         >  2:2534396/1‑100 (MQ=255)
                                                             ccAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGGTTTGTTTTGACTGGATTTTGAGCGTTTTCCCCAGCATTCAGTCGACATTTACCCg                                        >  2:589466/1‑100 (MQ=255)
                                                               aGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCggg                                      >  2:124494/1‑100 (MQ=255)
                                                               aGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCACTTTTTGCCCAGCATTCAGACGAAAATTGCCCggg                                      <  2:1544800/100‑1 (MQ=255)
                                                                 tATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGaa                                    <  1:836098/100‑1 (MQ=255)
                                                                   tGTTACTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAtt                                  >  2:118829/1‑100 (MQ=255)
                                                                    gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                                 <  1:1965456/100‑1 (MQ=255)
                                                                    gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                                 <  1:1718994/100‑1 (MQ=255)
                                                                    gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                                 <  1:249512/100‑1 (MQ=255)
                                                                      tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgtg                               <  1:1486481/100‑1 (MQ=255)
                                                                      tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATGTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgtg                               <  2:1598200/100‑1 (MQ=255)
                                                                       gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                              >  1:1333981/1‑100 (MQ=255)
                                                                       gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                              >  1:9143/1‑100 (MQ=255)
                                                                       gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                              >  1:304188/1‑100 (MQ=255)
                                                                       gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                              >  2:2117681/1‑100 (MQ=255)
                                                                        cTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaa                             >  1:391920/1‑100 (MQ=255)
                                                                        cTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaa                             >  1:644867/1‑100 (MQ=255)
                                                                        cTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                              <  1:124495/99‑1 (MQ=255)
                                                                         tAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTAAGACGAAAATTGCCCGGGAATTGTGaaa                            >  1:830066/1‑100 (MQ=255)
                                                                          aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaa                           >  1:2287249/1‑100 (MQ=255)
                                                                          aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaa                           >  1:397536/1‑100 (MQ=255)
                                                                          aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTATGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgtg                               >  2:2375320/1‑96 (MQ=255)
                                                                            aGATTTTATGAAAAACGGCCTGCGGGCCTTTTTGTTTTGTCTGGATTTTGCGCTTTTTTCCCAGCATTCAGACGAAAATTGCCCGGTAATTGTGaaaaaa                         <  2:2016489/100‑1 (MQ=255)
                                                                             gATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTTTTTATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCTGGGGTATTGTGAAAAAAt                        >  2:1457568/1‑100 (MQ=255)
                                                                               gacaaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAc                      >  1:1012302/9‑100 (MQ=38)
                                                                                tttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcg                     <  1:1269163/100‑1 (MQ=255)
                                                                                 caaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgc                    >  2:1340034/7‑100 (MQ=38)
                                                                                    tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAc                      >  2:826085/1‑95 (MQ=255)
                                                                                    tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc                 >  1:332231/1‑100 (MQ=255)
                                                                                    tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc                 >  1:1518871/1‑100 (MQ=255)
                                                                                    tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCAGGGAATTGTGAAACAATACGCGAc                 >  1:681176/1‑100 (MQ=255)
                                                                                    aacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc                 >  2:1188057/4‑100 (MQ=255)
                                                                                    aacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc                 >  1:174907/4‑100 (MQ=255)
                                                                                      aaaaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAg               <  2:707712/100‑1 (MQ=255)
                                                                                         aaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGTCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgc            <  1:442924/100‑1 (MQ=255)
                                                                                              ccTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAAt       <  1:1117192/100‑1 (MQ=255)
                                                                                                tGCGGCCCGTTTTGTTTTGTCTGGACTTTGCGCTTTTTGCCCAGCATTCAGACGCACATTGCCCGGGAATTGTGAAAAAATACCCGACAGCGCGCAATaa     >  2:1165335/1‑100 (MQ=255)
                                                                                                  cGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAcc   <  2:1297665/100‑1 (MQ=255)
                                                                                                  cGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAcc   <  1:1048738/100‑1 (MQ=255)
                                                                                                   gggCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCg  <  1:1167035/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCG  >  NC_000913/454757‑454955

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: