Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,820,538 A→G H171R (CAC→CGC)  spoT → (p)ppGpp synthetase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,820,5380AG99.6% 699.2 / ‑5.7 226H171R (CAC→CGC) spoT(p)ppGpp synthetase II
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (110/114);  total (111/115)
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.96e-01

Only 100 of 240 total aligned reads displayed.

CCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAG  >  NC_000913/3820448‑3820627
                                                                                          |                                                                                         
cccACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGcc                                                                                         >  2:2459255/1‑93 (MQ=255)
 ccACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAACCGTCGCCGCCTCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCg                                                                                        >  1:1007436/1‑93 (MQ=255)
 ccACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCg                                                                                        >  2:2607081/1‑93 (MQ=255)
 ccACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCg                                                                                        >  1:3344007/1‑93 (MQ=255)
   acaacaTGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGtt                                                                                      <  2:1207477/93‑1 (MQ=255)
   acaacaTGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGtt                                                                                      >  1:2326235/1‑93 (MQ=255)
   acaacaTGCGCACGCCGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCTGTGAAACTCTCGATATTTATAGCCCGCTGGCGCGCCGtt                                                                                      >  2:3638516/1‑93 (MQ=255)
    caacaTGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGttt                                                                                     >  2:1431380/1‑93 (MQ=255)
    caacaTGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGttt                                                                                     >  2:3767797/1‑93 (MQ=255)
    caacaTGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGttt                                                                                     >  2:629901/1‑93 (MQ=255)
      acaTGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAg                                                                                   >  1:3112015/1‑93 (MQ=255)
      acaTGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAg                                                                                   >  2:1401343/1‑93 (MQ=255)
       caTGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAgg                                                                                  <  1:1708864/93‑1 (MQ=255)
         tGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTa                                                                                <  2:4153557/93‑1 (MQ=255)
         tGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCATTTAGGTa                                                                                <  2:4569850/93‑1 (MQ=255)
          gcgcACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTAt                                                                               <  1:4546577/93‑1 (MQ=255)
          gcgcACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTAt                                                                               <  2:3997965/93‑1 (MQ=255)
          gcgcACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTAt                                                                               <  1:2748874/93‑1 (MQ=255)
            gcACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGGCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATcc                                                                             <  1:4414669/93‑1 (MQ=255)
            gcACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATcc                                                                             >  1:3974089/1‑93 (MQ=255)
            gcACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATcc                                                                             >  2:4638285/1‑93 (MQ=255)
               cGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATccacc                                                                          >  2:4249508/1‑93 (MQ=255)
                ccTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATccacca                                                                         <  1:4504191/92‑1 (MQ=255)
                  tGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACcaca                                                                       >  1:1586463/1‑93 (MQ=255)
                  tGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGGGCAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACcaca                                                                       >  2:3320267/1‑93 (MQ=255)
                   gggCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACAt                                                                      <  2:2118459/93‑1 (MQ=255)
                    ggCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACAtt                                                                     >  1:2458016/1‑93 (MQ=255)
                    ggCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACAtt                                                                     <  1:873142/93‑1 (MQ=255)
                      cTCACTTCGCCCGGACAGACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTaa                                                                   >  1:4173176/1‑93 (MQ=255)
                      cTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTaa                                                                   <  1:1228905/93‑1 (MQ=255)
                      cTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTaa                                                                   <  1:3440854/93‑1 (MQ=255)
                      cTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTaa                                                                   <  1:3005487/93‑1 (MQ=255)
                      cTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTaa                                                                   <  2:1053992/93‑1 (MQ=255)
                      cTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTaa                                                                   >  2:3213694/1‑93 (MQ=255)
                         aCTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAAc                                                                <  1:574744/93‑1 (MQ=255)
                         aCTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAAc                                                                <  2:4332908/93‑1 (MQ=255)
                          cTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAAcc                                                               <  1:147922/93‑1 (MQ=255)
                          cTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAAcc                                                               >  1:2074633/1‑93 (MQ=255)
                          cTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAAcc                                                               <  2:2267344/93‑1 (MQ=255)
                          cTTCGCCCAGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAAcc                                                               <  2:93190/93‑1 (MQ=255)
                           ttCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCCCCACATTAAAACCg                                                              >  2:3023550/1‑93 (MQ=255)
                           ttCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTGAAACCg                                                              >  2:4359578/1‑93 (MQ=255)
                           ttCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCg                                                              >  1:3791816/1‑93 (MQ=255)
                           ttCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCg                                                              >  2:1296323/1‑93 (MQ=255)
                            tCGCCCGGACAAACGTCGCCTCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGa                                                             <  2:1918655/93‑1 (MQ=255)
                              gCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAAc                                                           <  2:3173259/93‑1 (MQ=255)
                              gCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAAc                                                           <  2:2751555/93‑1 (MQ=255)
                                ccGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTc                                                         <  2:873339/93‑1 (MQ=255)
                                  ggACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGCACTCGa                                                       >  1:3233453/1‑93 (MQ=255)
                                  ggACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGa                                                       <  1:3660577/93‑1 (MQ=255)
                                  ggACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGAGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGa                                                       <  2:3141035/93‑1 (MQ=255)
                                   gACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGaa                                                      <  2:1668844/93‑1 (MQ=255)
                                   gACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGaa                                                      >  1:2189164/1‑93 (MQ=255)
                                   gACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGaa                                                      >  1:3630698/1‑93 (MQ=255)
                                    aCAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAag                                                     <  1:2890280/93‑1 (MQ=255)
                                     cAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAaga                                                    <  1:1868102/93‑1 (MQ=255)
                                     cAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAaga                                                    <  2:2715746/93‑1 (MQ=255)
                                      aaccGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTGAGGTATCCCCCACATTAAAACCGAACGCGAagag                                                   <  2:2850186/90‑1 (MQ=255)
                                       aaCGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGc                                                  >  1:664526/1‑93 (MQ=255)
                                          gTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTgg                                               <  2:2217270/93‑1 (MQ=255)
                                          gTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTgg                                               <  1:3573683/93‑1 (MQ=255)
                                           tCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTggg                                              <  2:738806/93‑1 (MQ=255)
                                           tCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTggg                                              <  2:3516363/93‑1 (MQ=255)
                                           tCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTggg                                              >  1:114136/1‑93 (MQ=255)
                                           tCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTggg                                              >  1:2356025/1‑93 (MQ=255)
                                               cgcATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGtttt                                          >  1:1624898/1‑93 (MQ=255)
                                                 cATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGa                                        <  2:2080138/93‑1 (MQ=255)
                                                  aTCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGat                                       >  1:4430609/1‑92 (MQ=255)
                                                  aTCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAg                                       <  1:2706102/93‑1 (MQ=255)
                                                  aTCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAg                                       <  1:2306899/93‑1 (MQ=255)
                                                    cGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGgc                                     >  2:3944412/1‑93 (MQ=255)
                                                      cccGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGgcgc                                   <  1:527490/93‑1 (MQ=255)
                                                       ccGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCt                                  >  2:3819878/1‑93 (MQ=255)
                                                       ccGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCt                                  >  1:2654694/1‑93 (MQ=255)
                                                       ccGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCt                                  <  1:4469832/93‑1 (MQ=255)
                                                       ccGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCt                                  >  1:1495633/1‑93 (MQ=255)
                                                        cGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTg                                 >  2:3581843/1‑93 (MQ=255)
                                                          tGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTa                               >  1:4649728/1‑93 (MQ=255)
                                                          tGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTa                               <  1:799376/93‑1 (MQ=255)
                                                          tGAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTa                               >  2:1790744/1‑93 (MQ=255)
                                                           gAAACTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTAt                              >  1:1611822/1‑93 (MQ=255)
                                                            aaaCTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTCAGGCGCTGTATc                             <  1:2533528/93‑1 (MQ=255)
                                                             aaCTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATcc                            >  1:4326235/1‑93 (MQ=255)
                                                             aaCTCTCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATcc                            >  2:2418271/1‑93 (MQ=255)
                                                               ctctCGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCa                          >  2:223745/1‑93 (MQ=255)
                                                                tctcGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCaa                         <  2:1856261/93‑1 (MQ=255)
                                                                 ctcGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAAc                        <  1:3924359/93‑1 (MQ=255)
                                                                  tcGAAATTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAAcc                       <  1:3709597/93‑1 (MQ=255)
                                                                     aaaTTTATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGtt                    >  1:4052369/1‑93 (MQ=255)
                                                                        tttATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATc                 >  1:2866368/1‑93 (MQ=255)
                                                                        tttATAGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATc                 >  2:4487486/1‑93 (MQ=255)
                                                                           ataGCCCGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATcgcg              <  1:1359891/93‑1 (MQ=255)
                                                                                ccGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATc         <  2:3416097/93‑1 (MQ=255)
                                                                                ccGCTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATc         >  1:1336630/1‑93 (MQ=255)
                                                                                  ccTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCaa       <  2:1720393/92‑1 (MQ=255)
                                                                                   gtGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCaaa      <  2:3656560/92‑1 (MQ=255)
                                                                                   cTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCaaa      >  1:1830774/1‑93 (MQ=255)
                                                                                   cTGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCaaa      >  2:3091571/1‑93 (MQ=255)
                                                                                    tGGCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAaag     >  1:4260968/1‑93 (MQ=255)
                                                                                     ggCGCGCCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAaaga    >  2:1591879/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCCACAACATGCGCACGCTGGGCTCACTTCGCCCGGACAAACGTCGCCGCATCGCCCGTGAAACTCTCGAAATTTATAGCCCGCTGGCGCACCGTTTAGGTATCCACCACATTAAAACCGAACTCGAAGAGCTGGGTTTTGAGGCGCTGTATCCCAACCGTTATCGCGTAATCAAAGAAG  >  NC_000913/3820448‑3820627

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