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% 401.0 / 67.5 178intergenic (+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 (82/67);  new base C (17/12);  total (99/79)
Fisher's exact test for biased strand distribution p-value = 8.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-02

Only 100 of 185 total aligned reads displayed.

TCTCCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATAA  >  NC_000913/454754‑454952
                                                                                                   |                                                                                                   
tctcCTCCCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCt                                                                                                     >  2:3714883/1‑100 (MQ=255)
   cctccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCg                                                                                                  >  2:125027/1‑100 (MQ=255)
     tccCTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCggg                                                                                                >  1:2903126/1‑100 (MQ=255)
        cTGTCGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCg                                                                                             <  1:2480198/100‑1 (MQ=255)
            cGTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGTTAAAGATTTTATGAAAAACGGCCTGCGGGCCgtttt                                                                                         <  2:3093445/100‑1 (MQ=255)
            cGTGGAGCAGGAAGCATAGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGAATGTGGCGAAAGACTTCATGAAAAACGGCCTGTGGGACgtgt                                                                                          >  1:3911717/1‑97 (MQ=255)
             ggggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgttttt                                                                                        <  1:2162998/98‑2 (MQ=255)
             gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttg                                                                                        >  1:3023766/1‑100 (MQ=255)
             gTGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTAAGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgtt                                                                                           >  1:848545/1‑97 (MQ=255)
              tGGAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttc                                                                                       <  1:1309751/100‑3 (MQ=255)
               ggAGCAGGAAGCATCGCGTAGAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                      >  1:2693996/1‑100 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtt                                                                                      >  2:3452711/1‑100 (MQ=255)
               ggAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttca                                                                                      <  2:2528284/100‑4 (MQ=255)
                gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgttt                                                                                     <  2:1862059/100‑1 (MQ=255)
                gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                     <  2:3016309/100‑5 (MQ=39)
                gAGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcat                                                                                     <  2:1007098/100‑5 (MQ=39)
                 ngCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                    >  2:2808777/2‑100 (MQ=255)
                 aGCAGGAAGCATCGGGTAAAAACGCAGTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                    <  2:2912167/100‑1 (MQ=255)
                 aGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCgttttgtttt                                                                                    >  2:3488687/1‑100 (MQ=255)
                 aGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcata                                                                                    <  1:688228/100‑6 (MQ=39)
                 aGCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcata                                                                                    <  2:1346438/100‑6 (MQ=39)
                  gCAGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                   >  2:3590522/1‑100 (MQ=255)
                    agcaaGCATCGCGTAAAAACGCATTTGCAACTGTCGGCTCTTTTCCATTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGttttgttttg                                                                                   <  1:1509292/95‑1 (MQ=255)
                    aGGAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCCGCAGGCCgtttttcataaaa                                                                                 <  1:3295560/100‑9 (MQ=38)
                    aGGAAGCATCACGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTCGCTACAGATTTTATGAAAAACCGCCTGCGGGCCGTTTTGTTTTGTc                                                                                 >  1:92933/1‑100 (MQ=255)
                     ggAAGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCt                                                                                >  2:3514857/1‑100 (MQ=255)
                      gAAGCATCGCGTAAAAACGCATTTGCACCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTg                                                                               >  2:115088/1‑100 (MQ=255)
                       aaGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTgg                                                                              <  2:3288839/100‑1 (MQ=255)
                       aaGCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTgg                                                                              <  2:747471/100‑1 (MQ=255)
                         gCATCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGa                                                                             <  2:3178800/99‑1 (MQ=255)
                           aTCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTTGa                                                                             >  2:2395520/1‑97 (MQ=255)
                            tCGCGTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGAtttt                                                                         >  2:1994515/1‑100 (MQ=255)
                             cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                        >  2:538938/1‑100 (MQ=255)
                             cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                        <  2:64786/100‑1 (MQ=255)
                             cgcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTg                                                                        >  2:3843671/1‑100 (MQ=255)
                              gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                       <  2:1014273/100‑1 (MQ=255)
                              gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                       <  2:1624583/100‑1 (MQ=255)
                              gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                       >  1:2333298/1‑100 (MQ=255)
                              gcgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgc                                                                       >  1:1850360/1‑100 (MQ=255)
                               cgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgcg                                                                      >  1:411778/1‑100 (MQ=255)
                               cgTAAAAACGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTgcg                                                                      >  2:4112518/1‑100 (MQ=255)
                                   aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGTGCCGTTTTGTTTTGTCTTGATTTTGCGCttt                                                                  <  2:2965426/100‑1 (MQ=255)
                                   aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttt                                                                  >  2:21215/1‑100 (MQ=255)
                                   aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCttt                                                                  >  2:21576/1‑100 (MQ=255)
                                   aaaaCGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCtt                                                                   >  2:2755112/1‑99 (MQ=255)
                                       cGCATTTGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGc                                                              >  2:1305087/1‑100 (MQ=255)
                                           tttGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAg                                                          >  2:3937261/1‑100 (MQ=255)
                                            ttGCAACTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGc                                                         <  2:398802/100‑1 (MQ=255)
                                              gCAACTGTCGGCGCTTTTCCAGTATTTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCAt                                                       <  1:868728/100‑1 (MQ=255)
                                                 aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGTTTTTTGCCCAGCATTCa                                                    <  1:1635715/100‑1 (MQ=255)
                                                 aCTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCa                                                    <  2:1934190/100‑1 (MQ=255)
                                                  cTGTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGACCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAg                                                   <  2:412145/100‑1 (MQ=255)
                                                    gTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGAc                                                 >  2:896037/1‑100 (MQ=255)
                                                    gTCGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGAc                                                 <  1:1186553/100‑1 (MQ=255)
                                                      cGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGa                                               >  1:3681038/1‑100 (MQ=255)
                                                      cGGCGCTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGACTGCGGGCCGTTTTGTTTTGTCTGGATTTTGAGCTCTTTGACCAGCGTTCAGACGa                                               >  2:3672372/1‑100 (MQ=255)
                                                        gcgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaaa                                             <  2:3387267/100‑1 (MQ=255)
                                                         cgcTTTTCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGaaaa                                            >  1:2466765/1‑100 (MQ=255)
                                                             tttCCAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTg                                         <  1:675449/99‑1 (MQ=255)
                                                             tttCCAGTATGTTACTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGc                                        >  2:2363771/1‑100 (MQ=255)
                                                                ccAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCg                                     >  2:1101142/1‑100 (MQ=255)
                                                                ccAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCg                                     <  1:1791368/100‑1 (MQ=255)
                                                                 cAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTTTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGccc                                      >  1:531587/1‑98 (MQ=255)
                                                                 cAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCgg                                    <  1:3156607/100‑1 (MQ=255)
                                                                 cAGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCgg                                    <  2:563015/100‑1 (MQ=255)
                                                                  aGTATGTTGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCggg                                   >  1:4240373/1‑100 (MQ=255)
                                                                         tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTTTGAATTTTGCACTTTTTGCCCAGCATTCAGCCGAAAATTGACCGGGAATtgt                             >  1:1077694/1‑99 (MQ=255)
                                                                         tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgtg                            <  2:383636/100‑1 (MQ=255)
                                                                         tGCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATtgtg                            >  2:3764962/1‑100 (MQ=255)
                                                                          gCTAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGa                           >  1:3861088/1‑100 (MQ=255)
                                                                            tAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaa                         >  2:2068351/1‑100 (MQ=255)
                                                                            tAAAGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaa                         <  2:871246/100‑1 (MQ=255)
                                                                               aGATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaaaa                      <  2:1390456/100‑1 (MQ=255)
                                                                                gATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGaaaaaa                      <  2:2938927/99‑1 (MQ=255)
                                                                                gATTTTATGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAAt                     >  1:2097384/1‑100 (MQ=255)
                                                                                cagacaaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCTGCATTTAGACGAAAATTGACCGGGTATTGTGAAAAAAt                     >  1:778374/11‑100 (MQ=25)
                                                                                 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)
                                                                                    caaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATAc                   >  1:3261061/7‑98 (MQ=38)
                                                                                    caaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCACGGGAATTGTGAAAAAATAcgc                 >  2:1571611/7‑100 (MQ=38)
                                                                                    caaaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCGTATTGCCCCGGATTCAGACGATAATTTCCCGGGAATTGTGAAAAAATAcgc                 >  1:2111653/7‑100 (MQ=25)
                                                                                      aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGa               >  2:1949469/5‑100 (MQ=39)
                                                                                      aaacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAACTGTGAAAAAATACGCGa               >  1:3999406/5‑100 (MQ=39)
                                                                                       tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTTTTCTTCTTTCAGACGAAAATTGCCCGGGAAGTGTGAGAAAATAcgcg                >  1:1917833/1‑98 (MQ=255)
                                                                                       tGAAAAACGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc              >  1:1717971/1‑100 (MQ=255)
                                                                                       aacaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc              >  2:4203271/4‑100 (MQ=255)
                                                                                        acaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc              >  2:1653635/3‑99 (MQ=255)
                                                                                        acaaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACa             >  2:3742200/3‑100 (MQ=255)
                                                                                         caaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATGGTGAAAAAATACGCAACa             >  1:3896314/2‑99 (MQ=255)
                                                                                         caaaaCGGCCCGCAGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTACCCGGGAATTGTGAAAAAATACGCGACAg            >  1:1880860/2‑100 (MQ=255)
                                                                                          aaaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGAc              >  1:3475458/1‑97 (MQ=255)
                                                                                           aaaCGGCCTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAgcg          >  1:1731280/1‑100 (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)
                                                                                                  cTGCGGGCCGTTTTGTTTTGTCTGGATTTTGCGCTTTTTGCCCAGCATTCAGACGAAAATTGCCCGGGAATTGTGAAAAAATACGCGACAGCGCGCAATa   <  2:1593359/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: