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.6% 880.3 / ‑5.5 243H125Y (CAC→TAC) rbsBD‑ribose periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (120/122);  total (120/123)
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.50e-01

Only 100 of 260 total aligned reads displayed.

GCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCCACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTG  >  NC_000913/3934186‑3934365
                                                                                          |                                                                                         
gCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCtac                                                                                         >  2:3269043/1‑93 (MQ=255)
gCAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCtac                                                                                         <  2:2073761/93‑1 (MQ=255)
 cAGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCtaca                                                                                        <  2:3630729/93‑1 (MQ=255)
  aGTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAt                                                                                       <  1:850338/93‑1 (MQ=255)
   gTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAtt                                                                                      <  1:110667/93‑1 (MQ=255)
   gTGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACAtt                                                                                      <  1:1590322/93‑1 (MQ=255)
    tGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTg                                                                                     >  2:743694/1‑93 (MQ=255)
    tGGGTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCACCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTg                                                                                     >  1:2547691/1‑93 (MQ=255)
     gggTAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGc                                                                                    <  1:3363834/93‑1 (MQ=255)
        tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTc                                                                                 <  2:3781001/93‑1 (MQ=255)
        tAATGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTc                                                                                 <  1:1908482/93‑1 (MQ=255)
         aaTGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCt                                                                                >  1:4185464/1‑93 (MQ=255)
          aTGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTg                                                                               <  1:650174/93‑1 (MQ=255)
          aTGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTg                                                                               >  2:567256/1‑93 (MQ=255)
           tGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGa                                                                              <  2:4184752/93‑1 (MQ=255)
           tGCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGa                                                                              >  2:1468103/1‑93 (MQ=255)
            gCTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGAt                                                                             <  1:3222856/93‑1 (MQ=255)
             cTGTGAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATa                                                                            <  1:2730536/93‑1 (MQ=255)
              tgtgAAGATGGCTAACCTGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATaa                                                                           <  1:3768960/93‑1 (MQ=255)
              tgtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATaa                                                                           <  2:4608630/93‑1 (MQ=255)
              tgtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATaa                                                                           >  2:67/1‑93 (MQ=255)
              tgtgAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATaa                                                                           >  2:2481783/1‑93 (MQ=255)
                 gAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGt                                                                        >  2:4039096/1‑93 (MQ=255)
                 gAAGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGt                                                                        <  2:2036989/93‑1 (MQ=255)
                   aGATGGCTAACCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTAc                                                                      >  2:3966668/1‑93 (MQ=255)
                           aaCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGggcggc                                                              >  2:2992545/1‑93 (MQ=255)
                           aaCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGggcggc                                                              >  2:436539/1‑93 (MQ=255)
                            aCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCa                                                             >  2:1243259/1‑93 (MQ=255)
                            aCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCa                                                             >  1:3038852/1‑93 (MQ=255)
                            aCCAGGCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCa                                                             >  1:4594578/1‑93 (MQ=255)
                                ggCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAAt                                                         <  1:2455929/93‑1 (MQ=255)
                                 gCGAACATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATc                                                        >  1:3107299/1‑93 (MQ=255)
                                    aaCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCt                                                     >  1:2378630/1‑93 (MQ=255)
                                    aaCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCt                                                     <  1:696752/93‑1 (MQ=255)
                                    aaCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCt                                                     <  1:4368965/93‑1 (MQ=255)
                                     aCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTg                                                    <  2:1825115/93‑1 (MQ=255)
                                     aCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTg                                                    >  2:1690802/1‑93 (MQ=255)
                                     aCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTg                                                    >  2:3824499/1‑93 (MQ=255)
                                     aCATCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTg                                                    >  2:4067927/1‑93 (MQ=255)
                                        tCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGTGGCAAAATCGCTGGTg                                                 >  2:2869201/1‑93 (MQ=255)
                                        tCCCGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTg                                                 >  1:4111625/1‑93 (MQ=255)
                                         cccGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGa                                                <  2:1791526/93‑1 (MQ=255)
                                          ccGGTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGAt                                               <  1:1230648/93‑1 (MQ=255)
                                            ggTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTa                                             >  1:3616528/1‑93 (MQ=255)
                                            ggTTATCACTCTTGACCGCCAGGCAACGAAAGGTAAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTa                                             <  2:268458/93‑1 (MQ=255)
                                             gTTATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTAc                                            >  1:2598125/1‑93 (MQ=255)
                                              ttATCACTCGTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACa                                           <  1:1415625/93‑1 (MQ=255)
                                               tATCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACAt                                          <  1:4334291/93‑1 (MQ=255)
                                                aTCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATc                                         <  1:2232080/93‑1 (MQ=255)
                                                aTCACTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATc                                         <  2:3483621/93‑1 (MQ=255)
                                                   aCTCTTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATcgcg                                      >  2:3650460/1‑93 (MQ=255)
                                                     tctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCgaa                                    >  2:4424344/1‑93 (MQ=255)
                                                     tctTGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCgaa                                    >  2:3423768/1‑93 (MQ=255)
                                                        tGACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCgaagaa                                 <  1:776137/93‑1 (MQ=255)
                                                         gACCGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGaaa                                <  1:3957639/93‑1 (MQ=255)
                                                            cGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCg                             <  2:1283303/93‑1 (MQ=255)
                                                            cGCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCg                             <  1:3325487/93‑1 (MQ=255)
                                                             gCCAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAACTCGCTGGTGATTACATCGCGAAGAAAGCgg                            >  1:1306977/1‑93 (MQ=255)
                                                               cAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCgggg                          >  1:2006066/1‑92 (MQ=255)
                                                               cAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGt                          <  2:1540784/93‑1 (MQ=255)
                                                               cAGGCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGt                          <  1:291586/93‑1 (MQ=255)
                                                                aGGCAACTAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTg                         <  2:3864003/93‑1 (MQ=255)
                                                                aGGCAACGAAAGATGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTg                         <  1:815543/93‑1 (MQ=255)
                                                                  gCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGaa                       >  1:944875/1‑93 (MQ=255)
                                                                  gCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGaa                       <  1:2805688/93‑1 (MQ=255)
                                                                  gCAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGaa                       >  1:3148982/1‑93 (MQ=255)
                                                                   cAACGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAg                      <  1:2971530/93‑1 (MQ=255)
                                                                    aaCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAgg                     >  1:3406248/1‑93 (MQ=255)
                                                                    aaCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAgg                     <  1:3918449/93‑1 (MQ=255)
                                                                    aaCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAgg                     <  2:4229565/93‑1 (MQ=255)
                                                                    aaCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAgg                     <  2:3406911/93‑1 (MQ=255)
                                                                    aaCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAgg                     >  1:485932/1‑93 (MQ=255)
                                                                     aCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAggg                    >  2:332564/1‑92 (MQ=255)
                                                                     aCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGt                    >  1:3575979/1‑93 (MQ=255)
                                                                     aCGAAAGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGt                    >  1:470467/1‑93 (MQ=255)
                                                                        aaaGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGcc                 >  1:4521813/1‑93 (MQ=255)
                                                                         aaGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCa                <  1:2844737/93‑1 (MQ=255)
                                                                          aGGTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGCGCCaa               >  2:3198541/1‑93 (MQ=255)
                                                                           ggTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCaaa              <  2:2792751/93‑1 (MQ=255)
                                                                           ggTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCaaa              <  2:3098658/93‑1 (MQ=255)
                                                                            gTGAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAg             >  2:1054986/1‑93 (MQ=255)
                                                                              gAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGtt           <  2:2336794/93‑1 (MQ=255)
                                                                              gAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGtt           >  1:1831939/1‑93 (MQ=255)
                                                                              gAAGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGtt           >  1:4059762/1‑93 (MQ=255)
                                                                              gAAGTGGTGAGCTACATTGCTTCTGATAACGGACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGAGAAGGTGCCAAAggt           >  2:528210/1‑91 (MQ=255)
                                                                               aaGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGATGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTa          <  2:704879/93‑1 (MQ=255)
                                                                                aGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTAt         <  1:2708356/93‑1 (MQ=255)
                                                                                aGTGGTGAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTAt         <  1:3846665/93‑1 (MQ=255)
                                                                                 gtggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc        <  2:952011/93‑1 (MQ=255)
                                                                                 gtggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc        <  2:2669827/93‑1 (MQ=255)
                                                                                 gtggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATc        >  1:169407/1‑93 (MQ=255)
                                                                                   ggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa      >  1:424462/1‑93 (MQ=255)
                                                                                   ggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa      <  2:796426/93‑1 (MQ=255)
                                                                                   ggtgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa      <  1:1128343/93‑1 (MQ=255)
                                                                                   ggggAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGa      <  2:1131763/90‑1 (MQ=255)
                                                                                     tgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGc    >  2:239541/1‑93 (MQ=255)
                                                                                     tgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGc    >  2:1072006/1‑93 (MQ=255)
                                                                                     tgAGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGc    <  1:2658929/93‑1 (MQ=255)
                                                                                     tgAGCGACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGc    <  2:3061468/93‑1 (MQ=255)
                                                                                       aGCTACATTGCTTCTGATAACGTACTGGGCGGCAAAATCGCTGGTGATTACATCGCGAAGAAAGCGGGTGAAGGTGCCAAAGTTATCGAGCTg  >  1:3191254/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: