Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,551,174 G→C 7.1% L15V (CTG→GTG)  yfeY ← RpoE‑regulated lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,551,1740GC 7.1% 192.4 / 6.4 56L15V (CTG→GTG) yfeYRpoE‑regulated lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (29/23);  new base C (2/2);  total (31/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.76e-01

CAATGGCTTGTTCTTGCAGTGGTGTGGACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACTTCATAAAAACAAGCGCGTAAACGTCG  >  NC_000913/2551036‑2551315
                                                                                                                                          |                                                                                                                                             
cAATGGCTTGTTCTTGCAGTGGTGTGGACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATgg                                                                                                                                      <  8:60577/148‑1 (MQ=255)
   tGGCTTGTTCTTGCAGTGGTGTGGACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCa                                                                                                                                    >  8:394044/1‑147 (MQ=255)
         gTTCTTGCAGTGGTGTGGACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAgcg                                                                                                                               <  5:123402/146‑1 (MQ=255)
               gCAGTGGTGTGGACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAAc                                                                                                                         <  3:288505/146‑1 (MQ=255)
                  gtggtgTGGACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTaa                                                                                                                    <  7:394044/148‑1 (MQ=255)
                    ggtgTGGACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCa                                                                                                                                             <  1:132780/121‑1 (MQ=255)
                    ggtgTGGACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCa                                                                                                                                             >  2:132780/1‑121 (MQ=255)
                    ggtgTGGACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAAc                                                                                                                  >  7:25884/1‑148 (MQ=255)
                    ggtgTGGACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCACCATCAATGGCATAGCGCATAACATTAAAc                                                                                                                  >  7:259881/1‑148 (MQ=255)
                    ggtgTGGACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCACCATCAATGGCATAGCGCATAACATTAAAc                                                                                                                  >  7:266243/1‑148 (MQ=255)
                        tgGACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAg                                                                                                              >  4:243643/1‑148 (MQ=255)
                         ggACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCa                                                                                                                                          <  1:194378/119‑1 (MQ=255)
                         ggACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCa                                                                                                                                          >  2:194378/1‑119 (MQ=255)
                               cGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGGATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTc                                                                                                       <  2:5439/148‑1 (MQ=255)
                               cGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTc                                                                                                       <  3:243643/148‑1 (MQ=255)
                                gTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTc                                                                                                       <  3:253461/147‑1 (MQ=255)
                                gTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTc                                                                                                       >  4:253461/1‑147 (MQ=255)
                                  tAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTc                                                                                                       >  3:338044/1‑145 (MQ=255)
                                  tAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTc                                                                                                       <  4:338044/145‑1 (MQ=255)
                                  tAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATg                                                                                                    >  8:85463/1‑148 (MQ=255)
                                                   gctcgctcACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTaa                                                                                                                    <  5:127281/115‑1 (MQ=255)
                                                   gctcgctcACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTaa                                                                                                                    >  6:127281/1‑115 (MQ=255)
                                                   gctcgctcACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTa                                                                                                                     >  5:108818/1‑114 (MQ=255)
                                                   gctcgctcACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTa                                                                                                                     <  6:108818/114‑1 (MQ=255)
                                                   gctcgctcACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTAt                                                                                   >  1:357032/1‑148 (MQ=255)
                                                   gctcgctcACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTAt                                                                                   >  6:93331/1‑148 (MQ=255)
                                                    ctcgctcACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTAtt                                                                                  >  7:172447/1‑148 (MQ=255)
                                                    ctcgctcACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTAt                                                                                   <  7:85463/147‑1 (MQ=255)
                                                               tggtggACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGt                                                                         <  8:25884/146‑1 (MQ=255)
                                                                ggtggACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTAt                                                                                   >  8:186027/1‑135 (MQ=255)
                                                                ggtggACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTAt                                                                                   <  7:186027/135‑1 (MQ=255)
                                                                                 aGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGACCAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCAc                                                     >  8:192868/1‑148 (MQ=255)
                                                                                    tCCACGGGTTAGCTGCAGACCAGTTAACACAAGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATAACATGTAAATTCCTTTGGATATTCAATAACTTAGCTTGGAGTGCAAATTCACCa                                                   >  1:35115/1‑147 (MQ=255)
                                                                                    tCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTgg                                                                  >  4:23233/1‑132 (MQ=255)
                                                                                    tCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTgg                                                                  <  3:23233/132‑1 (MQ=255)
                                                                                               gCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCATCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGt                                        <  3:45106/147‑1 (MQ=255)
                                                                                                   cAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTc                                   >  7:137428/1‑148 (MQ=255)
                                                                                                     gACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCtt                                 >  6:326555/1‑148 (MQ=255)
                                                                                                       ccAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCttt                                >  4:47757/1‑147 (MQ=255)
                                                                                                          gTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACt                            <  8:137428/148‑1 (MQ=255)
                                                                                                            tAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCGTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACtt                           >  4:250346/1‑147 (MQ=255)
                                                                                                            tAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACtt                           >  8:359380/1‑147 (MQ=255)
                                                                                                            tAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAAAGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACtt                           >  4:50044/1‑147 (MQ=255)
                                                                                                                 gacgacATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACTTCATaa                      >  4:154625/1‑147 (MQ=255)
                                                                                                                     acATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACTTCATAAAAAc                  >  5:339730/1‑147 (MQ=255)
                                                                                                                        tCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACTTCATAAAAACAAgc              >  4:100926/1‑148 (MQ=255)
                                                                                                                         cGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTc                                                                                                  <  7:263508/63‑1 (MQ=255)
                                                                                                                         cGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTc                                                                                                  >  8:263508/1‑63 (MQ=255)
                                                                                                                         cGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACTTCATAAAAACAAgcg             >  5:198463/1‑148 (MQ=255)
                                                                                                                            ggAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACTTCATAAAAACAAGcgcg           >  4:344090/1‑147 (MQ=255)
                                                                                                                            ggAACAGCCGGTCACCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACTTCATAAAAACAAGCGCGt          <  8:266243/148‑1 (MQ=255)
                                                                                                                             gAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACTTCATAAAAACAAGCGCGt          <  8:172447/147‑1 (MQ=255)
                                                                                                                             gAACAGCCGGTCACCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACTTCATAAAAACAAGCGCGt          <  8:259881/147‑1 (MQ=255)
                                                                                                                              aaCAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACTTCATAAAAACAAGCGCGTaa        <  2:357032/148‑1 (MQ=255)
                                                                                                                              aaCAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACTTCATAAAAACAAGCGCGTaa        <  6:198463/148‑1 (MQ=255)
                                                                                                                                    cGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACTTCATAAAAACAAGCGCGTAAACGTCg  <  3:250346/148‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                             
CAATGGCTTGTTCTTGCAGTGGTGTGGACGCCGTTAATTCACCCACGCCCTGCTCGCTCACTTTGGTGGACGACCCAAACCAGTTCCACGGGTTAGCGGCAGACCAGTTAACTGACGACATCGTGGAACAGCCGGTCAGCATCAATGGCATAGCGCATAACATTAAACGCAGCGATTTCATGTCACTTCCTTTGGTTATTCAATAACGTTGCTTGGAGTGCAAATTCACCAAAAAGTGCCGTTATTCTTTTACTTCATAAAAACAAGCGCGTAAACGTCG  >  NC_000913/2551036‑2551315

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