Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 89,815 C→T S61L (TCG→TTG)  yabB → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091389,8150CT97.6% 579.3 / ‑0.4 167S61L (TCG→TTG) yabBhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (3/1);  new base T (87/76);  total (90/77)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 177 total aligned reads displayed.

CCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGAT  >  NC_000913/89725‑89907
                                                                                          |                                                                                            
ccGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATtgc                                                                                            >  1:4182474/1‑93 (MQ=255)
ccGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATtgc                                                                                            >  1:3873576/1‑93 (MQ=255)
ccGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATtgc                                                                                            >  2:2465026/1‑93 (MQ=255)
 cGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATtgcg                                                                                           <  1:972694/93‑1 (MQ=255)
  ggTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGt                                                                                          <  1:1747673/93‑1 (MQ=255)
  ggTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGt                                                                                          >  2:361313/1‑93 (MQ=255)
   gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                                         <  2:4361105/93‑1 (MQ=255)
   gTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTc                                                                                         <  2:3287848/93‑1 (MQ=255)
    tCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCt                                                                                        >  1:1488889/1‑93 (MQ=255)
    tCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCt                                                                                        <  1:2108564/93‑1 (MQ=255)
      aaaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGt                                                                                      >  1:1077621/1‑93 (MQ=255)
      aaaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGt                                                                                      >  2:2287434/1‑93 (MQ=255)
       aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTTCCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTc                                                                                     >  1:1066940/1‑93 (MQ=255)
       aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGGAATTATCAAGCAAAAATTATTGCGTCTGTc                                                                                     >  2:3851555/1‑93 (MQ=255)
       aaTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTc                                                                                     <  2:2359463/93‑1 (MQ=255)
        aTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCg                                                                                    >  2:4039804/1‑93 (MQ=255)
        aTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCg                                                                                    >  1:1636904/1‑93 (MQ=255)
        aTGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCg                                                                                    <  1:2901897/93‑1 (MQ=255)
         tGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGa                                                                                   <  2:1304398/93‑1 (MQ=255)
         tGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGa                                                                                   >  1:3017054/1‑93 (MQ=255)
         tGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGa                                                                                   <  1:3043971/93‑1 (MQ=255)
            tttGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCa                                                                                >  1:4374384/1‑93 (MQ=255)
             ttGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCAt                                                                               <  2:4194475/93‑1 (MQ=255)
                 aCCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAc                                                                           >  2:1716304/1‑93 (MQ=255)
                  ccATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAcc                                                                          >  1:2575963/1‑93 (MQ=255)
                   cATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAccc                                                                         <  2:3659986/93‑1 (MQ=255)
                   cATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAAccc                                                                         <  2:605971/93‑1 (MQ=255)
                    aTTGACATTTATCACCCGTGCCTGCTGATTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCg                                                                        >  1:89107/1‑93 (MQ=255)
                     ttGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCTTCTGTCGAGCATGAACCCgg                                                                       <  1:849319/93‑1 (MQ=255)
                     ttGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCgg                                                                       >  2:1892883/1‑93 (MQ=255)
                     ttGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCgg                                                                       >  1:4693420/1‑93 (MQ=255)
                      tGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGt                                                                      <  2:2315741/93‑1 (MQ=255)
                         cATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGa                                                                   <  1:241444/93‑1 (MQ=255)
                         cATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGa                                                                   <  2:1568496/93‑1 (MQ=255)
                          aTTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAg                                                                  >  2:3789435/1‑93 (MQ=255)
                          aTTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAg                                                                  >  2:3789434/1‑93 (MQ=255)
                            ttATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAgcg                                                                <  2:2202859/93‑1 (MQ=255)
                              aTCACCCGTGCCTGCTGCTTTACACCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATCAACCCGGTTGAGCGcc                                                              >  1:2839724/1‑93 (MQ=255)
                               tCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCg                                                             <  2:4243120/93‑1 (MQ=255)
                                 ccccGTTCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgtg                                                           <  2:537214/92‑1 (MQ=255)
                                  cccGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCgtgt                                                          >  2:3627776/1‑93 (MQ=255)
                                    cGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGc                                                        <  2:648734/93‑1 (MQ=255)
                                       gCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgc                                                     >  1:4012872/1‑93 (MQ=255)
                                       gCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgc                                                     >  2:383035/1‑93 (MQ=255)
                                        ccTGCTGCTTTACCCCCTGCCTGCATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGATGAGCGCCGTGCGCAgcg                                                    >  2:146712/1‑93 (MQ=255)
                                         ctgctgCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgcgc                                                   >  2:4313132/1‑93 (MQ=255)
                                         ctgctgCTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAgcgc                                                   >  2:377414/1‑93 (MQ=255)
                                          tgctgcTTTCCCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGcc                                                  <  1:1857497/93‑1 (MQ=255)
                                          tgctgcTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGcc                                                  >  2:3941383/1‑93 (MQ=255)
                                          tgctgcTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGcc                                                  >  2:1130540/1‑93 (MQ=255)
                                          tgctgcTTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGcc                                                  >  2:1037939/1‑93 (MQ=255)
                                           gctgctTTACCCCCTGCCTGACTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCt                                                 >  2:4658827/1‑93 (MQ=255)
                                             tgctTTACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTAc                                               <  1:936198/93‑1 (MQ=255)
                                                tttACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGt                                            <  1:1217508/93‑1 (MQ=255)
                                                tttACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTAGTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGt                                            <  2:2107311/93‑1 (MQ=255)
                                                 ttACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGtt                                           <  2:703575/93‑1 (MQ=255)
                                                  tACCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTa                                          <  2:2833837/93‑1 (MQ=255)
                                                   tcccccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAg                                         <  2:2213346/92‑1 (MQ=255)
                                                   aCCCCCTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAg                                         <  2:1927550/93‑1 (MQ=255)
                                                    cccccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAgg                                        <  1:961799/93‑1 (MQ=255)
                                                    cccccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAgg                                        <  1:3521907/93‑1 (MQ=255)
                                                     ccccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGt                                       >  2:2663572/1‑93 (MQ=255)
                                                     ccccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGt                                       >  1:3048501/1‑93 (MQ=255)
                                                     ccccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGATGAACGCCGTGTGCAGCGCCTACTGTTAggc                                       >  2:54973/1‑92 (MQ=255)
                                                      cccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTc                                      >  2:2634639/1‑93 (MQ=255)
                                                      cccTGCCTGAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTc                                      <  1:3670234/93‑1 (MQ=255)
                                                       cctgcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCa                                     <  2:824874/93‑1 (MQ=255)
                                                       cctgcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCa                                     >  1:3866399/1‑93 (MQ=255)
                                                          gcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGc                                  >  1:785489/1‑93 (MQ=255)
                                                          gcctgAATGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGc                                  <  1:4110206/93‑1 (MQ=255)
                                                              gAAGGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGc                              <  2:2774810/93‑1 (MQ=255)
                                                               aaTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCg                             >  1:2473028/1‑93 (MQ=255)
                                                                aTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGa                            >  1:140709/1‑93 (MQ=255)
                                                                aTGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGCCATGCCAGCGa                            >  1:745860/1‑93 (MQ=255)
                                                                 tGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGaa                           >  1:1714613/1‑93 (MQ=255)
                                                                 tGGGAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGaa                           <  1:3284671/93‑1 (MQ=255)
                                                                  gggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAAt                          >  1:1410982/1‑93 (MQ=255)
                                                                  gggAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAAt                          <  1:609220/93‑1 (MQ=255)
                                                                    gAAATTATCAAGCAAAAATTATTGCGTCTGTCGAGAATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGt                        >  2:2146492/1‑93 (MQ=255)
                                                                     aaaTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTc                       <  2:4413404/93‑1 (MQ=255)
                                                                      aaTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCa                      >  1:2657091/1‑93 (MQ=255)
                                                                      aaTTATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCa                      >  2:1332835/1‑93 (MQ=255)
                                                                       aTTATCAAGCAAAAATTATTGCGTCTGTCGAGCCTGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAg                     >  1:1632296/1‑93 (MQ=255)
                                                                        ttATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAga                    <  2:2807673/93‑1 (MQ=255)
                                                                        ttATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAga                    <  2:2358013/93‑1 (MQ=255)
                                                                        atATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAga                    <  2:2807674/92‑1 (MQ=255)
                                                                         tATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgat                   <  1:3374001/93‑1 (MQ=255)
                                                                         tATCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgat                   <  1:926071/93‑1 (MQ=255)
                                                                          aTCAAGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatg                  >  2:2756045/1‑93 (MQ=255)
                                                                           tcaaGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatgg                 >  2:3091492/4‑93 (MQ=255)
                                                                             aaGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatggat               >  2:2416631/2‑93 (MQ=255)
                                                                              aGCAAAAATTATTGCGTCTGTCGAGCATGAACGCGGTTGAGCGCCGTGTGCAGCGCCTAATGTTAGGTCATGCCAGCGACTGTCAgatggatg              >  1:378844/1‑93 (MQ=255)
                                                                              aGCAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAgatggatg              <  1:1322805/93‑1 (MQ=255)
                                                                                cAAAAATTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGgc            <  2:2013537/93‑1 (MQ=255)
                                                                                 aaaaaTTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGgcg           >  1:3450244/1‑93 (MQ=255)
                                                                                  aaaaTTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGgcgc          <  1:3154559/93‑1 (MQ=255)
                                                                                   aaaTTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCa         >  1:4483142/1‑93 (MQ=255)
                                                                                    aaTTATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCCGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAg        <  1:56027/93‑1 (MQ=255)
                                                                                       tATTGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTc     >  1:2825663/1‑93 (MQ=255)
                                                                                          cgcgTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGAt  <  1:3843698/93‑1 (MQ=255)
                                                                                          |                                                                                            
CCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATCGAGCAAAAATTATCGCGTCTGTCGAGCATGAACCCGGTTGAGCGCCGTGTGCAGCGCCTACTGTTAGGTCATGCCAGCGAATGTCAGATGGATGGCGCAGGTCGAT  >  NC_000913/89725‑89907

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: