Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,362,984 A→G I217T (ATC→ACC)  yfaY ← inactive PncC family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,362,9840AG98.1% 173.6 / ‑4.3 53I217T (ATC→ACC) yfaYinactive PncC family protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (27/24);  total (28/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 = 7.83e-01

CACGCTCTGTCCGGCAACACGTTTAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGATAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGCTTTCCGAACGACCAAAAGTA  >  NC_000913/2362893‑2363079
                                                                                           |                                                                                               
cACGCTCTGTCCGGCAACACGTTTAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCAtt                                                                                         >  1:1704451/1‑100 (MQ=255)
 aCGCTCTGTCCGGCAACACGTTTAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTg                                                                                        <  2:583642/100‑1 (MQ=255)
       tGTCCGGCAACACGTTTAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGc                                                                                  >  1:1218139/1‑100 (MQ=255)
         tCCGGCAACACGTTTAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGa                                                                                >  1:1424534/1‑100 (MQ=255)
         tCCGGCAACACGTTTAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGa                                                                                >  1:718727/1‑100 (MQ=255)
           cGGCAACACGTTTAACATCCAGCCACAGTTTTTCGATGGCCTGTTGCTCGGTTGCCGGGCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGc                                                                                  <  1:1432527/96‑1 (MQ=255)
             ccgatctgTTTAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGc                                                                            <  2:1484584/93‑1 (MQ=255)
              cAACACGTTTAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGcc                                                                           <  1:163600/100‑1 (MQ=255)
                    gTTTAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGcagatcgg                                                                     >  1:1484573/1‑93 (MQ=255)
                    gTTTAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGt                                                                     >  1:489895/1‑100 (MQ=255)
                    gTTTAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGt                                                                     <  2:1537676/100‑1 (MQ=255)
                    gTTTAACATCCAGCCACAGTTTTTCCATCGCCGGTGGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGt                                                                     <  2:891354/100‑1 (MQ=255)
                     tttAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGtt                                                                    <  1:1027198/100‑1 (MQ=255)
                      ttAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTa                                                                   >  1:201774/1‑100 (MQ=255)
                       tAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTAc                                                                  <  1:1643119/100‑1 (MQ=255)
                        aaCATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCCTCCCGGTGCGGCGAGATTCAGATCGATGGTAGTCCTGGAGGAGCGATAGCCCATTGTTACg                                                                 >  1:1589868/1‑100 (MQ=255)
                          cATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGcc                                                               >  2:1411951/1‑100 (MQ=255)
                          cATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGcc                                                               >  2:1271272/1‑100 (MQ=255)
                            tCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCcg                                                             >  2:819836/1‑100 (MQ=255)
                            tCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCcg                                                             >  1:305322/1‑100 (MQ=255)
                             ccAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCcg                                                             >  1:1158631/1‑99 (MQ=255)
                               aGCCACATTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCcggc                                                           <  2:1168907/99‑1 (MQ=255)
                                gCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCcggcgg                                                         <  2:660794/100‑1 (MQ=255)
                                 ccACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCggcggc                                                        >  2:1563127/1‑100 (MQ=255)
                                 ccACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCggcggc                                                        <  2:1218150/100‑1 (MQ=255)
                                 ccACAGTGTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTGCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCggcggc                                                        <  2:1798797/100‑1 (MQ=255)
                                      gTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGTGGAGCGAAAGCCTATTGTTATGCCCGGCGGCAGt                                                     >  1:1062502/1‑98 (MQ=255)
                                       tttttCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGt                                                  <  2:1344236/100‑1 (MQ=255)
                                           tCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTagag                                              <  2:1424545/100‑1 (MQ=255)
                                                 gCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCa                                        >  2:1180289/1‑100 (MQ=255)
                                                        gCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTtg                                 >  1:1002198/1‑100 (MQ=255)
                                                         cTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGAGTCCAGGCTTtgt                                >  2:993770/1‑100 (MQ=255)
                                                          tCGCTTGCCTGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTtgtg                               <  2:625276/100‑1 (MQ=255)
                                                            gCTTGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGcc                             <  1:723687/100‑1 (MQ=255)
                                                              ttGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAg                           <  2:1442154/100‑1 (MQ=255)
                                                              ttGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAg                           <  2:718731/100‑1 (MQ=255)
                                                               tGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGGCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGa                          >  1:437181/1‑100 (MQ=255)
                                                               tGCCGGTCCGGTGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCTCGACGGCAGTTGTAGAATGTCCAGGCTTTGTGCCAGa                          >  2:1498470/1‑100 (MQ=255)
                                                                 ccggtccggtGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATc                        >  1:1363078/1‑100 (MQ=255)
                                                                   ggtccggtGAGTTTCAGTTCGATGGTAGGCATTGATGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGc                      <  2:497321/100‑1 (MQ=255)
                                                                    gtccggtGAGTTTCATTTCGATTTTATGCATTGAGTATCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGc                      <  2:419988/99‑1 (MQ=255)
                                                                      ccggtGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAATGTCCAGGCTTTGAGCAAGATCGCtct                   >  2:142342/1‑98 (MQ=255)
                                                                        tgtgtgtTTCAGTTCGATTGTAGGTATTGAGGAGCGATAGCCCATTGTAACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGCTTTcc                 <  2:1132587/99‑1 (MQ=255)
                                                                        ggtGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGCTTTcc                 <  2:608548/100‑1 (MQ=255)
                                                                        ggtGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGCTTTcc                 <  1:399432/100‑1 (MQ=255)
                                                                        ggtGAGTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGCTTTc                  <  2:295304/99‑1 (MQ=255)
                                                                             gTTTCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGCTTTCCGAACg            >  2:1086692/1‑100 (MQ=255)
                                                                             gTTTCAGTTCGATGATAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGCTTTCCGAACg            >  2:638102/1‑100 (MQ=255)
                                                                               ttCAGTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGACTTTGTGCCAGATCGCTTTCCGAACGAc          >  2:599611/1‑100 (MQ=255)
                                                                                   gTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGCTTTCCGAACGACCaaa      >  2:1460300/1‑100 (MQ=255)
                                                                                   gTTCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGCTTTCCGAACGACCaaa      >  2:1417082/1‑100 (MQ=255)
                                                                                     tCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGCTTTCCGAACGACCAAAag    <  2:656440/100‑1 (MQ=255)
                                                                                     tCGATGGTAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGCTTTCCGAACGACCAAAag    >  1:7426/1‑100 (MQ=255)
                                                                                       gatggtAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGCTTTCCGAACGACCAAAagta  <  1:1615326/100‑1 (MQ=255)
                                                                                           |                                                                                               
CACGCTCTGTCCGGCAACACGTTTAACATCCAGCCACAGTTTTTCCATCGCCTGTTGCTCGCTTGCCGGTCCGGTGAGTTTCAGTTCGATGATAGGCATTGAGGAGCGATAGCCCATTGTTACGCCCGGCGGCAGTTGTAGAGTGTCCAGGCTTTGTGCCAGATCGCTTTCCGAACGACCAAAAGTA  >  NC_000913/2362893‑2363079

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