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.C97.6% 139.8 / ‑2.9 42intergenic (‑304/‑177)yieP/BW25113_RS19520FadR/GntR family transcriptional regulator/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (18/23);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGAT  >  NZ_CP009273/3934909‑3935075
                                                                                   |                                                                                    
gTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACAcccccgc                                                                               <  2:532203/90‑4 (MQ=18)
              ccTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGaaaaa                                                                  >  2:957058/1‑90 (MQ=255)
              ccTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAaaag                                                                 >  2:535898/1‑90 (MQ=33)
                 gagaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagca                                                               >  2:748409/1‑90 (MQ=255)
                   gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa                                                             >  2:968924/1‑90 (MQ=255)
                   gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa                                                             >  2:940202/1‑90 (MQ=255)
                   gaACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaa                                                             >  2:207606/1‑90 (MQ=255)
                     aCTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaagc                                                           <  2:552209/90‑1 (MQ=255)
                     aCTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaagc                                                           <  1:174421/90‑1 (MQ=255)
                     aCTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAaaagcaaagc                                                           <  1:511101/90‑1 (MQ=255)
                      cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCg                                                          >  1:1041163/1‑90 (MQ=255)
                      cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCg                                                          >  1:361579/1‑90 (MQ=255)
                      cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACAcccccgc                                                                               >  1:606091/1‑65 (MQ=25)
                      cTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACAcccccgc                                                                               <  2:606091/68‑4 (MQ=25)
                        ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGc                                                        >  1:1002833/1‑90 (MQ=255)
                        ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGc                                                        >  2:925974/1‑90 (MQ=255)
                        ccGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGc                                                        >  2:544190/1‑90 (MQ=255)
                           gCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACt                                                     <  2:1026277/90‑1 (MQ=255)
                           gCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACt                                                     <  1:1209614/90‑1 (MQ=255)
                                     cAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAAc                                           <  1:637943/90‑1 (MQ=255)
                                                gCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGa                                <  2:1012560/90‑1 (MQ=255)
                                                 cTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGAc                               <  1:789475/90‑1 (MQ=255)
                                                 cTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGAc                               <  2:141198/90‑1 (MQ=255)
                                                 cTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGAc                               <  2:753318/90‑1 (MQ=255)
                                                    gACTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAAt                            <  1:1090441/90‑1 (MQ=255)
                                                     aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc                           <  2:901329/90‑1 (MQ=255)
                                                     aCTCTGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATc                           <  1:764795/90‑1 (MQ=255)
                                                      ctctGTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCt                          <  2:847084/90‑1 (MQ=255)
                                                          gTAGCGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCtgtgt                      >  2:717406/1‑90 (MQ=255)
                                                              cGGGAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGc                  <  1:1114497/90‑1 (MQ=255)
                                                               gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGAc                               <  1:85139/76‑1 (MQ=255)
                                                               gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGAc                               >  2:85139/1‑76 (MQ=255)
                                                               gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                 >  2:735197/1‑90 (MQ=255)
                                                               gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                 >  1:781988/1‑90 (MQ=255)
                                                               gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                 >  2:1090754/1‑90 (MQ=255)
                                                               gggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCa                 >  1:860336/1‑90 (MQ=255)
                                                                ggAAAGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCAc                <  2:831100/90‑1 (MQ=255)
                                                                    aGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGa            >  1:260771/1‑90 (MQ=255)
                                                                    aGCGTATTATGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGa            <  2:882513/90‑1 (MQ=255)
                                                                          ttatGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATAc      <  1:484449/90‑1 (MQ=255)
                                                                           tatGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACg     <  2:577447/90‑1 (MQ=255)
                                                                           tatGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACg     <  1:150423/90‑1 (MQ=255)
                                                                           tatGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACg     <  1:171634/90‑1 (MQ=255)
                                                                             tGCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGa   <  2:1026504/90‑1 (MQ=255)
                                                                              gCACACCCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGAt  >  1:771940/1‑90 (MQ=255)
                                                                                   |                                                                                    
GTCAGCGGGGTTCTCCTGAGAACTCCGGCAGAGAAAGCAAAAATAAATGCTTGACTCTGTAGCGGGAAAGCGTATTATGCACACCCCGCGCCGCTGAGAAAAAGCAAAGCGGCACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGAT  >  NZ_CP009273/3934909‑3935075

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