Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 454,856 G→A 16.1% 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,8560GA16.1% 427.5 / 63.1 179intergenic (+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 (81/68);  new base A (17/12);  total (99/80)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

Only 100 of 185 total aligned reads displayed.

CCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCG  >  NC_000913/454757‑454955
                                                                                                   |                                                                                                   
cctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCg                                                                                                     >  1:4332545/1‑100 (MQ=255)
  tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCggg                                                                                                   >  1:2903126/1‑100 (MQ=255)
    ccTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGcc                                                                                                 <  2:1034213/100‑1 (MQ=255)
    ccTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGcc                                                                                                 >  2:3781200/1‑100 (MQ=255)
        tCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttt                                                                                             >  1:1637412/1‑100 (MQ=255)
         cGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGTTAAAGATTTTATGAAAAACGGCCTGCGGGCCgtttt                                                                                            <  2:3093445/100‑1 (MQ=255)
         cGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgtttt                                                                                            >  2:1600683/1‑100 (MQ=255)
          ggggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttt                                                                                            <  1:1917648/97‑1 (MQ=255)
          gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttg                                                                                           >  1:3023766/1‑100 (MQ=255)
            ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                         >  2:3452711/1‑100 (MQ=255)
            ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                         <  2:2528284/100‑4 (MQ=255)
             gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                        <  2:1007098/100‑5 (MQ=39)
             gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                        <  2:3016309/100‑5 (MQ=39)
              ngCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                       >  2:2808777/2‑100 (MQ=255)
              aGCAGGAAGCATCGGGTAAAAACGCAGTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                       <  2:2912167/100‑1 (MQ=255)
              aGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                       >  1:2247038/1‑100 (MQ=255)
              aGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                       >  2:3488687/1‑100 (MQ=255)
              aGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcata                                                                                       <  2:1346438/100‑6 (MQ=39)
               gCAGGAAGCATCGCGTACAAACGCAGTTGCAACCAACGGCGCGTTGCCCTTATGTGGTGAAAGATTTTATCAAAAACGGCTTGAGGGCCGttatggtttg                                                                                      >  1:172240/1‑100 (MQ=255)
               gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                      >  2:3590522/1‑100 (MQ=255)
                 aGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTGCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTc                                                                                    <  1:1849841/100‑1 (MQ=255)
                  ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                                   >  2:3514857/1‑100 (MQ=255)
                  ggAAGCATCGCGTAAAAACGCATTTGCAACTGGCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                                   <  1:3621130/100‑1 (MQ=255)
                   gAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTg                                                                                  <  1:588720/100‑1 (MQ=255)
                    aaGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTgg                                                                                 <  2:747471/100‑1 (MQ=255)
                    aaGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTgg                                                                                 <  2:3288839/100‑1 (MQ=255)
                      gCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGa                                                                                <  2:3178800/99‑1 (MQ=255)
                        aTCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTTGa                                                                                >  2:2395520/1‑97 (MQ=255)
                        aTCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAttt                                                                             >  1:211995/1‑100 (MQ=255)
                        aTCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAttt                                                                             >  1:3961711/1‑100 (MQ=255)
                         tCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAtttt                                                                            >  2:1994412/1‑100 (MQ=255)
                         tCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAtttt                                                                            <  2:2102592/100‑1 (MQ=255)
                          ctcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                           <  1:2676598/98‑1 (MQ=255)
                          cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                           <  2:1584807/100‑1 (MQ=255)
                          cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                           <  2:64786/100‑1 (MQ=255)
                          cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                           <  2:3731118/100‑1 (MQ=255)
                          cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                           >  2:538938/1‑100 (MQ=255)
                           gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                          <  2:1013931/100‑1 (MQ=255)
                            cgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgcg                                                                         >  2:4112518/1‑100 (MQ=255)
                                aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGTGCCGTTTTGTTTTGTCTTGATTTTGCGCttt                                                                     <  2:2965426/100‑1 (MQ=255)
                                aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttt                                                                     >  2:21576/1‑100 (MQ=255)
                                aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCtt                                                                      >  2:2755112/1‑99 (MQ=255)
                                      cATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                               >  1:2383757/1‑100 (MQ=255)
                                       aTTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGccc                                                               >  1:2472476/1‑99 (MQ=255)
                                       aTTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGc                                                                 <  1:3458933/97‑1 (MQ=255)
                                        tttGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAg                                                             <  1:1548861/100‑1 (MQ=255)
                                         ttGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGc                                                            <  2:398802/100‑1 (MQ=255)
                                          tGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCa                                                           >  2:3755141/1‑100 (MQ=255)
                                              aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                       <  2:2205205/100‑1 (MQ=255)
                                               cTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGACCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAg                                                      <  2:412145/100‑1 (MQ=255)
                                                 gTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGAc                                                    >  2:896037/1‑100 (MQ=255)
                                                   cGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGa                                                  >  1:144067/1‑100 (MQ=255)
                                                   cGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGACTGCGGGCCGTTTTGTTTTGTCTGGATTTTGAGCTCTTTGACCAGCGTTCAGACGa                                                  >  2:3672372/1‑100 (MQ=255)
                                                    ggCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaa                                                 <  1:3909608/100‑1 (MQ=255)
                                                    ggCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTGTGTTTTGTTTGGATTTTGCGCTTTTTGCCCAGCATTCAGAc                                                    >  1:695694/1‑97 (MQ=255)
                                                     gcgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaaa                                                <  2:3387267/100‑1 (MQ=255)
                                                      cgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGCCCGTTTTGATTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaaaa                                               >  1:915803/1‑100 (MQ=255)
                                                       gcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAAt                                              >  2:126456/1‑100 (MQ=255)
                                                         ttttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTg                                            <  1:675812/100‑1 (MQ=255)
                                                             ccAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCg                                        >  2:1924807/1‑100 (MQ=255)
                                                             ccAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCg                                        >  2:1101142/1‑100 (MQ=255)
                                                             ccAGTATGATGCTAAAGATTTTATAAAAAGAGGCCTGCTGGACGTTTTGTTTTGTGTGGATTTTCCGCTTTTTGACCAGCATTCAGACGAAAAt                                              >  1:3282544/1‑94 (MQ=255)
                                                              cAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCgg                                       <  2:563015/100‑1 (MQ=255)
                                                              cAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCgg                                       <  1:3156607/100‑1 (MQ=255)
                                                                   tGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAAtt                                  >  2:1182030/1‑100 (MQ=255)
                                                                    gTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtg                                 <  1:84880/100‑1 (MQ=255)
                                                                      tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTTTGAATTTTGCACTTTTTGCCCAGCATTCAGCCGAAAATTGACCGGGAATtgt                                >  1:1077694/1‑99 (MQ=255)
                                                                       gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                              <  1:412345/100‑1 (MQ=255)
                                                                       gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                              >  1:3861088/1‑100 (MQ=255)
                                                                        cTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaa                             <  2:1947058/100‑1 (MQ=255)
                                                                        cTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaa                             <  1:816692/100‑1 (MQ=255)
                                                                         tAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaa                            <  2:871246/100‑1 (MQ=255)
                                                                          aaaGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGACAATTGCCCGGGAATTGTGaaaa                           >  1:119011/1‑100 (MQ=255)
                                                                          aaaGATTTTATGAAAAACGGCCTGCGGGCCGTATTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaa                           >  1:1986836/1‑100 (MQ=255)
                                                                            aGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaaaa                         <  2:1390456/100‑1 (MQ=255)
                                                                             gATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaaaa                         <  2:2938927/99‑1 (MQ=255)
                                                                              aTTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATa                       <  2:4292824/100‑1 (MQ=255)
                                                                               ttttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAc                      <  2:1447882/100‑1 (MQ=255)
                                                                               ttttATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAc                      >  1:567843/1‑100 (MQ=255)
                                                                                 caaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgc                    >  2:442887/7‑100 (MQ=38)
                                                                                 caaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCGTATTGCCCCGGATTCAGACGATAATTTCCCGGGAATTGTGAAAAAATAcgc                    >  1:2111653/7‑100 (MQ=25)
                                                                                  tATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAcgcg                   <  2:1815734/100‑1 (MQ=255)
                                                                                    aacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc                 >  2:4203271/4‑100 (MQ=255)
                                                                                     gAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa                >  2:1997082/1‑100 (MQ=255)
                                                                                     acaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc                 >  2:1653635/3‑99 (MQ=255)
                                                                                      caaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAg               >  2:3854735/2‑100 (MQ=255)
                                                                                        aaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcg             >  1:1731924/1‑100 (MQ=255)
                                                                                        aaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcg             >  1:1859287/1‑100 (MQ=255)
                                                                                        aaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCACGGGAATTGTGAAAAAATACGCGACAg               >  1:387072/1‑98 (MQ=255)
                                                                                          aCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcgcg           <  2:1486601/100‑1 (MQ=255)
                                                                                            ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCa         <  2:2694011/100‑1 (MQ=255)
                                                                                            ggCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCa         <  1:1069217/100‑1 (MQ=255)
                                                                                              ccTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAAt       >  2:405453/1‑100 (MQ=255)
                                                                                              ccTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAAt       <  1:1720496/100‑1 (MQ=255)
                                                                                                tGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATaa     <  2:3837190/100‑1 (MQ=255)
                                                                                                tGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATaa     >  1:4031438/1‑100 (MQ=255)
                                                                                                 gCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAc    >  2:4034318/1‑100 (MQ=255)
                                                                                                 gCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAc    >  2:238245/1‑100 (MQ=255)
                                                                                                  cGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAAcc   >  2:151277/1‑100 (MQ=255)
                                                                                                   gggCCGTTTTGTTGTGTCTGGATTTTGCGCTTTTTGCCAAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAACCg  >  1:2717000/1‑100 (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: