Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,709,561 (C)6→7 coding (915/1692 nt) eptB ← KDO phosphoethanolamine transferase, Ca(2+)‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,709,5551.C98.0% 183.7 / ‑3.3 51coding (921/1692 nt)eptBKDO phosphoethanolamine transferase, Ca(2+)‑inducible
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (20/30);  total (21/30)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

AAGCAATGTTGTCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGAAAA  >  NC_000913/3709416‑3709682
                                                                                                                                            |                                                                                                                               
aaGCAATGTTGTCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCcgcg                                                                                                                                  <  1:1263169/139‑1 (MQ=255)
  gCAATGTTGTCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCgcgc                                                                                                                                >  1:928528/1‑139 (MQ=255)
       gttgtCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGcccccc                                                                                                                           >  2:1098976/1‑139 (MQ=255)
       gttgtCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGcccccc                                                                                                                           >  2:1159495/1‑139 (MQ=255)
           tCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccctga                                                                                                                       >  1:668864/1‑135 (MQ=255)
             cGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACg                                                                                                                      >  2:384977/1‑139 (MQ=255)
                 aTCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACa                                                                                                                  <  1:497669/139‑1 (MQ=255)
                           tGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAAc                                                                                                        <  1:1454749/139‑1 (MQ=255)
                            gTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTTTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACg                                                                                                       <  1:764676/139‑1 (MQ=255)
                            gTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACg                                                                                                       <  1:1138046/139‑1 (MQ=255)
                            gTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACg                                                                                                       <  1:50127/139‑1 (MQ=255)
                                  ccACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTga                                                                                                 >  1:139214/1‑139 (MQ=255)
                                  ccACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTga                                                                                                 >  1:379774/1‑139 (MQ=255)
                                       tttCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGtt                                                                                            <  1:495046/139‑1 (MQ=255)
                                           gCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGt                                                                                        <  2:379774/139‑1 (MQ=255)
                                               tGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccctg                                                                                                                        <  2:258507/102‑4 (MQ=255)
                                               tGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccctg                                                                                                                        >  1:258507/1‑99 (MQ=255)
                                                gCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                                   <  2:1408655/139‑1 (MQ=255)
                                                gCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                                   <  1:457597/139‑1 (MQ=255)
                                                    aGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATc                                                                               >  1:151816/1‑139 (MQ=255)
                                                     gCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTGATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGAGGTATca                                                                              <  2:956882/139‑1 (MQ=255)
                                                     gCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATca                                                                              <  2:151816/139‑1 (MQ=255)
                                                         agagGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAgg                                                                          >  2:1434660/1‑139 (MQ=255)
                                                           agGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACg                                                                                                       <  2:441901/108‑1 (MQ=255)
                                                           agGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACg                                                                                                       >  1:441901/1‑108 (MQ=255)
                                                                aGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                                   <  2:905822/123‑1 (MQ=255)
                                                                aGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                                   >  1:905822/1‑123 (MQ=255)
                                                                              aCTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACg                                                                                                       <  1:275690/89‑1 (MQ=255)
                                                                              aCTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCa                                                     >  1:1203166/1‑139 (MQ=255)
                                                                              aCTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCa                                                     <  2:1203166/139‑1 (MQ=255)
                                                                              aCTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTGCCGCGCCCCCCCTGACGTACAAACATGCAACg                                                                                                       >  2:275690/1‑89 (MQ=255)
                                                                                    tCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                     >  1:614081/1‑101 (MQ=255)
                                                                                    tCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                     <  2:614081/101‑1 (MQ=255)
                                                                                          ccGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTcc                                          <  2:519306/138‑1 (MQ=255)
                                                                                          ccGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTcc                                          >  1:519306/1‑138 (MQ=255)
                                                                                            gcgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTgg                                       >  2:1209098/1‑139 (MQ=255)
                                                                                              gAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGc                                     <  2:308543/139‑1 (MQ=255)
                                                                                                    gttctgttctGTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTACCCACGGATCGCCGCCAGATTTTTCTCCTGGGCCAGttt                               <  2:167333/139‑1 (MQ=255)
                                                                                                      tctgttctTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGTACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCg                             <  2:928528/139‑1 (MQ=255)
                                                                                                           tctTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTg                        <  1:300623/139‑1 (MQ=255)
                                                                                                           tctTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTg                        <  2:139214/139‑1 (MQ=255)
                                                                                                               ttAATGTGCGCTGCGGATTATCTTCCGCGccccccct                                                                                                                         <  1:742484/37‑3 (MQ=255)
                                                                                                               ttAATGTGCGCTGCGGATTATCTTCCGCGccccccct                                                                                                                         >  2:742484/1‑35 (MQ=255)
                                                                                                                    gtgCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                                   >  2:1409762/1‑71 (MQ=255)
                                                                                                                    gtgCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                                   <  1:1409762/71‑1 (MQ=255)
                                                                                                                    gtgCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATca                                                                              <  2:1239175/76‑1 (MQ=255)
                                                                                                                    gtgCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATca                                                                              >  1:1239175/1‑76 (MQ=255)
                                                                                                                    gtgCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGc               >  1:809393/1‑139 (MQ=255)
                                                                                                                       cgcTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGt                                                                       <  2:1417232/80‑1 (MQ=255)
                                                                                                                       cgcTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGt                                                                       >  1:1417232/1‑80 (MQ=255)
                                                                                                                       cgcTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAcgccg                                                        <  1:477641/95‑1 (MQ=255)
                                                                                                                       cgcTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAcgccg                                                        >  2:477641/1‑95 (MQ=255)
                                                                                                                         cTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCa            <  2:249555/137‑1 (MQ=255)
                                                                                                                         cTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCa            >  1:249555/1‑137 (MQ=255)
                                                                                                                            cGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCaa                                                                                                         >  2:827833/1‑41 (MQ=255)
                                                                                                                            cGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCaa                                                                                                         <  1:827833/41‑1 (MQ=255)
                                                                                                                            cGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGACTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGcc       >  2:953560/1‑139 (MQ=255)
                                                                                                                            cGGATTATCTTCCGCGCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCg      >  1:182192/1‑139 (MQ=255)
                                                                                                                                 tATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGaaaa  <  2:789012/139‑1 (MQ=255)
                                                                                                                                 tATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGaaaa  <  1:224176/139‑1 (MQ=255)
                                                                                                                                            |                                                                                                                               
AAGCAATGTTGTCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCCGAAAA  >  NC_000913/3709416‑3709682

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