Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,931,938 T→C S179P (TCC→CCC)  rbsA → ATP‑binding component of D‑ribose high‑affinity transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,931,9380TC99.6% 1444.4 / ‑4.2 456S179P (TCC→CCC) rbsAATP‑binding component of D‑ribose high‑affinity transport system
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (205/249);  total (206/250)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

Only 100 of 487 total aligned reads displayed.

ATGGTTGAAATCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAATCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGATTTGCGA  >  NC_000913/3931848‑3932029
                                                                                          |                                                                                           
aTGGTTGAAATCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAAccc                                                                                           >  1:861432/1‑93 (MQ=255)
   gTTGAAATCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTg                                                                                        <  1:6149681/93‑1 (MQ=255)
   gTTGAAATCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTg                                                                                        <  2:6207229/93‑1 (MQ=255)
    ttGAAATCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGt                                                                                       <  1:9567990/93‑1 (MQ=255)
    ttGAAATCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGt                                                                                       >  2:9131330/1‑93 (MQ=255)
     tGAAATCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGtt                                                                                      <  1:4199782/93‑1 (MQ=255)
      gAAATCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTc                                                                                     <  1:3878028/93‑1 (MQ=255)
        aaTCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCcg                                                                                   >  2:6460437/1‑93 (MQ=255)
        aaTCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCcg                                                                                   >  1:4713745/1‑93 (MQ=255)
         aTCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCcgc                                                                                  <  1:8063430/93‑1 (MQ=255)
          tCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCcgcg                                                                                 >  1:200592/1‑93 (MQ=255)
          tCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCcgcg                                                                                 <  2:8479619/93‑1 (MQ=255)
             ccAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCa                                                                              >  2:5648532/1‑93 (MQ=255)
             ccAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCa                                                                              <  1:7756066/93‑1 (MQ=255)
              cAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCAt                                                                             >  1:7101428/1‑93 (MQ=255)
              cAAAGTCCTGACCTTTGAGTCGAAAGTCATCATTATGGATCAACCGCCCGATGCGCTGACCAATACCGAAACCGAACCCCTGATCCGCGTCAt                                                                             <  2:147701/93‑1 (MQ=255)
               aaaGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATc                                                                            >  2:3659327/1‑93 (MQ=255)
                aaGGGCTGAGCTTTGAGTCGAAAGTCATCATTATGGGTGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATcc                                                                           <  1:6516467/93‑1 (MQ=255)
                   tGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCcgcg                                                                        <  2:9747554/93‑1 (MQ=255)
                    gCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGa                                                                       >  2:1149035/1‑93 (MQ=255)
                     cTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAg                                                                      <  2:6094935/93‑1 (MQ=255)
                      tGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGc                                                                     <  2:1728578/93‑1 (MQ=255)
                      tGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGc                                                                     >  1:3207965/1‑93 (MQ=255)
                       gAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCt                                                                    <  1:2530994/93‑1 (MQ=255)
                         gCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGa                                                                  <  1:4371847/93‑1 (MQ=255)
                            ttGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAACCCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAAt                                                               <  2:4450944/93‑1 (MQ=255)
                            ttGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAAt                                                               <  2:3774521/93‑1 (MQ=255)
                             tGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATc                                                              >  1:1027653/1‑93 (MQ=255)
                                gTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCa                                                           >  2:7824678/1‑93 (MQ=255)
                                 tCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGCCCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCaa                                                          <  2:8979001/93‑1 (MQ=255)
                                 tCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCaa                                                          <  1:6228767/93‑1 (MQ=255)
                                 tCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCaa                                                          <  2:2526567/93‑1 (MQ=255)
                                  cGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAg                                                         <  1:5390177/93‑1 (MQ=255)
                                  cGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAg                                                         <  1:3250828/93‑1 (MQ=255)
                                   gAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAgg                                                        >  2:8369867/1‑93 (MQ=255)
                                    aaaGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGc                                                       <  1:2642215/93‑1 (MQ=255)
                                    aaaGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGc                                                       <  2:4007275/93‑1 (MQ=255)
                                    aaaGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGc                                                       >  2:4205456/1‑93 (MQ=255)
                                     aaGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGcc                                                      >  2:9291820/1‑93 (MQ=255)
                                       gTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCcgc                                                    >  1:3746802/1‑93 (MQ=255)
                                        tcatcaTTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCcgcg                                                   >  2:7089916/1‑93 (MQ=255)
                                        tcatcaTTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCcgcg                                                   <  2:7368174/93‑1 (MQ=255)
                                        tcatcaTTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCcgcg                                                   >  1:3680663/1‑93 (MQ=255)
                                          atcatTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGt                                                 >  1:2980505/1‑93 (MQ=255)
                                           tcatTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGtt                                                >  1:3402129/1‑92 (MQ=255)
                                           tcatTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTa                                                >  2:5996355/1‑93 (MQ=255)
                                           tcatTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTa                                                <  1:1302010/93‑1 (MQ=255)
                                           tcatTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTa                                                >  2:4777231/1‑93 (MQ=255)
                                            catTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTAt                                               <  2:8116915/93‑1 (MQ=255)
                                             atTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTAtt                                              >  2:3345622/1‑93 (MQ=255)
                                              ttATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTg                                             >  1:6463480/1‑93 (MQ=255)
                                              ttATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTg                                             >  2:2197035/1‑93 (MQ=255)
                                                aTGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTc                                           <  1:7838823/93‑1 (MQ=255)
                                                  ggATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCta                                         <  2:3003805/93‑1 (MQ=255)
                                                  ggATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCta                                         >  1:4881015/1‑93 (MQ=255)
                                                  ggATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCta                                         >  2:4814152/1‑93 (MQ=255)
                                                   gATGACCCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCtat                                        <  2:7925472/93‑1 (MQ=255)
                                                   gATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCGGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCtat                                        <  2:2731432/93‑1 (MQ=255)
                                                   gATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCtat                                        <  1:6923265/93‑1 (MQ=255)
                                                    aTGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCtata                                       <  2:123913/93‑1 (MQ=255)
                                                     tGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTatat                                      <  2:5338030/93‑1 (MQ=255)
                                                      gAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATAtc                                     <  1:7263346/93‑1 (MQ=255)
                                                      gAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATAtc                                     >  2:500913/1‑93 (MQ=255)
                                                        accgaccgATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATAtctc                                   >  1:939696/1‑93 (MQ=255)
                                                         ccgaccgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTcc                                  <  1:7538409/93‑1 (MQ=255)
                                                         ccgaccgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTcc                                  >  1:8429030/1‑93 (MQ=255)
                                                         ccgaccgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGCCTATATCTcc                                  >  1:1398448/1‑93 (MQ=255)
                                                          cgaccgaTGCGCTGCCCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTccc                                 <  2:2803672/93‑1 (MQ=255)
                                                          cgaccgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGCAATCGCAAGGCCGCGGTATTGTCTATATCTccc                                 >  1:8899149/1‑93 (MQ=255)
                                                           gaccgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCa                                <  2:1870145/93‑1 (MQ=255)
                                                           gaccgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCa                                >  1:9870571/1‑93 (MQ=255)
                                                           gaccgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCa                                <  1:3997592/93‑1 (MQ=255)
                                                           gaccgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCa                                <  2:7644736/93‑1 (MQ=255)
                                                            accgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCAc                               >  1:2756178/1‑93 (MQ=255)
                                                              cgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCg                             >  1:5822512/1‑93 (MQ=255)
                                                               gaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGc                            >  2:5949715/1‑93 (MQ=255)
                                                               gaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCCGTATTGTCTATATCTCCCACCGc                            >  2:5600588/1‑93 (MQ=255)
                                                                 tGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCAt                          <  2:20517/93‑1 (MQ=255)
                                                                  gcgcTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATg                         <  1:6776594/93‑1 (MQ=255)
                                                                    gcTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATgaa                       >  2:9476342/1‑93 (MQ=255)
                                                                     cTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATgaaa                      <  1:7590400/93‑1 (MQ=255)
                                                                     cTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATgaaa                      <  1:7496624/93‑1 (MQ=255)
                                                                      tGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATgaaag                     <  1:8257152/93‑1 (MQ=255)
                                                                         ccGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATgaaagaaa                  >  1:5544672/1‑93 (MQ=255)
                                                                          cGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAAt                 <  1:9142833/93‑1 (MQ=255)
                                                                           gATCCCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATc                <  2:3517987/93‑1 (MQ=255)
                                                                           gATACCTATACTGACCCCCTGTTCCGCGTCATCCGCCGGCTGAAATCGCAAGGCGGCGTATTGTCTATATCTCCCACCGCATGAAAGAAATCt               <  1:4541743/93‑1 (MQ=255)
                                                                           gATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATc                >  1:1156771/1‑93 (MQ=255)
                                                                           gATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATc                >  1:10005869/1‑93 (MQ=255)
                                                                            aTACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCt               >  1:8754076/1‑93 (MQ=255)
                                                                              accgaacccgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTc             <  2:3207337/93‑1 (MQ=255)
                                                                              accgaaaccgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTc             <  2:240145/93‑1 (MQ=255)
                                                                               ccgaaaccgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCg            <  1:357806/93‑1 (MQ=255)
                                                                               ccgaaaccgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCg            <  2:8837490/93‑1 (MQ=255)
                                                                               ccgaaaccgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCg            <  1:5094451/93‑1 (MQ=255)
                                                                                cgaaaccgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCga           <  1:8820403/93‑1 (MQ=255)
                                                                                 gaaaccgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCgag          >  2:6885180/1‑93 (MQ=255)
                                                                                    accgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGAtt       <  2:8214043/93‑1 (MQ=255)
                                                                                    accgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGAtt       >  1:4940310/1‑93 (MQ=255)
                                                                                      cgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGATTTg     <  1:2212461/93‑1 (MQ=255)
                                                                                          |                                                                                           
ATGGTTGAAATCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAATCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGATTTGCGA  >  NC_000913/3931848‑3932029

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