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.5% 601.7 / ‑4.8 189S179P (TCC→CCC) rbsAATP‑binding component of D‑ribose high‑affinity transport system
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (83/105);  total (83/106)
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.54e-01

Only 100 of 204 total aligned reads displayed.

ATGGTTGAAATCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAATCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGATTTGC  >  NC_000913/3931848‑3932027
                                                                                          |                                                                                         
 tGGTTGAAATCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAAcccc                                                                                        <  2:4500186/93‑1 (MQ=255)
       aaaTCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTcc                                                                                  >  2:4011746/1‑93 (MQ=255)
         aTCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCcgc                                                                                >  1:2108157/1‑93 (MQ=255)
            gCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTc                                                                             >  2:2240582/1‑93 (MQ=255)
            gCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTc                                                                             >  1:3335821/1‑93 (MQ=255)
             ccAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCa                                                                            <  2:4749160/93‑1 (MQ=255)
              cAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCAt                                                                           <  1:1830300/93‑1 (MQ=255)
               aaaGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATc                                                                          <  2:1400374/93‑1 (MQ=255)
               aaaGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATc                                                                          >  1:1195171/1‑93 (MQ=255)
                aaGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATcc                                                                         <  2:220407/93‑1 (MQ=255)
                  gcgCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCcgc                                                                       >  2:2897100/3‑93 (MQ=255)
                  gTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCcgc                                                                       <  2:3437382/93‑1 (MQ=255)
                    gCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGa                                                                     >  1:29483/1‑93 (MQ=255)
                     cTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAg                                                                    >  2:2421378/1‑93 (MQ=255)
                     cTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAg                                                                    >  2:4564509/1‑93 (MQ=255)
                     cTGAGCTTTGAGTCGAAACTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAg                                                                    <  1:2666698/93‑1 (MQ=255)
                      tGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGc                                                                   <  2:2994099/93‑1 (MQ=255)
                       gAGCTTTGAGTCGAAAGTCATCATTATGGATTAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCt                                                                  >  1:2194028/1‑93 (MQ=255)
                        aGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGCCCGATACCGAAACCGAACCCCTGTCCCGCGTGATCCgcgcgcgg                                                                 >  1:878092/1‑91 (MQ=255)
                        aGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTg                                                                 <  2:3668952/93‑1 (MQ=255)
                         gCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGa                                                                >  1:369246/1‑93 (MQ=255)
                         gCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGa                                                                >  2:566951/1‑93 (MQ=255)
                         gCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGa                                                                >  1:2674527/1‑93 (MQ=255)
                          cTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGaa                                                               <  1:3962928/93‑1 (MQ=255)
                           tttGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGaaa                                                              <  2:763322/93‑1 (MQ=255)
                           tttGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGaaa                                                              >  2:2222870/1‑93 (MQ=255)
                           tttGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGaaa                                                              >  1:2017521/1‑93 (MQ=255)
                           tttGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGaaa                                                              >  1:1908068/1‑93 (MQ=255)
                            ttGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAAt                                                             <  2:4097816/93‑1 (MQ=255)
                             tGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATc                                                            >  2:3502088/1‑93 (MQ=255)
                             tGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATc                                                            <  2:1189711/93‑1 (MQ=255)
                                 tCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCaa                                                        >  2:3153762/1‑93 (MQ=255)
                                 tCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCaa                                                        >  1:4498176/1‑93 (MQ=255)
                                    aaaGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGc                                                     <  2:1697226/93‑1 (MQ=255)
                                    aaaGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGc                                                     >  1:3902530/1‑93 (MQ=255)
                                     aaGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCATGGcc                                                    <  2:2126343/93‑1 (MQ=255)
                                     aaGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGcc                                                    <  2:4250860/93‑1 (MQ=255)
                                     aaGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGcc                                                    <  1:3768255/93‑1 (MQ=255)
                                      aGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTACTCCGCGAGCTGAAATCGCCAGGCcg                                                   >  1:4430318/1‑93 (MQ=255)
                                       gTCATCATGATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCcgc                                                  <  1:2917532/93‑1 (MQ=255)
                                          atcatTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGt                                               >  1:53868/1‑93 (MQ=255)
                                            catTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTAt                                             <  1:10907/93‑1 (MQ=255)
                                            catTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTAt                                             <  1:655794/93‑1 (MQ=255)
                                            catTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTAt                                             <  1:4320897/93‑1 (MQ=255)
                                             atTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTAtt                                            <  1:1998559/93‑1 (MQ=255)
                                             atTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTAtt                                            <  1:2370765/93‑1 (MQ=255)
                                              ttATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTg                                           >  1:1437334/1‑93 (MQ=255)
                                              ttATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTg                                           >  1:1524961/1‑93 (MQ=255)
                                                 tGGATGAGCCGCCCGATGCGCTGACCGATACCGAATCCGACCCCCTGTCCCTCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCt                                        <  2:3268975/93‑1 (MQ=255)
                                                  ggaggaCCCGAACGATGCGCTGACTGATACCAAATCCGACCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTAATGTCta                                       <  2:604013/93‑1 (MQ=255)
                                                  ggATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCta                                       >  1:3318277/1‑93 (MQ=255)
                                                  ggATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCta                                       <  1:1176320/93‑1 (MQ=255)
                                                  ggATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCta                                       <  1:490294/93‑1 (MQ=255)
                                                  ggATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCATGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCta                                       >  1:3208352/1‑93 (MQ=255)
                                                   gATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCtat                                      >  2:2732892/1‑93 (MQ=255)
                                                   gATGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCtat                                      >  1:4510757/1‑93 (MQ=255)
                                                    gtGAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCtata                                     <  1:3606024/92‑1 (MQ=255)
                                                     tgcaccgactgATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAGTCTCAAGGCTGCGGTATTGTTTatat                                    <  2:4649632/90‑1 (MQ=255)
                                                      gAACCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATAtc                                   <  2:1055879/93‑1 (MQ=255)
                                                       aaCCGACCGATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGCCCGCGGTATTGTCTATAtct                                  >  1:1087301/1‑93 (MQ=255)
                                                        accgaccgATGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATAtctc                                 <  2:4462557/93‑1 (MQ=255)
                                                         ccgaccgaTGCGCTGACCGATACCGAAACCGACCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTcc                                <  2:3056331/93‑1 (MQ=255)
                                                         ccgaccgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTcc                                <  1:405290/93‑1 (MQ=255)
                                                          cgaccgaTGCGCTGACCTATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTccc                               <  1:4124975/93‑1 (MQ=255)
                                                           gaccgaTGCGCTGACCGATACCGAACCCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCa                              <  2:4346833/93‑1 (MQ=255)
                                                           gaccgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCa                              >  1:2297283/1‑93 (MQ=255)
                                                            accgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCAc                             >  2:1601607/1‑93 (MQ=255)
                                                            accgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCAc                             >  1:1412856/1‑93 (MQ=255)
                                                              cgaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCg                           >  1:3393947/1‑93 (MQ=255)
                                                               gaTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGc                          >  1:2070139/1‑93 (MQ=255)
                                                                aTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCa                         <  2:1142589/93‑1 (MQ=255)
                                                                aTGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCa                         <  2:3415381/93‑1 (MQ=255)
                                                                 tGCGCTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCAt                        <  1:427692/93‑1 (MQ=255)
                                                                  gcgcTGACCGATACCGAAACCGAACCCCTGTGCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATg                       <  1:807768/93‑1 (MQ=255)
                                                                   cgcTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATga                      <  2:3661489/93‑1 (MQ=255)
                                                                   cgcTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATga                      <  1:241931/93‑1 (MQ=255)
                                                                     cggACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATgaaa                    <  1:3471100/91‑1 (MQ=255)
                                                                     cTGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATgaaa                    >  1:1298523/1‑93 (MQ=255)
                                                                      tGACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATgaaag                   >  1:272065/1‑93 (MQ=255)
                                                                       gACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATgaaaga                  >  1:4633561/1‑93 (MQ=255)
                                                                       gACCGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATgaaaga                  <  1:364296/93‑1 (MQ=255)
                                                                         ccGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATgaaagaaa                <  2:1741738/93‑1 (MQ=255)
                                                                         ccGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATgaaagaaa                >  2:4255258/1‑93 (MQ=255)
                                                                          cGATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCGGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAAt               <  2:2833441/93‑1 (MQ=255)
                                                                           gATACCGCAACCGAACCCCTGTCCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATc              <  2:3848013/93‑1 (MQ=255)
                                                                           gATACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATc              <  1:2833925/93‑1 (MQ=255)
                                                                            aTACCGAAACCGACCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCt             <  1:1774707/93‑1 (MQ=255)
                                                                             tACCGAAACCGAACCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCtt            >  1:756191/1‑93 (MQ=255)
                                                                              accgaaaccgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTc           <  2:4177646/93‑1 (MQ=255)
                                                                               ccgaaaccgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCg          <  2:294282/93‑1 (MQ=255)
                                                                               ccgaaaccgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCg          <  2:2647522/93‑1 (MQ=255)
                                                                                cgaaaccgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCga         <  2:849161/93‑1 (MQ=255)
                                                                                  aaaccgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCgaga       <  2:3974494/93‑1 (MQ=255)
                                                                                    accgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGAtt     >  1:1370052/1‑93 (MQ=255)
                                                                                    accgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGAtt     <  1:2480623/93‑1 (MQ=255)
                                                                                    accgaaCCCCTGTCCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGAtt     <  2:2619030/93‑1 (MQ=255)
                                                                                     ccgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGAttt    <  1:920629/93‑1 (MQ=255)
                                                                                      cgaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGATTTg   <  2:3745696/93‑1 (MQ=255)
                                                                                       gaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGATTTGc  <  1:2595255/93‑1 (MQ=255)
                                                                                       gaaCCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGATTTGc  >  1:3071174/1‑93 (MQ=255)
                                                                                          |                                                                                         
ATGGTTGAAATCGCCAAAGTGCTGAGCTTTGAGTCGAAAGTCATCATTATGGATGAACCGACCGATGCGCTGACCGATACCGAAACCGAATCCCTGTTCCGCGTCATCCGCGAGCTGAAATCGCAAGGCCGCGGTATTGTCTATATCTCCCACCGCATGAAAGAAATCTTCGAGATTTGC  >  NC_000913/3931848‑3932027

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