Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,021,684 Δ1 bp coding (542/609 nt)|coding (1/243 nt) b4556 →
tatD →
hypothetical protein
Mg‑dependent DNase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,021,6830C.98.7% 696.0 / ‑1.4 155coding (541/609 nt)b4556hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (85/68);  total (86/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

GGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGC  >  NC_000913/4021608‑4021754
                                                                           |                                                                       
 gCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcg                                                                        <  2:2068222/76‑2 (MQ=255)
 gCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATGACCGGTTGGGTTTGCGATGAAcg                                                                        <  1:1055760/76‑2 (MQ=255)
   acacGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAAcgac                                                                      <  1:1614959/76‑4 (MQ=255)
     acGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGACTGCGCTGAACGAcgc                                                                   >  2:1808673/1‑76 (MQ=255)
      cgcgAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGAcgcg                                                                  <  2:599954/76‑1 (MQ=255)
      cgcgAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATAGCGGCATTACCGGTTGGGTTTGCGATGAACGAcgcg                                                                  <  1:1484301/76‑1 (MQ=255)
       gcgAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCgg                                                                 >  2:4079166/1‑76 (MQ=255)
         gAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGAc                                                               <  2:1592472/76‑1 (MQ=255)
         gAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGAc                                                               >  1:2260450/1‑76 (MQ=255)
           agagATGCGGGGGGGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTg                                                             <  2:76300/76‑1 (MQ=255)
           agagATGCGGGCGTGCGGGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTg                                                             <  2:5782552/76‑1 (MQ=255)
           agagATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTg                                                             >  1:491623/1‑76 (MQ=255)
             agagGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGa                                                           <  2:911153/76‑1 (MQ=255)
             agaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGa                                                           <  2:129656/76‑1 (MQ=255)
              gaTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAg                                                          <  1:1828221/76‑1 (MQ=255)
              gaTGCAGGCGGGCGTGGCGCATGGAATTTATATCGGCCTTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAg                                                          >  2:5100368/1‑76 (MQ=255)
               aTGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGc                                                         <  1:2680553/76‑1 (MQ=255)
                tGCGGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCt                                                        <  1:4135051/76‑1 (MQ=255)
                 gCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTg                                                       >  1:1687098/1‑76 (MQ=255)
                  cAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGc                                                      <  1:5814475/76‑1 (MQ=255)
                  cAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGc                                                      <  1:4292241/76‑1 (MQ=255)
                  cAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGc                                                      <  1:4259253/76‑1 (MQ=255)
                   aGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCg                                                     >  2:4168328/1‑76 (MQ=255)
                    ggCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCgg                                                    <  2:3678735/76‑1 (MQ=255)
                    ggCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCgg                                                    <  2:1359002/76‑1 (MQ=255)
                    ggCGTGCGTGGCGCATGGAATTTATATCGGCATTACACGTTGGGCTTGCGATGAACGACGCGGACTGGCGCAGCgg                                                    >  1:5733583/1‑76 (MQ=255)
                     gcgtgcgtGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGCGCTGCggg                                                   >  2:543054/1‑76 (MQ=255)
                      cgtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGa                                                  <  2:4303321/76‑1 (MQ=255)
                      cgtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGCTTTGCGATGAACGACGCGGACTGGAGCTGCGGGa                                                  >  2:1983053/1‑76 (MQ=255)
                       gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGaa                                                 >  1:71338/1‑76 (MQ=255)
                       gtgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGaa                                                 >  1:5695768/1‑76 (MQ=255)
                        tgcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAAc                                                <  1:1429803/76‑1 (MQ=255)
                         gcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACt                                               <  2:2619150/76‑1 (MQ=255)
                         gcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACt                                               >  2:4416184/1‑76 (MQ=255)
                         gcgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACt                                               >  1:1186378/1‑76 (MQ=255)
                          cgtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACtt                                              >  1:2598685/1‑76 (MQ=255)
                           gtgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACttt                                             <  2:5674792/76‑1 (MQ=255)
                            tgGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACtttt                                            >  1:3293362/1‑76 (MQ=255)
                              gcggATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGc                                          <  2:1906666/76‑1 (MQ=255)
                              gcgcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGCGCTGCGGGAACtttttc                                          >  2:5177094/1‑74 (MQ=255)
                              gcgcATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGc                                          >  1:4588587/1‑76 (MQ=255)
                              gcgcAGGGAATTTAGATCGGCATTACCGGTTGGGTGTACGATGAACGACGCGGACTGGAGGTGCGGGAACATGTGc                                          >  1:5746570/1‑76 (MQ=255)
                                gcATGGAATTTAGATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCg                                        <  2:1267204/76‑1 (MQ=255)
                                 cATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGt                                       >  2:1144044/1‑76 (MQ=255)
                                  aTGGAATTTATATCGGCATTACTGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGtt                                      >  1:1392046/1‑76 (MQ=255)
                                  aTGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAATTTTTGCCGtt                                      >  2:3982131/1‑76 (MQ=255)
                                  aTGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGtt                                      >  2:2703695/1‑76 (MQ=255)
                                  aTGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGtt                                      >  1:3453565/1‑76 (MQ=255)
                                   tGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTg                                     >  1:2055326/1‑76 (MQ=255)
                                    ggAGTTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGa                                    <  2:1536671/76‑1 (MQ=255)
                                     gattttATATCGGCATGACCGGTTGGGTTTGCGATGGACGACGTGGACTGGAGCTGCGGGAACTTTTGCCGTCGAt                                   >  1:651845/4‑76 (MQ=255)
                                     gAATTTATATCGGCATTACCGGTTGGGTTTGCGCTGAACGACGCGGACTGGGGCTGCGGGAACTGTTGCAGTTgct                                   >  1:4591413/1‑74 (MQ=255)
                                      taTTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGATCTTTTGCCGTTGAtt                                  <  2:1059547/75‑1 (MQ=255)
                                       aTTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTc                                 >  1:3044664/1‑76 (MQ=255)
                                       aTTTATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTc                                 <  2:3743665/76‑1 (MQ=255)
                                        tttATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTcc                                >  2:4809133/1‑76 (MQ=255)
                                        tttATATCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTcc                                <  2:359789/76‑1 (MQ=255)
                                           atatCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACGTTTGCCGTTGATTCcggc                             <  2:1016622/76‑1 (MQ=255)
                                           atatCGGCATTACCGGCTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCCGGAACATCAGCCGTTGATTCcggc                             >  2:1054294/1‑76 (MQ=255)
                                            tatCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGCACTTTTGCCGTTGATTCcggcg                            >  1:4276036/1‑76 (MQ=255)
                                            tatCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcg                            >  2:1307773/1‑76 (MQ=255)
                                            tatCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcg                            >  1:2887661/1‑76 (MQ=255)
                                            tatCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcg                            <  1:581945/76‑1 (MQ=255)
                                             atCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcgg                           <  1:4168329/76‑1 (MQ=255)
                                             atCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcgg                           >  1:2517073/1‑76 (MQ=255)
                                             atCGGCATTACCGGGTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCcggcgg                           <  2:796182/76‑1 (MQ=255)
                                              tCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGa                          >  2:5394789/1‑76 (MQ=255)
                                              tCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGa                          >  1:4975829/1‑76 (MQ=255)
                                              tCGGCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGa                          <  1:4682575/76‑1 (MQ=255)
                                                ggCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaa                        >  1:5452759/1‑76 (MQ=255)
                                                 gCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaa                       >  1:550162/1‑76 (MQ=255)
                                                 gCATTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGaaaa                       >  2:3253888/1‑76 (MQ=255)
                                                   aTTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGGAAAAt                     >  2:700023/1‑76 (MQ=255)
                                                   aTTACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAt                     >  1:5130426/1‑76 (MQ=255)
                                                    ttACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAtt                    <  1:2758365/76‑1 (MQ=255)
                                                    ttACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAtt                    >  2:2339644/1‑76 (MQ=255)
                                                    ttACCGGTTGGGGTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAAtt                    >  1:5260189/1‑76 (MQ=255)
                                                     tACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTa                   >  1:1877507/1‑76 (MQ=255)
                                                     tACCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTa                   >  2:4375254/1‑76 (MQ=255)
                                                      tccGGGTGGGTTTGCGAGGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTAc                  <  2:1632730/75‑1 (MQ=255)
                                                      aCCGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTAc                  >  2:4940985/1‑76 (MQ=255)
                                                       ccGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACt                 <  1:2505701/76‑1 (MQ=255)
                                                        cGGTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTg                <  1:1472857/76‑1 (MQ=255)
                                                         ggTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCGGGCGGAAAAATTACTGa               <  2:1412269/76‑1 (MQ=255)
                                                         ggTTGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGa               <  2:5663970/76‑1 (MQ=255)
                                                          gtggggTGTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGAt              <  2:3341291/73‑1 (MQ=255)
                                                           ttGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATc             >  2:3709297/1‑76 (MQ=255)
                                                           ttGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATc             <  2:3075443/76‑1 (MQ=255)
                                                            tGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCg            <  2:5120673/76‑1 (MQ=255)
                                                            tGGGTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCg            <  1:4401255/76‑1 (MQ=255)
                                                              ggTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGaa          >  2:4464991/1‑76 (MQ=255)
                                                              ggTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGaa          <  2:2598685/76‑1 (MQ=255)
                                                              ggTTTGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGaa          <  1:4375255/76‑1 (MQ=255)
                                                                tttGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGCAAAATTACTGATCGAAAc        <  1:1307773/76‑1 (MQ=255)
                                                                 ttGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACt       >  1:1989559/1‑76 (MQ=255)
                                                                  tGCGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTg      >  2:4875473/1‑76 (MQ=255)
                                                                    cGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGAt    <  2:2992832/76‑1 (MQ=255)
                                                                    cGATGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGAt    >  1:4393827/1‑76 (MQ=255)
                                                                      aTGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgc  >  2:2379184/1‑76 (MQ=255)
                                                                      aTGAACGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATgc  <  2:3274234/76‑1 (MQ=255)
                                                                           |                                                                       
GGCACACGCGAAGAGATGCAGGCGTGCGTGGCGCATGGAATTTATATCGGCATTACCGGTTGGGTTTGCGATGAACCGACGCGGACTGGAGCTGCGGGAACTTTTGCCGTTGATTCCGGCGGAAAAATTACTGATCGAAACTGATGC  >  NC_000913/4021608‑4021754

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: