Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,227,332 G→A 12.0% Y47Y (TAC→TAT pliG ← periplasmic inhibitor of g‑type lysozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,227,3320GA12.0% 145.8 / 10.8 50Y47Y (TAC→TATpliGperiplasmic inhibitor of g‑type lysozyme
Reads supporting (aligned to +/- strand):  ref base G (22/22);  new base A (3/3);  total (25/25)
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.52e-01

AGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGATTTT  >  NC_000913/1227195‑1227469
                                                                                                                                         |                                                                                                                                         
aGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTatcgtaatcgta                                                                                                                                         <  2:1435586/140‑1 (MQ=255)
aGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTatcgtaatcgta                                                                                                                                         <  2:114076/140‑1 (MQ=255)
           atatCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTAttt                                                                                                                              <  1:893872/140‑1 (MQ=255)
               cTGGATAATTCAACTGAATCATCGATTCGTGGCCCGAAGAGGTAAGTATCAGCGCGTTCATTAGAGATACTCTCATGCACTTTCTGTCGTTTTTTGGCATAGTAGGTATATGTATCGTAATCGTATCCCTTTATTTCGcc                                                                                                                          >  2:441226/1‑140 (MQ=255)
                tGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCg                                                                                                                         <  1:1212213/140‑1 (MQ=255)
                    tAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGt                                                                                                                     >  2:558133/1‑140 (MQ=255)
                    tAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGt                                                                                                                     >  2:558117/1‑140 (MQ=255)
                    tAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGt                                                                                                                     >  2:557473/1‑140 (MQ=255)
                    tAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCACTTTATTTCGCCGGAGt                                                                                                                     >  2:940871/1‑140 (MQ=255)
                          aaCTGAATCATCGATTCCTGGCCCGAAGAGGTACGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCATTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAg                                                                                                               <  2:65514/140‑1 (MQ=255)
                          aaCTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTATTCGTATCCCTTTATTTCGCCGGAGTACTGAg                                                                                                               <  1:1501709/140‑1 (MQ=255)
                          aaCTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAg                                                                                                               <  1:405392/140‑1 (MQ=255)
                          aaCTGAATCATCGAGTCCTGGCCCGAAGAGGTAAGTATCAGCGCCGTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAg                                                                                                               <  2:529271/140‑1 (MQ=255)
                            cTGAATCATCGATTCCTGGCCCGCAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAgcg                                                                                                             <  2:959223/140‑1 (MQ=255)
                            cTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAgcg                                                                                                             <  1:783322/140‑1 (MQ=255)
                                    tcGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATg                                                                                                     >  2:326494/1‑140 (MQ=255)
                                                           aGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGa                                                                              <  1:634554/140‑1 (MQ=255)
                                                             tATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACa                                                                            <  1:1017188/140‑1 (MQ=255)
                                                                  gcgcCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCtt                                                                       >  2:1194735/1‑140 (MQ=255)
                                                                  gcgcCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCtt                                                                       >  2:1207804/1‑140 (MQ=255)
                                                                  gcgcCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCtt                                                                       >  2:1194796/1‑140 (MQ=255)
                                                                     ccTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTAcc                                                                    <  2:349072/140‑1 (MQ=255)
                                                                       ttCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCg                                                                  >  2:514997/1‑140 (MQ=255)
                                                                       ttCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCg                                                                  <  2:269870/140‑1 (MQ=255)
                                                                       ttCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCg                                                                  <  2:269755/140‑1 (MQ=255)
                                                                       ttCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCg                                                                  >  2:141909/1‑140 (MQ=255)
                                                                             agagATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGCATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCaa                                                            >  1:247911/1‑140 (MQ=255)
                                                                                 aTACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGATTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACAATGACATTCTTACCCGCAGCAAAGcg                                                        >  2:1291711/1‑139 (MQ=255)
                                                                                 aTACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCt                                                        >  2:129022/1‑140 (MQ=255)
                                                                                   aCTCACATGCACTTTCTGACCTTTTTTGGCATTGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCg                                                      >  2:834769/1‑140 (MQ=255)
                                                                                   aCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCg                                                      <  1:1019670/140‑1 (MQ=255)
                                                                                   aCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCg                                                      >  1:1335638/1‑140 (MQ=255)
                                                                                   aCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCg                                                      <  1:1454592/140‑1 (MQ=255)
                                                                                   aCTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCg                                                      <  1:1456675/140‑1 (MQ=255)
                                                                                   aCTCACATGCACTTTATGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCg                                                      <  2:264497/140‑1 (MQ=255)
                                                                                        caTGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAg                                                 <  2:747751/140‑1 (MQ=255)
                                                                                         aTGCACTTTCTGAGCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGtggtgg                                                >  1:480344/1‑140 (MQ=255)
                                                                                         aTGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAgg                                                >  1:1500690/1‑140 (MQ=255)
                                                                                            cACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTa                                             >  2:41145/1‑140 (MQ=255)
                                                                                                 tCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGt                                        >  1:1148007/1‑140 (MQ=255)
                                                                                                      ccTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAag                                   <  1:362545/140‑1 (MQ=255)
                                                                                                         tttttGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACACGCAGCAAAGATCTTGGAGGTTAGCAGTGCAagcag                                >  2:750997/1‑140 (MQ=255)
                                                                                                              ggCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATa                           <  1:81640/140‑1 (MQ=255)
                                                                                                                    gAAGGTATATGTATCGTAATCATATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCg                                                      >  1:673501/1‑107 (MQ=255)
                                                                                                                    gAAGGTATATGTATCGTAATCATATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCg                                                      >  1:673502/1‑107 (MQ=255)
                                                                                                                    gAAGGTATATGTATCGTAATCATATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCg                                                      <  2:673529/107‑1 (MQ=255)
                                                                                                                    gAAGGTATATGTATCGTAATCATATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCg                                                      <  2:673502/107‑1 (MQ=255)
                                                                                                                    gAAGGTATATGTATCGTAATCATATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCg                                                      <  2:673501/107‑1 (MQ=255)
                                                                                                                    gAAGGTATATGTATCGTAATCATATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCg                                                      >  1:673529/1‑107 (MQ=255)
                                                                                                                         tataTGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTg                <  2:492955/140‑1 (MQ=255)
                                                                                                                          atatGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGATATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGa               <  1:441226/140‑1 (MQ=255)
                                                                                                                               tATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGctc          >  1:1100302/1‑140 (MQ=255)
                                                                                                                                  cgtaatcgtaTCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGAGATTCTTACCCGCAGCAAAGCTCGTGGAGGTTTGCTGTGCAAGCAGCTATACAGCCTTCCTGATGGTGTTg       >  1:1442024/1‑140 (MQ=255)
                                                                                                                                  cgtaatcgtaTCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTg       >  1:1139668/1‑140 (MQ=255)
                                                                                                                                  cgtaatcgtaTCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTg       >  2:766493/1‑140 (MQ=255)
                                                                                                                                  cgtaatcgtaTCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTg       >  2:217603/1‑140 (MQ=255)
                                                                                                                                  cgtaatcgtaTCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGAAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTg       >  1:1139658/1‑140 (MQ=255)
                                                                                                                                   gtaatcgtaTCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTgt      >  1:832943/1‑139 (MQ=255)
                                                                                                                                   gtaatcgtaTCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGa      >  1:784403/1‑140 (MQ=255)
                                                                                                                                   gtaatcgtaTCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGa      >  1:1033997/1‑140 (MQ=255)
                                                                                                                                   gtaatcgtaTCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGa      >  1:392697/1‑140 (MQ=255)
                                                                                                                                       tcgtaTCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGAtttt  >  2:244147/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                         
AGTTCTGGGGAATATCTGGATAAATCAACTGAATCATCGATTCCTGGCCCGAAGAGGTAAGTATCAGCGCCTTCATTAGAGATACTCACATGCACTTTCTGACCTTTTTTGGCATAGAAGGTATATGTATCGTAATCGTATCCCTTTATTTCGCCGGAGTACTGAGCGCTGCTATGGCCTTTTCTGAACTCAACATTGACATTCTTACCCGCAGCAAAGCTCGTGGAGGTTAGCAGTGCAAGCAGCAATACAGCCTTCCTGATGCTCTTGATTTT  >  NC_000913/1227195‑1227469

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