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,2760CT98.6% 506.6 / ‑3.9 145H125Y (CAC→TAC) rbsBD‑ribose periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (62/81);  total (63/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 150 total aligned reads displayed.

GCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCT  >  NC_000913/3934186‑3934364
                                                                                          |                                                                                        
gCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCtac                                                                                        >  1:3258865/1‑93 (MQ=255)
gCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCtac                                                                                        <  1:2842466/93‑1 (MQ=255)
  aGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAt                                                                                      >  2:837736/1‑93 (MQ=255)
  aGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAt                                                                                      >  1:708453/1‑93 (MQ=255)
  aGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAt                                                                                      >  1:3041213/1‑93 (MQ=255)
   gTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAtt                                                                                     >  1:829275/1‑93 (MQ=255)
   gTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTATTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAtt                                                                                     >  2:599117/1‑93 (MQ=255)
    tGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTg                                                                                    <  1:164494/93‑1 (MQ=255)
     gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGc                                                                                   >  2:2206283/1‑93 (MQ=255)
      ggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCt                                                                                  <  2:1884488/93‑1 (MQ=255)
       gTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCtt                                                                                 <  1:4261755/93‑1 (MQ=255)
       gTAATGCTGTGAAGAGAGCTAACGAGGCTAACAGACCGGTTGTCAATCGTGACCTCCAGGCACCGCAAAGTGAAGTCGTGAGCGACATTGGtt                                                                                 <  2:3991123/93‑1 (MQ=255)
        tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTc                                                                                <  1:3761112/93‑1 (MQ=255)
        tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTc                                                                                <  1:2247135/93‑1 (MQ=255)
         aaTGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCt                                                                               >  2:2385138/1‑93 (MQ=255)
            gCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGAt                                                                            >  1:2560946/1‑93 (MQ=255)
              tgtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATaa                                                                          >  2:824466/1‑93 (MQ=255)
              tgtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATaa                                                                          >  2:294551/1‑93 (MQ=255)
                tgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACg                                                                        <  2:1443620/93‑1 (MQ=255)
                 gAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGt                                                                       <  2:3712291/93‑1 (MQ=255)
                 gAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGt                                                                       >  1:1706802/1‑93 (MQ=255)
                  aaGATGGCTAACCAGGCGACCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTa                                                                      <  2:4267900/93‑1 (MQ=255)
                   aGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTAc                                                                     <  2:3558933/93‑1 (MQ=255)
                    gATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACt                                                                    <  2:3023204/93‑1 (MQ=255)
                    gATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACt                                                                    <  1:3656946/93‑1 (MQ=255)
                    gATGGCTAACCAGGCGAACATACCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACt                                                                    >  2:1262284/1‑93 (MQ=255)
                     aTGGCTAACCAGGCGATCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTg                                                                   <  1:4225009/93‑1 (MQ=255)
                     aTGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTg                                                                   <  1:3247278/93‑1 (MQ=255)
                       ggCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTggg                                                                 <  2:2276407/93‑1 (MQ=255)
                           aaCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGggcggc                                                             <  2:2173747/93‑1 (MQ=255)
                           aaCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACCTCGCTTCTGATAACGAACTGggaggc                                                             >  2:1585476/1‑93 (MQ=255)
                           aaCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGGGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGggcggc                                                             >  1:2477209/1‑93 (MQ=255)
                             ccAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCaa                                                           <  1:1807878/93‑1 (MQ=255)
                               aGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCaaaa                                                         >  2:984250/1‑93 (MQ=255)
                                ggCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAAt                                                        <  1:1001822/93‑1 (MQ=255)
                                 gCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATc                                                       <  2:2165986/93‑1 (MQ=255)
                                   gAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGc                                                     >  1:2683367/1‑93 (MQ=255)
                                   gAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGc                                                     <  1:2832303/93‑1 (MQ=255)
                                   gAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGc                                                     <  1:90371/93‑1 (MQ=255)
                                     aCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTg                                                   >  1:2608492/1‑93 (MQ=255)
                                     aCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTg                                                   >  1:4412531/1‑93 (MQ=255)
                                       aTCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTggg                                                 >  1:3332890/1‑92 (MQ=255)
                                       aTCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGt                                                 >  1:4387619/1‑93 (MQ=255)
                                        tCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTg                                                >  1:381985/1‑93 (MQ=255)
                                        tCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTg                                                <  1:4452580/93‑1 (MQ=255)
                                         cccGGTTATCACTCTTGACCGCCAGGCAACGAACGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGa                                               <  1:4593531/93‑1 (MQ=255)
                                          ccGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGAt                                              >  1:3005209/1‑93 (MQ=255)
                                           cGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGAtt                                             >  1:2301592/1‑93 (MQ=255)
                                             gTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTAc                                           <  2:1578296/93‑1 (MQ=255)
                                             gTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTAc                                           >  1:1813238/1‑93 (MQ=255)
                                             gTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGGGATTaa                                           >  2:3217415/1‑92 (MQ=255)
                                               tATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACAt                                         <  2:2903410/93‑1 (MQ=255)
                                               tATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACAt                                         <  1:3575816/93‑1 (MQ=255)
                                                aTCCCTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATc                                        <  2:3700258/93‑1 (MQ=255)
                                                aTCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATc                                        <  2:1963662/93‑1 (MQ=255)
                                                  cACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcgc                                      <  2:265998/93‑1 (MQ=255)
                                                  cACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcgc                                      >  1:3823762/1‑93 (MQ=255)
                                                   aCTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcgcg                                     >  2:244289/1‑93 (MQ=255)
                                                    ctctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCga                                    >  1:1006575/1‑93 (MQ=255)
                                                     tctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCgaa                                   >  1:4652076/1‑93 (MQ=255)
                                                     tctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCgaa                                   <  1:2676411/93‑1 (MQ=255)
                                                         gACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGaaa                               <  1:3403402/93‑1 (MQ=255)
                                                         gACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGaaa                               <  1:2144565/93‑1 (MQ=255)
                                                         gACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGaaa                               >  2:4577453/1‑93 (MQ=255)
                                                           ccGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGc                             <  2:2019650/93‑1 (MQ=255)
                                                           ccGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGc                             <  2:3857459/93‑1 (MQ=255)
                                                            cGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCg                            <  1:2047405/93‑1 (MQ=255)
                                                            cGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCg                            >  1:3142335/1‑93 (MQ=255)
                                                             gCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCgg                           <  1:4481404/93‑1 (MQ=255)
                                                             gCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCgg                           <  1:4680939/93‑1 (MQ=255)
                                                              ccAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCggg                          <  1:3759027/93‑1 (MQ=255)
                                                              ccAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCggg                          <  2:1422279/93‑1 (MQ=255)
                                                               cAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCgggg                         >  1:563658/1‑92 (MQ=255)
                                                               cAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCgggg                         >  1:313542/1‑92 (MQ=255)
                                                               cAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGt                         <  2:1294424/93‑1 (MQ=255)
                                                               cAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGAGCGGCAAAATCGCTGGTGATTACATCGCGCAGAAAGCGGGt                         >  2:2050977/1‑93 (MQ=255)
                                                                 ggaaaCGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGa                       <  2:3836057/90‑1 (MQ=255)
                                                                 ggCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGa                       <  2:1520917/93‑1 (MQ=255)
                                                                 ggCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGa                       >  1:4474690/1‑93 (MQ=255)
                                                                    aaCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAgg                    <  1:2246772/93‑1 (MQ=255)
                                                                       gAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGc                 >  2:4301605/1‑93 (MQ=255)
                                                                       gAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGc                 >  2:980786/1‑93 (MQ=255)
                                                                        aaaGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGcc                <  2:2739910/93‑1 (MQ=255)
                                                                        aaaGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGcc                <  2:1690608/93‑1 (MQ=255)
                                                                           ggTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCAGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCaaa             <  1:1335566/93‑1 (MQ=255)
                                                                            gTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAg            <  1:3659937/93‑1 (MQ=255)
                                                                             tGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGt           >  1:1433216/1‑93 (MQ=255)
                                                                             tGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGt           >  2:2304976/1‑93 (MQ=255)
                                                                              gAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGtt          <  2:3955652/93‑1 (MQ=255)
                                                                               aaGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTa         >  1:3572610/1‑93 (MQ=255)
                                                                               aaGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTa         >  2:1806155/1‑93 (MQ=255)
                                                                                aGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTAt        <  2:4526523/93‑1 (MQ=255)
                                                                                aGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTAt        <  2:703987/93‑1 (MQ=255)
                                                                                 gtggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc       <  1:234469/93‑1 (MQ=255)
                                                                                   ggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa     <  1:2409566/93‑1 (MQ=255)
                                                                                   ggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa     <  2:1496214/93‑1 (MQ=255)
                                                                                    gtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAGGAAAGCGGGTGAAGGTGCCAAAGTTATCGAg    <  2:624536/93‑1 (MQ=255)
                                                                                     tgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGc   <  2:4092630/93‑1 (MQ=255)
                                                                                      gAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCt  >  2:321806/1‑93 (MQ=255)
                                                                                      gAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCt  >  2:782277/1‑93 (MQ=255)
                                                                                          |                                                                                        
GCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCT  >  NC_000913/3934186‑3934364

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