Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,934,276 C→T H125Y (CAC→TAC)  rbsB → D‑ribose periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,934,2760CT99.3% 1649.7 / ‑3.6 452H125Y (CAC→TAC) rbsBD‑ribose periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (234/215);  total (236/216)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-01

Only 100 of 474 total aligned reads displayed.

GCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTG  >  NC_000913/3934186‑3934365
                                                                                          |                                                                                         
gCAGTGGGTAATGCTGTGAAGATGGCTAGCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCtac                                                                                         <  1:9733288/93‑1 (MQ=255)
gCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCtac                                                                                         >  1:4988349/1‑93 (MQ=255)
 cAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCtaca                                                                                        <  2:4843933/93‑1 (MQ=255)
  aGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAt                                                                                       <  2:9140821/93‑1 (MQ=255)
    tGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTg                                                                                     >  2:58408/1‑93 (MQ=255)
     gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGc                                                                                    <  1:5399796/93‑1 (MQ=255)
       gTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCtt                                                                                  <  1:8566105/93‑1 (MQ=255)
        tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTc                                                                                 >  1:6918776/1‑93 (MQ=255)
        tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTc                                                                                 <  1:9561334/93‑1 (MQ=255)
        tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTc                                                                                 >  1:3997857/1‑93 (MQ=255)
        tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTc                                                                                 <  2:2511852/93‑1 (MQ=255)
              tgtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATaa                                                                           <  1:4775451/93‑1 (MQ=255)
               gtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGTGGAAGTGGGGAGCTACATTGCGTCTGATAAc                                                                          >  1:2396433/1‑93 (MQ=255)
               gtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAAc                                                                          <  2:7771875/93‑1 (MQ=255)
               gtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAAc                                                                          <  1:415997/93‑1 (MQ=255)
                tgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACg                                                                         >  2:2221084/1‑93 (MQ=255)
                 gAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGt                                                                        >  2:5988940/1‑93 (MQ=255)
                 gAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGt                                                                        >  2:7550243/1‑93 (MQ=255)
                  aaGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTa                                                                       <  2:8450215/93‑1 (MQ=255)
                  aaGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTa                                                                       <  2:5019952/93‑1 (MQ=255)
                   aGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTAc                                                                      >  2:4002191/1‑93 (MQ=255)
                   aGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTAc                                                                      >  1:4885623/1‑93 (MQ=255)
                    gATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACt                                                                     <  1:2246566/93‑1 (MQ=255)
                    gATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACt                                                                     <  1:7937194/93‑1 (MQ=255)
                    gATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACt                                                                     >  2:7163369/1‑93 (MQ=255)
                    gATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATAGCTTCTGATAACGTACt                                                                     >  1:8953051/1‑93 (MQ=255)
                     aTGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTg                                                                    <  2:57450/93‑1 (MQ=255)
                     aTGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTg                                                                    <  1:2123983/93‑1 (MQ=255)
                         cTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGggcg                                                                <  2:4706426/93‑1 (MQ=255)
                          tAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGggcgg                                                               <  1:1596267/93‑1 (MQ=255)
                           aaCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGggcggc                                                              >  2:10006856/1‑93 (MQ=255)
                             ccAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCaa                                                            >  2:1528742/1‑93 (MQ=255)
                              cAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGGGAGCTACATTGCTTCTGATAACGTACTGGGCGGCaaa                                                           >  1:5096496/1‑93 (MQ=255)
                               cggCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCaaaa                                                          <  1:3534661/92‑1 (MQ=255)
                               aGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCaaaa                                                          >  1:5747838/1‑93 (MQ=255)
                                ggCGAACGACCCGGTGATCGCTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAAt                                                         <  1:1979795/93‑1 (MQ=255)
                                ggCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAAt                                                         <  2:3401057/93‑1 (MQ=255)
                                ggCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAAt                                                         >  1:1018424/1‑93 (MQ=255)
                                 gCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATc                                                        >  1:5016315/1‑93 (MQ=255)
                                    aaCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCt                                                     >  2:9366592/1‑93 (MQ=255)
                                     aCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTg                                                    <  1:768251/93‑1 (MQ=255)
                                       aTCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGt                                                  >  2:159467/1‑93 (MQ=255)
                                        tCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTg                                                 <  2:7684898/93‑1 (MQ=255)
                                        tCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTg                                                 >  2:1292362/1‑93 (MQ=255)
                                        tCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGGGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTg                                                 >  1:5941493/1‑93 (MQ=255)
                                            ggTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTa                                             >  1:3856669/1‑93 (MQ=255)
                                             gTTATCACTCTTGAGCGCCAGGCAACGAAAGGGGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTAc                                            >  2:9053832/1‑93 (MQ=255)
                                             gTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTAc                                            >  1:8790069/1‑93 (MQ=255)
                                             gTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTAc                                            <  2:7273795/93‑1 (MQ=255)
                                              ttATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACa                                           <  2:4430686/93‑1 (MQ=255)
                                              ttATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGCGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACa                                           >  1:6357766/1‑93 (MQ=255)
                                               tATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACAt                                          <  2:6325240/93‑1 (MQ=255)
                                                ttCACTCTTGACCGCCAGGCAACGAAAGGTGGAGTAGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATc                                         <  1:6824480/92‑1 (MQ=255)
                                                 tCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcg                                        <  2:188781/93‑1 (MQ=255)
                                                  cACTCTTTACCGCCAGTCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcgc                                       <  1:9188035/93‑1 (MQ=255)
                                                  cACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcgc                                       <  1:2530147/93‑1 (MQ=255)
                                                  cACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTAAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcgc                                       <  1:6483267/93‑1 (MQ=255)
                                                   aCTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcgcg                                      >  2:7936613/1‑93 (MQ=255)
                                                      ctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCgaag                                   >  2:520138/1‑93 (MQ=255)
                                                       ttGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCgaaga                                  <  1:1116784/93‑1 (MQ=255)
                                                       ttGACCGCCAGGCAACGAAAGGTGAAGTGGAGAGCTACATTGCTTCTGATAACATACTGGGCGGCAAAATCGCTGGAGATTACATCGCtaaga                                  >  2:5574199/1‑93 (MQ=255)
                                                        tGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCgaagaa                                 <  1:7216610/93‑1 (MQ=255)
                                                           ccGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGc                              <  2:1672794/93‑1 (MQ=255)
                                                             gCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCgg                            >  2:6745070/1‑93 (MQ=255)
                                                             gCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCgg                            >  1:2714981/1‑93 (MQ=255)
                                                              ccAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCggg                           <  2:8239078/93‑1 (MQ=255)
                                                              ccAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCggg                           <  2:6527719/93‑1 (MQ=255)
                                                              ccAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCggg                           <  2:4131461/93‑1 (MQ=255)
                                                              ccAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCggg                           >  1:1403113/1‑93 (MQ=255)
                                                               cAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCgggg                          >  2:1965627/1‑92 (MQ=255)
                                                                aGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTg                         <  1:4403882/93‑1 (MQ=255)
                                                                aGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTg                         >  1:5308422/1‑93 (MQ=255)
                                                                 ggCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGa                        <  2:3691085/93‑1 (MQ=255)
                                                                 ggCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGa                        >  1:3218761/1‑93 (MQ=255)
                                                                  gCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGaa                       >  2:5369724/1‑93 (MQ=255)
                                                                  gCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGaa                       <  1:5522192/93‑1 (MQ=255)
                                                                    aaCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAgg                     <  2:3125332/93‑1 (MQ=255)
                                                                    aaCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAgg                     <  2:9299570/93‑1 (MQ=255)
                                                                     aCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGt                    >  2:9489546/1‑93 (MQ=255)
                                                                      cGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTg                   <  1:3670148/93‑1 (MQ=255)
                                                                      cGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTg                   >  2:8302063/1‑93 (MQ=255)
                                                                       gAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGc                  >  1:6107670/1‑93 (MQ=255)
                                                                          aGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCaa               >  1:4603540/1‑93 (MQ=255)
                                                                          aGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCaa               >  1:3032369/1‑93 (MQ=255)
                                                                            gTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAg             <  2:8027688/93‑1 (MQ=255)
                                                                            gTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGGGAAGGTGCCAAAg             >  1:7450471/1‑93 (MQ=255)
                                                                             tGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGt            <  1:8269200/93‑1 (MQ=255)
                                                                             ggAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGt            <  2:1415377/92‑1 (MQ=255)
                                                                              gAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGtt           >  2:5653738/1‑93 (MQ=255)
                                                                               taGTGGTGTGCTACATTGCTTCTGATTACGTACTGGGCGGCAAAATCGCTGGTGCTTACATCGCGAGGAAAGCGGGTGAAGGTGCCAAAGTTa          <  2:8981129/92‑1 (MQ=255)
                                                                               aaGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTa          <  1:2919713/93‑1 (MQ=255)
                                                                               aaGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTa          <  2:333046/93‑1 (MQ=255)
                                                                                 gtggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc        <  2:870886/93‑1 (MQ=255)
                                                                                 gtggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc        <  2:4796908/93‑1 (MQ=255)
                                                                                 gtggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc        >  2:3036701/1‑93 (MQ=255)
                                                                                 gtggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc        <  2:9657888/93‑1 (MQ=255)
                                                                                   ggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa      <  2:355538/93‑1 (MQ=255)
                                                                                     tgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGc    <  1:9478592/93‑1 (MQ=255)
                                                                                       aGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTg  >  1:4247266/1‑93 (MQ=255)
                                                                                       aGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTg  >  1:5820584/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTG  >  NC_000913/3934186‑3934365

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