Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 454,853 T→C 16.2% 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,8530TC16.2% 392.2 / 62.3 177intergenic (+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 (73/75);  new base C (6/23);  total (79/98)
Fisher's exact test for biased strand distribution p-value = 4.43e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

Only 100 of 179 total aligned reads displayed.

TCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAA  >  NC_000913/454754‑454952
                                                                                                   |                                                                                                   
tctcCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCt                                                                                                     >  1:2781373/1‑100 (MQ=255)
         tGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTAGGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgt                                                                                            >  2:905659/1‑100 (MQ=255)
          gTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtt                                                                                           <  1:1161429/100‑1 (MQ=255)
          gTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACCGCCTGCGGGCCgtt                                                                                           >  2:2859751/1‑100 (MQ=255)
            ggTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttt                                                                                         <  2:2740656/99‑1 (MQ=255)
             gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgttttt                                                                                        <  1:396343/100‑2 (MQ=255)
             gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgttttt                                                                                        <  1:468384/100‑2 (MQ=255)
             gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGAGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttt                                                                                         <  2:3743673/99‑1 (MQ=255)
              tGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgt                                                                                       >  2:2966166/1‑100 (MQ=255)
              tGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttc                                                                                       <  1:1801104/100‑3 (MQ=255)
              tGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttc                                                                                       <  1:394763/100‑3 (MQ=255)
              tGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgttttcc                                                                                       <  2:821048/100‑3 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                      >  1:3802263/1‑100 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                      <  2:1340426/100‑1 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                      <  2:244325/100‑1 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                      <  2:244339/100‑1 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgt                                                                                       >  1:2080906/1‑99 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                      <  1:302235/100‑4 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                      <  2:358788/100‑4 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGAAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                      <  2:881400/100‑4 (MQ=255)
                gAGCAGGAAGCATCGCGTAAAAACGTATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTACAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                     <  2:2229271/100‑5 (MQ=38)
                gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgttt                                                                                     <  2:549221/100‑1 (MQ=255)
                 aGCTGGAAGTATCGCGTAAATACGCATTTGTAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                      <  1:2523480/98‑4 (MQ=37)
                  gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                   <  1:4191743/100‑1 (MQ=255)
                    aGGGAGCATCGCGTAAAAACGCATTTGCAACTGTCTGCGCTTTTCCAGTATGTGGCTAAAGATTTTATGAAGAACGGCCCGCAGGCCgtttttcat                                                                                     <  2:3404786/96‑5 (MQ=255)
                    aGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTc                                                                                 <  1:1589447/100‑1 (MQ=255)
                    aGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTc                                                                                 <  1:3124081/100‑1 (MQ=255)
                    aGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcataaaa                                                                                 <  1:2135894/100‑9 (MQ=38)
                      gAAGCATCGCGTAAAAACGCAGTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTg                                                                               <  1:249719/100‑1 (MQ=255)
                       aaGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTgg                                                                              <  1:2723500/100‑1 (MQ=255)
                         gCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTgg                                                                              >  1:433419/1‑98 (MQ=255)
                           aTCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAttt                                                                          >  1:3357114/1‑100 (MQ=255)
                              gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGACAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                       >  2:851095/1‑100 (MQ=255)
                              gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                       <  2:3645355/100‑1 (MQ=255)
                              gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                       >  1:3493716/1‑100 (MQ=255)
                               cgTACAAACGCATTTGCAACTGTCGGCGCATTTCCAGTATGTTGCTAAAGATTTTATGAAAAAAGGCATGCGGGCGGTTTTGTTTTGTATGGCTTTTg                                                                        >  2:2477731/1‑98 (MQ=255)
                               cgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGTCGTTTTGTTTTGTCTGGATTTTgc                                                                       <  2:433418/99‑1 (MQ=255)
                                 tAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCt                                                                    >  1:3324644/1‑100 (MQ=255)
                                 tAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCt                                                                    >  1:3639712/1‑100 (MQ=255)
                                 tAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCt                                                                    >  2:3068600/1‑100 (MQ=255)
                                   aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttt                                                                  >  2:1916599/1‑100 (MQ=255)
                                   aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttt                                                                  <  2:96940/100‑1 (MQ=255)
                                    aaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAATGGCCTGAGGGCCGTTGTGTTTTGTCTGGATTTTGCGCtt                                                                   >  2:2817698/1‑98 (MQ=255)
                                    aaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCtttt                                                                 >  1:399872/1‑100 (MQ=255)
                                       cGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGc                                                              >  1:3548219/1‑100 (MQ=255)
                                       cGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGc                                                              >  1:1679699/1‑100 (MQ=255)
                                       cGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGc                                                              <  1:22186/100‑1 (MQ=255)
                                       cGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGc                                                              >  2:4069100/1‑100 (MQ=255)
                                       cGCATGTACAACCGTCGGCGACTCTCCAGTATGCTGATAAAGATTGTATAAAAATCGGACTACAGCCAGTTTTATTTTGTCTGGATTTTGCGCTTTTAGc                                                              <  2:3477006/100‑1 (MQ=255)
                                         cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                            <  2:1718365/100‑1 (MQ=255)
                                           tttGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAg                                                          <  1:1466456/100‑1 (MQ=255)
                                              gCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAt                                                       <  2:643328/100‑1 (MQ=255)
                                               cAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGACAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAtt                                                      >  2:513086/1‑100 (MQ=255)
                                                 aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                    >  1:338485/1‑100 (MQ=255)
                                                  cTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAg                                                   >  2:1589875/1‑100 (MQ=255)
                                                     tCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACg                                                <  1:3750604/100‑1 (MQ=255)
                                                       ggCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTTTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaa                                              >  2:1661537/1‑100 (MQ=255)
                                                       ggCGCTTTTCCAGAATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaa                                              >  2:2596816/1‑100 (MQ=255)
                                                        gcgcTTTTCCAGTATGTTGTTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTTTTTTGTTTGTATATTGCGCTTTTTGCCCAGCATTCAGACg                                                >  2:2500580/1‑97 (MQ=255)
                                                          gcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGACTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAAt                                           <  2:1410098/100‑1 (MQ=255)
                                                            ttttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTg                                         <  2:543557/100‑1 (MQ=255)
                                                            ttttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTg                                         >  1:2899366/1‑100 (MQ=255)
                                                             tttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGc                                        >  2:4178566/1‑100 (MQ=255)
                                                              ttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGcc                                       >  1:1835941/1‑100 (MQ=255)
                                                              ttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGcc                                       >  1:2586080/1‑100 (MQ=255)
                                                               tCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGccc                                      <  1:201540/100‑1 (MQ=255)
                                                                caaGTATGTTGCTAAAGATTTTATGAAAAACCGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCg                                     >  1:2826666/3‑100 (MQ=255)
                                                                  aGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCggg                                   <  1:2178847/100‑1 (MQ=255)
                                                                      tGTTGCTAACGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAtt                               >  2:1514787/1‑100 (MQ=255)
                                                                      tGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAtt                               >  2:1036625/1‑100 (MQ=255)
                                                                      tGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAtt                               >  2:1139120/1‑100 (MQ=255)
                                                                      tGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAtt                               >  1:1963391/1‑100 (MQ=255)
                                                                       gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGTCCTGCTTTCAGACGAAAATTGCCCGGGAAtt                               >  2:2992998/1‑99 (MQ=255)
                                                                       gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                              >  1:2657036/1‑100 (MQ=255)
                                                                       gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                              <  2:2113160/100‑1 (MQ=255)
                                                                        ttGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgt                             >  1:3759654/1‑100 (MQ=255)
                                                                         tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgtg                            <  1:1139120/100‑1 (MQ=255)
                                                                         tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgtg                            <  1:2913081/100‑1 (MQ=255)
                                                                          gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                           <  1:696950/100‑1 (MQ=255)
                                                                           cTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaa                          >  1:1471/1‑100 (MQ=255)
                                                                              aaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaaa                       <  1:1264555/100‑1 (MQ=255)
                                                                                gATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAAt                     >  1:2105822/1‑100 (MQ=255)
                                                                                 aTTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATa                    <  1:1916843/100‑1 (MQ=255)
                                                                                     tATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgcg                >  2:3105592/1‑100 (MQ=255)
                                                                                     tATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgcg                <  1:3737402/100‑1 (MQ=255)
                                                                                     tATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgcg                <  2:3847440/100‑1 (MQ=255)
                                                                                      aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa               >  2:3970776/5‑100 (MQ=39)
                                                                                      aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa               >  1:968721/5‑100 (MQ=39)
                                                                                       tGAAAGACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc              >  1:532687/1‑100 (MQ=255)
                                                                                       tGAAAAACGGCCTGCGGGCCGTCTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCATCATTCATACGAAAATTGGCCGGGAATTGTGAAAAAATTCGCGAc              >  2:1733079/1‑100 (MQ=255)
                                                                                        gAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa             >  1:1859811/1‑100 (MQ=255)
                                                                                          aaatCGGCCTGCGGGCCGTTTTTTTTTGTCTGGACTTTGCGCTTTTTGACCAGCATTCAGACGAAAATTCCCCGGGAATTCGGAAAAAATACGCGACAg            <  2:3322339/99‑1 (MQ=255)
                                                                                          aaaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa             >  1:284425/1‑98 (MQ=255)
                                                                                          aaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCCGCATTCAGACGAAAATTGCACGGGGATTGTGAAAAAATACGCGACAgc           >  2:795601/1‑100 (MQ=255)
                                                                                            aaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgc         >  1:2043040/1‑100 (MQ=255)
                                                                                              cGGCCTGCGGGCCATTTTGTTTTGTCTAGATATTGCGCAAATTAAACAGCATTCAGACGAAAATTGCCAAAAAAAAAAAAAAAAATACGCGACAgcgcgc       <  2:1171939/100‑1 (MQ=255)
                                                                                               ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGGATTGTGAAAAAATACGCGACAGCGCGCa      >  2:1034788/1‑100 (MQ=255)
                                                                                               ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCa      <  1:7562/100‑1 (MQ=255)
                                                                                                gCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGGCCGGGAATTGTGAAACAATACGCGACAgcg          >  2:3467181/1‑95 (MQ=255)
                                                                                                   tGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGATAGCGCGCAATaa  <  1:2750736/100‑1 (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 < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: