Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 454,853 T→C 17.0% intergenic (+64/‑280) 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,8530TC17.0% 204.5 / 35.7 94intergenic (+64/‑280)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 T (45/33);  new base C (4/12);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 2.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

TCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAA  >  NC_000913/454754‑454952
                                                                                                   |                                                                                                   
tctcCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCt                                                                                                     >  2:653183/1‑100 (MQ=255)
  tcctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGc                                                                                                   >  1:218769/1‑100 (MQ=255)
   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:360593/1‑100 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                      >  1:1030197/1‑100 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                      <  2:2057021/100‑4 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                      <  1:2232534/100‑4 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttc                                                                                       <  1:1659901/99‑3 (MQ=255)
                gggCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTGATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                     <  1:1244108/98‑5 (MQ=38)
                gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                     <  1:1408183/100‑5 (MQ=39)
                gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                     <  1:2362178/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                                                                                   <  2:1732902/100‑1 (MQ=255)
                  gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                   >  1:2085251/1‑100 (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:1913454/100‑1 (MQ=255)
                                         cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                            >  2:383109/1‑100 (MQ=255)
                                             tGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGc                                                         >  2:927034/1‑99 (MQ=255)
                                              gCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAt                                                       <  2:2574747/100‑1 (MQ=255)
                                               cAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAtt                                                      >  1:1821025/1‑100 (MQ=255)
                                               cAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAtt                                                      >  2:243213/1‑100 (MQ=255)
                                                 aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                    <  1:1532575/100‑1 (MQ=255)
                                                 aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                    <  2:833511/100‑1 (MQ=255)
                                                  cTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAg                                                   <  2:218767/100‑1 (MQ=255)
                                                    gTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGAc                                                 <  1:850784/100‑1 (MQ=255)
                                                         cgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaaaa                                            >  2:2042912/1‑100 (MQ=255)
                                                         cgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaaaa                                            <  1:23813/100‑1 (MQ=255)
                                                             tttCCAGTATGTTGCTAAAGATGTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTg                                         <  2:2140841/99‑1 (MQ=255)
                                                               tCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGccc                                      >  2:2534396/1‑100 (MQ=255)
                                                               tCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGccc                                      >  1:691993/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:1718994/100‑1 (MQ=255)
                                                                       gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                              <  1:1965456/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                           >  2:2117681/1‑100 (MQ=255)
                                                                          gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                           >  1:304188/1‑100 (MQ=255)
                                                                          gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                           >  1:9143/1‑100 (MQ=255)
                                                                           cTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaa                          >  1:644867/1‑100 (MQ=255)
                                                                           cTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaa                          >  1:391920/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:1518871/1‑100 (MQ=255)
                                                                                       tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc              >  1:332231/1‑100 (MQ=255)
                                                                                       tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCAGGGAATTGTGAAACAATACGCGAc              >  1:681176/1‑100 (MQ=255)
                                                                                       aacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc              >  1:174907/4‑100 (MQ=255)
                                                                                       aacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc              >  2:1188057/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)
                                                                                                   |                                                                                                   
TCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAA  >  NC_000913/454754‑454952

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: