Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 107,553 (T)9→8 intergenic (+79/‑152) lpxC → / → secM UDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913107,5450T.92.3% 251.9 / 15.4 65intergenic (+71/‑160)lpxC/secMUDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base . (28/32);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAC  >  NC_000913/107448‑107640
                                                                                                 |                                                                                               
aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAACTCTGGCACTCTCTCCGGCCCGGTGCACCAGTCGttt                                                                                               >  2:930878/1‑100 (MQ=255)
aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttt                                                                                               >  2:1071368/1‑100 (MQ=255)
  aGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttt                                                                                             >  1:1126343/1‑100 (MQ=255)
    cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttt                                                                                             <  1:291218/98‑1 (MQ=255)
      ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                         <  1:937759/100‑2 (MQ=255)
       cTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttga                                                                                        >  1:636703/1‑98 (MQ=255)
       cTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTg                                                                                        <  2:809146/100‑1 (MQ=255)
           aGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAttt                                                                                   >  1:524700/1‑100 (MQ=255)
            gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAtttt                                                                                  >  1:1283754/1‑100 (MQ=255)
            gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAtttt                                                                                  >  1:132847/1‑100 (MQ=255)
            gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAtttt                                                                                  <  2:130347/100‑1 (MQ=255)
            gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAtttt                                                                                  >  1:631810/1‑100 (MQ=255)
            gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAAtttn                                                                                  <  1:1441042/100‑2 (MQ=255)
               gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTtata                                                                               >  1:374563/1‑100 (MQ=255)
                  cTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAaga                                                                            <  2:694057/100‑1 (MQ=255)
                    ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGc                                                                          >  2:755336/1‑100 (MQ=255)
                    ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGACAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGc                                                                          >  2:859078/1‑100 (MQ=255)
                       aTACCGCCATTCATACTGTCGTATAAAATTGTACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGc                                                                          <  2:1486139/97‑1 (MQ=255)
                       aTAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCtat                                                                       >  2:3433/1‑100 (MQ=255)
                        tAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCgata                                                                      <  1:167361/100‑1 (MQ=255)
                        gaaCGCCAGTTAAAATGCCCTACAAAATTGCACTGCCAATTTTTGCACTTCCTCCGGCCAGGGGAACCAGTCGTCTTTTTTGAATTTTATAAGAGCta                                                                        <  2:1308125/97‑1 (MQ=255)
                          ncGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATaaa                                                                    >  2:610925/2‑100 (MQ=255)
                            gACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGCACCAGTCGTTTTTTTTGAATTTTCTAAGTGGTATacaaa                                                                  >  2:1359207/1‑100 (MQ=255)
                              cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                                <  2:945162/100‑1 (MQ=255)
                                 ttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCCCTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGACTTTTATAATAGCTATAAAAAACGGt                                                             >  2:913128/1‑100 (MQ=255)
                                      cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                        >  1:766534/1‑100 (MQ=255)
                                            tataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGt                                                  >  2:135675/1‑100 (MQ=255)
                                              taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                                >  1:159639/1‑100 (MQ=255)
                                               aaaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtttt                                               <  2:766532/100‑1 (MQ=255)
                                                aaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGtttt                                               >  2:1091712/1‑99 (MQ=255)
                                                 aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             <  1:1499981/100‑1 (MQ=255)
                                                 aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             <  1:900154/100‑1 (MQ=255)
                                                 aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                             <  1:185761/100‑1 (MQ=255)
                                                  aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                            >  1:875745/1‑100 (MQ=255)
                                                  aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                            <  2:1480100/100‑1 (MQ=255)
                                                   ttCGACTGGCAAATCTGGCTCTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           <  1:1368742/100‑1 (MQ=255)
                                                   ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                           >  1:386927/1‑100 (MQ=255)
                                                    tCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                          <  1:859083/100‑1 (MQ=255)
                                                    tCGACTGGCAAATCTGGCACTCGCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTaa                                          <  1:638796/100‑1 (MQ=255)
                                                     cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTATTTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         >  2:435989/1‑100 (MQ=255)
                                                     cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGGTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                         >  2:1046912/1‑100 (MQ=255)
                                                        cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      >  1:1449055/1‑100 (MQ=255)
                                                        cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      >  1:1435398/1‑100 (MQ=255)
                                                        cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      <  2:242054/100‑1 (MQ=255)
                                                        cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      >  2:284111/1‑100 (MQ=255)
                                                        cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                      >  1:1405801/1‑100 (MQ=255)
                                                          ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTTTTTTTTTAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  2:1459225/1‑100 (MQ=255)
                                                          ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                    >  1:489417/1‑100 (MQ=255)
                                                             aaaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTTTTTTTTTAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCTCTTTTc                                 >  2:549794/1‑100 (MQ=255)
                                                              aaTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTcc                                <  1:672248/100‑1 (MQ=255)
                                                               aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                               <  1:755338/100‑1 (MQ=255)
                                                               aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                               <  2:132847/100‑1 (MQ=255)
                                                               aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCTTTTTTTTTAAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                               <  1:328942/100‑1 (MQ=255)
                                                               aTCTGGCACTCTCTCCGGACAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                               <  2:964805/100‑1 (MQ=255)
                                                                tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                              <  1:1188549/100‑1 (MQ=255)
                                                                 cTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGca                             <  2:1283749/100‑1 (MQ=255)
                                                                  tGGTACTTTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTcc                                <  2:689035/96‑1 (MQ=255)
                                                                    gCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACaa                          >  1:241666/1‑100 (MQ=255)
                                                                     cACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACaat                         >  1:192703/1‑99 (MQ=255)
                                                                     cACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAAc                         <  1:1164149/100‑1 (MQ=255)
                                                                      aCTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACt                        >  1:104288/1‑100 (MQ=255)
                                                                           ctcCGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                   <  2:484366/100‑1 (MQ=255)
                                                                              cGGCCAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCa                >  2:882364/1‑100 (MQ=255)
                                                                                  cAGGTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCg            <  2:170170/100‑1 (MQ=255)
                                                                                     gTGAACCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc         <  1:976058/100‑1 (MQ=255)
                                                                                        aaCCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgt      <  2:270741/100‑1 (MQ=255)
                                                                                         aCCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtn     <  1:610927/100‑2 (MQ=255)
                                                                                         aCCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgtg     <  2:524699/100‑1 (MQ=255)
                                                                                         aCCAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtgt      >  2:751102/1‑99 (MQ=255)
                                                                                          ccAGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTgg    <  2:631808/100‑1 (MQ=255)
                                                                                            aGTCGTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAc  <  1:694566/100‑1 (MQ=255)
                                                                                                 |                                                                                               
AAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGAC  >  NC_000913/107448‑107640

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