Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 454,853 T→C 15.6% 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,8530TC15.6% 397.8 / 58.5 172intergenic (+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 (76/69);  new base C (15/12);  total (91/81)
Fisher's exact test for biased strand distribution p-value = 8.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 172 total aligned reads displayed.

CTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATA  >  NC_000913/454755‑454951
                                                                                                  |                                                                                                  
 tcctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGc                                                                                                  >  1:1073041/1‑100 (MQ=255)
    tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCggg                                                                                               <  2:586727/100‑1 (MQ=255)
     cccTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGc                                                                                              >  2:581756/1‑100 (MQ=255)
      ccTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGcc                                                                                             >  1:588584/1‑100 (MQ=255)
         gTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgtt                                                                                          >  2:3367326/1‑100 (MQ=255)
         gTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgtt                                                                                          >  1:2648611/1‑100 (MQ=255)
            gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttg                                                                                       <  2:3668774/100‑1 (MQ=255)
              ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTTCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                     >  1:1700222/1‑100 (MQ=255)
              ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                     >  2:2133908/1‑100 (MQ=255)
              ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                     <  1:2880758/100‑1 (MQ=255)
              ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                     <  2:559304/100‑4 (MQ=255)
               gAGCAGGAATGATCGCCTAAAAACGCATTTGCAACTGTCGGCGCTTTTGCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                     <  1:3328001/99‑4 (MQ=38)
               gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgttt                                                                                    <  1:1500118/100‑1 (MQ=255)
               gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                    <  2:1285281/100‑5 (MQ=39)
                 gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                  <  2:674058/100‑1 (MQ=255)
                 gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcataa                                                                                  <  2:1180905/100‑7 (MQ=38)
                  cAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcataaa                                                                                 <  2:1364417/100‑8 (MQ=38)
                   aGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTc                                                                                <  1:712139/100‑1 (MQ=255)
                    ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                               <  1:1487609/100‑1 (MQ=255)
                    ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                               >  1:283189/1‑100 (MQ=255)
                     gAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTg                                                                              <  2:764339/100‑1 (MQ=255)
                     gAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTg                                                                              <  2:3685121/100‑1 (MQ=255)
                     gAAGCATCGAGTAAAACTGCATTTGCAACTGTCGGCTCCTTTCAAGTATGGTGCTAACGATTTTATGCAAAACGGCCCGCAGGCCgtttttcata                                                                                   <  1:2307625/95‑6 (MQ=19)
                      aaGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTgg                                                                             <  2:3700111/100‑1 (MQ=255)
                          aTCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCGGGCCGTTTTGTTTTGTCTGGAttt                                                                         >  1:617383/1‑100 (MQ=255)
                            cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                       <  2:2907703/100‑1 (MQ=255)
                            cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                       <  2:1983638/100‑1 (MQ=255)
                            cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                       >  2:1286298/1‑100 (MQ=255)
                             gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                      >  1:2480487/1‑100 (MQ=255)
                             gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                      >  2:1740308/1‑100 (MQ=255)
                             gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTCGTCTgaatgctgg                                                                      >  2:1399099/1‑92 (MQ=255)
                              cgTAAAAACGCAGTTGGAACTGTCGGGGCTTTTTCAGTATGTTGCTAAAGATTGTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgcg                                                                     <  1:1357994/100‑1 (MQ=255)
                                tAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCt                                                                   >  1:180249/1‑100 (MQ=255)
                                  aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttt                                                                 >  2:578710/1‑100 (MQ=255)
                                   aaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCtttt                                                                >  2:2369173/1‑100 (MQ=255)
                                       gCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGcc                                                            <  2:3689845/100‑1 (MQ=255)
                                        cATTTTCAACTGTCGGCGCTTTTCCAGGATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                           <  1:934189/100‑1 (MQ=255)
                                        cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                           >  1:2378313/1‑100 (MQ=255)
                                        cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                           <  2:2975536/100‑1 (MQ=255)
                                        cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                           >  2:3592627/1‑100 (MQ=255)
                                         aTTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCa                                                          >  2:2516839/1‑100 (MQ=255)
                                            tGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCa                                                       >  2:354908/1‑100 (MQ=255)
                                            tGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCa                                                       >  2:441817/1‑100 (MQ=255)
                                             gCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAt                                                      <  1:341032/100‑1 (MQ=255)
                                               aaCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTc                                                    <  1:3698934/100‑1 (MQ=255)
                                                   gTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGAc                                                <  1:218870/100‑1 (MQ=255)
                                                     cGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGa                                              >  1:1265333/1‑100 (MQ=255)
                                                      ggCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaa                                             >  2:2397202/1‑100 (MQ=255)
                                                        cgcTTTTCCAGTATGTTTATAAAGATTTTATTAAAATCGGCCTGCGGGCCGCTTTGTTTTGTCTGAATTATACACTTTATGCCCAGCATTCAGAc                                                <  1:1209860/95‑1 (MQ=255)
                                                         gcTTTTCCAGTCTGTTGCTAAAGATTTTATGAAACACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAAt                                          >  1:1994878/1‑100 (MQ=255)
                                                          cTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAAtt                                         >  1:2772008/1‑100 (MQ=255)
                                                          cTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAAtt                                         >  1:3101890/1‑100 (MQ=255)
                                                           ttttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAAtt                                         >  1:351841/1‑99 (MQ=255)
                                                           ttttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTg                                        >  2:3299191/1‑100 (MQ=255)
                                                                cAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGTCCAGCATTCAGACGAAAATTGCCCgg                                   <  1:1286290/100‑1 (MQ=255)
                                                                cAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCgg                                   <  2:3525883/100‑1 (MQ=255)
                                                                   tATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGTCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCgg                                   >  1:1104107/1‑97 (MQ=255)
                                                                    aTGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAt                               <  1:1975324/100‑1 (MQ=255)
                                                                      gTTGCTAATGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                             <  1:519773/100‑1 (MQ=255)
                                                                      gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAtt                              >  1:3226548/1‑99 (MQ=255)
                                                                      gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                             <  2:336186/100‑1 (MQ=255)
                                                                       ttGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgt                            <  2:2461054/100‑1 (MQ=255)
                                                                        tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgtg                           <  2:3118540/100‑1 (MQ=255)
                                                                         gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                          >  1:1058544/1‑100 (MQ=255)
                                                                          cTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGGCGGGATGTTACGCTTTTTGAAAAGAATTCAGACGAAAATAGCAAGGAAATtg                             >  2:3428123/1‑96 (MQ=255)
                                                                           tAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaa                        <  2:3549649/100‑1 (MQ=255)
                                                                            aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaa                       <  1:941472/100‑1 (MQ=255)
                                                                              aGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaaaa                     >  2:1722560/1‑100 (MQ=255)
                                                                               gATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGCATTGTGAAAAAAt                    >  1:1157896/1‑100 (MQ=255)
                                                                               gATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAAt                    <  2:2439684/100‑1 (MQ=255)
                                                                               gATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAAt                    >  2:674380/1‑100 (MQ=255)
                                                                               cagacaaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTTTCCAGCATTCAGACGAAAATTGCCCGGGCATTGTGAAAAAAt                    >  1:3614971/11‑100 (MQ=25)
                                                                                aTTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATa                   <  2:858791/100‑1 (MQ=255)
                                                                                 ttttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAc                  >  1:1923263/1‑100 (MQ=255)
                                                                                 ttttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACAAAAATTGCCCGGGAATTGTGAAAAAATa                   >  1:2597758/1‑99 (MQ=255)
                                                                                 ttttATGAAAAACGGCCTGCGGGACGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGAATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATa                   >  1:2335383/1‑99 (MQ=255)
                                                                                  acaaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCGAGTATTCAGACGAAAATTGCCCGTGCATTGTGAAAAAATAcg                 >  1:882653/8‑100 (MQ=37)
                                                                                   caaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAagc                >  1:3757618/7‑100 (MQ=38)
                                                                                    tATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgcg               <  1:1845562/100‑1 (MQ=255)
                                                                                    aaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGAAATTGTGAAAAAATAcgcg               >  2:2220110/6‑100 (MQ=38)
                                                                                     aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa              >  2:1532063/5‑100 (MQ=39)
                                                                                     aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa              >  2:2707605/5‑100 (MQ=39)
                                                                                     aTGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTTTTTTTTTTTTAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgcg               >  1:1198990/1‑99 (MQ=255)
                                                                                     aTGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTATGAAAAAATACGCGa              >  2:3756470/1‑100 (MQ=255)
                                                                                      tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgcg               <  1:772474/98‑1 (MQ=255)
                                                                                      tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc             >  2:198669/1‑100 (MQ=255)
                                                                                       gAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa            >  2:3133537/1‑100 (MQ=255)
                                                                                        caaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAg           >  1:2592222/2‑100 (MQ=255)
                                                                                        aaaaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAg           <  1:1892504/100‑1 (MQ=255)
                                                                                         aaaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgc          >  2:1831824/1‑100 (MQ=255)
                                                                                           aaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgc        <  2:1670650/100‑1 (MQ=255)
                                                                                             cGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgcg       <  1:1329704/99‑1 (MQ=255)
                                                                                             cGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCCCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgcgc      <  2:2676759/100‑1 (MQ=255)
                                                                                              ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCa     >  1:3237230/1‑100 (MQ=255)
                                                                                              ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCa     <  1:2070283/100‑1 (MQ=255)
                                                                                              ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAACACGCGACAGCGCGCa     >  2:3766420/1‑100 (MQ=255)
                                                                                               gCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCaa    >  2:1394525/1‑100 (MQ=255)
                                                                                               gCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCaa    >  2:2538596/1‑100 (MQ=255)
                                                                                                ccTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAAt   >  2:236178/1‑100 (MQ=255)
                                                                                                 cTGCGGGCCGTTTTGGTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATa  <  1:345491/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
CTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATA  >  NC_000913/454755‑454951

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: