Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,934,995 +C intergenic (‑308/‑173) yieP ← / → BW25113_RS19520 FadR/GntR family transcriptional regulator/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,934,9911.C93.2% 123.2 / 3.3 44intergenic (‑304/‑177)yieP/BW25113_RS19520FadR/GntR family transcriptional regulator/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (24/17);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 8.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

AGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTT  >  NZ_CP009273/3934912‑3935079
                                                                                  |                                                                                        
aGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGcgcc                                                                                 >  1:1438301/1‑90 (MQ=255)
      gTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTgag                                                                           >  1:1175514/1‑90 (MQ=255)
                  aCTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGTGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaagc                                                               >  1:263560/1‑90 (MQ=11)
                  aCTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaagc                                                               <  1:513220/90‑1 (MQ=255)
                  aCTCCGGCAGAGAAAGCAAAAATAAATGCNTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGACAAaaagcaaagc                                                               <  2:546192/90‑1 (MQ=14)
                   cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCg                                                              >  2:362498/1‑90 (MQ=255)
                   cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCg                                                              >  1:769/1‑90 (MQ=255)
                     ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGc                                                            >  2:1116623/1‑90 (MQ=255)
                       ggCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACt                                                         <  1:1519299/90‑1 (MQ=39)
                          agagAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGc                                                       <  1:1083780/90‑1 (MQ=255)
                             gaAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGctct                                                    >  1:742302/1‑90 (MQ=255)
                               aaGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCttt                                                  >  1:1117749/1‑90 (MQ=255)
                               aaGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCttt                                                  >  1:622408/1‑90 (MQ=255)
                                 gCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTaa                                                >  1:229772/1‑90 (MQ=255)
                                    aaaaTAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACaa                                             <  2:1405707/90‑1 (MQ=255)
                                    aaaaTAAATGCTAGACTCTGTAGCGGGAAATGGTATTATGCACACCCCCGTGCCGCTGAGAAAAAGCAAAGCGGCGCTGCTCTTTAACaa                                             <  2:804465/90‑1 (MQ=11)
                                          aatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATc                                       <  1:407886/90‑1 (MQ=255)
                                          aatGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATc                                       <  1:1314090/90‑1 (MQ=255)
                                            tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg                                     <  1:1229032/90‑1 (MQ=255)
                                            tGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg                                     <  2:224869/90‑1 (MQ=255)
                                               ttGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACa                                  >  2:537319/1‑90 (MQ=255)
                                                tGACTCTGTAGCGGGAAAGGGTATTATGCACCCCCCCGCGCCGCCGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACaa                                 <  2:1320846/90‑1 (MQ=14)
                                                tGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACaa                                 >  1:599780/1‑90 (MQ=255)
                                                     ctGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgt                            >  2:115736/1‑90 (MQ=255)
                                                           cGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGc                      <  2:1422293/90‑1 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg                                     >  1:138793/1‑74 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAg                                     <  2:138793/74‑1 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc                               <  2:1213696/80‑1 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc                               >  1:1213696/1‑80 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                     >  2:1419379/1‑90 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                     >  2:802260/1‑90 (MQ=255)
                                                            gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGCGTGGGCa                     >  1:839649/1‑90 (MQ=11)
                                                               aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                  >  1:931688/1‑90 (MQ=255)
                                                               aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                  >  1:1315625/1‑90 (MQ=255)
                                                               aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                  >  2:1558873/1‑90 (MQ=255)
                                                               aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                  >  1:654303/1‑90 (MQ=255)
                                                               aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                  >  2:357853/1‑90 (MQ=255)
                                                               aaaGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTc                  <  1:537319/90‑1 (MQ=255)
                                                                 aGCGTATTATGCACACCCCCGCGCCGCTGAGaaaaa                                                                      <  2:292856/36‑1 (MQ=11)
                                                                 aGCGTATTATGCACACCCCCGCGCCGCTGAGaaaaa                                                                      >  1:292856/1‑36 (MQ=11)
                                                                  gCGTATTATGCAACCCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGac               <  1:293550/90‑2 (MQ=11)
                                                                        tatGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACg         <  1:1419379/90‑1 (MQ=255)
                                                                          tGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGa       <  1:237035/90‑1 (MQ=255)
                                                                          tGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGa       <  2:1035197/90‑1 (MQ=255)
                                                                                cacCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCGTTAACAATTTATCAGACAATCTATGTGGGCACTCGAAGATACGGATTCtt  >  2:1421758/1‑90 (MQ=255)
                                                                                  |                                                                                        
AGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTT  >  NZ_CP009273/3934912‑3935079

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