Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 661,369 T→C N4D (AAT→GAT)  ybeF ← LysR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913661,3690TC98.3% 210.2 / NA 60N4D (AAT→GAT) ybeFLysR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (36/23);  minor base A (1/0);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGC  >  NC_000913/661275‑661463
                                                                                              |                                                                                              
gTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTAt                                                                                           >  1:1364754/1‑100 (MQ=255)
    gAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCAc                                                                                       >  1:794195/1‑100 (MQ=255)
      aGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGc                                                                                     <  2:1493492/100‑1 (MQ=255)
       ggTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCa                                                                                    <  1:1161471/100‑1 (MQ=255)
         ttAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACt                                                                                  <  2:1005060/100‑1 (MQ=255)
            aGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTAc                                                                               >  1:799824/1‑100 (MQ=255)
             gATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACt                                                                              >  2:990897/1‑100 (MQ=255)
             gATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACt                                                                              <  1:57583/100‑1 (MQ=255)
             gATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACt                                                                              >  2:455844/1‑100 (MQ=255)
              aTCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTc                                                                             <  2:504784/100‑1 (MQ=255)
               tCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTcc                                                                            >  1:1263347/1‑100 (MQ=255)
                cAATATTTCGTAAGGTTGTAAATATTTGTGCTTTACCTTCTGACGACTTTCTGCTTAAACAGGGGTCAATTTGATTATCACTATCCACGCACTTAc                                                                               <  2:1381654/96‑1 (MQ=255)
                   tatTTCGTAAGGTTGTAAATAATTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATAACTATCCACGCACTTACTCCAAtt                                                                        >  2:478533/1‑100 (MQ=255)
                    atTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAAttt                                                                       <  1:1541165/100‑1 (MQ=255)
                     tttCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAAtttt                                                                      <  2:979508/100‑1 (MQ=255)
                       tCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTAt                                                                    >  1:706210/1‑100 (MQ=255)
                         gTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTc                                                                  <  1:1140905/100‑1 (MQ=255)
                             ggTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATgg                                                              <  2:651338/100‑1 (MQ=255)
                             ggTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATg                                                               >  2:845210/1‑99 (MQ=255)
                                 gTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGaaaa                                                          >  2:7257/1‑100 (MQ=255)
                                     atatTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAaataa                                                      <  1:777746/100‑1 (MQ=255)
                                      tatTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAaataat                                                     <  2:1168187/100‑1 (MQ=255)
                                       atTTGTGGTTTTCATTCTGACTACTTTCTGCTTAAACAGGTTTCAATTTTATTGTCAAGATCCACGAAATTACTCCAATTTTATGCATGGAAAAataata                                                    <  2:1521406/100‑1 (MQ=255)
                                         ttGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATAtt                                                  >  1:981268/1‑100 (MQ=255)
                                         ttGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATAtt                                                  >  2:320712/1‑100 (MQ=255)
                                          tgtgGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATAttt                                                 >  1:846876/1‑100 (MQ=255)
                                           gtgGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTa                                                <  2:44402/100‑1 (MQ=255)
                                              gTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTaaaa                                             >  2:318632/1‑100 (MQ=255)
                                                 tACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAAt                                          >  2:299526/1‑100 (MQ=255)
                                                 tACCTTCTGACGACTTTCTGCTTAAACAGGGTGCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAAt                                          <  1:1514576/100‑1 (MQ=255)
                                                   ccTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAAtaa                                        <  1:1071773/100‑3 (MQ=255)
                                                   ccTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTa                                        >  1:1018871/1‑100 (MQ=255)
                                                               tttCTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTAt                            >  1:310487/1‑100 (MQ=255)
                                                                ttCTGCTTAAACATGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATg                           >  2:689147/1‑100 (MQ=255)
                                                                ttCTGCTTAAACAGAGTTCAATTTGAATATCAATATACAAGCAATTACTAAAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTa                             <  2:170773/98‑1 (MQ=255)
                                                                  cTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTc                         >  2:1070368/1‑100 (MQ=255)
                                                                  cTGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTc                         >  1:412298/1‑100 (MQ=255)
                                                                   tGCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCt                        >  2:767829/1‑100 (MQ=255)
                                                                    gCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTa                       >  2:1364620/1‑100 (MQ=255)
                                                                    gCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTa                       >  2:1449653/1‑100 (MQ=255)
                                                                    gCTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTa                       >  2:250196/1‑100 (MQ=255)
                                                                     cTTAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTaa                      >  2:436797/1‑100 (MQ=255)
                                                                      ttAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAAt                     <  2:940649/100‑1 (MQ=255)
                                                                      ttAAACAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAAt                     >  1:728080/1‑100 (MQ=255)
                                                                        aaaCAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTg                   >  2:593999/1‑100 (MQ=255)
                                                                        aaaCAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGCCTAATTg                   >  2:137098/1‑100 (MQ=255)
                                                                         aaCAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTgg                  >  1:1100378/1‑100 (MQ=255)
                                                                         aaCAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTgg                  >  1:97157/1‑100 (MQ=255)
                                                                           catggTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTgga                 <  2:492704/96‑1 (MQ=255)
                                                                           cAGGGTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTggaa                >  2:329084/1‑100 (MQ=255)
                                                                             gggTTCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTggaacg              >  1:224680/1‑100 (MQ=255)
                                                                                ttCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTggaacggaa           <  1:625688/100‑1 (MQ=255)
                                                                                 tCAATTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTggaacggaac          <  1:254056/100‑1 (MQ=255)
                                                                                 tCAATTTGATGATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAACTTCCAATAATCTTATGTCTAATTggaacggaac          <  2:1041696/100‑1 (MQ=255)
                                                                                   aaTTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGc        >  1:1557019/1‑100 (MQ=255)
                                                                                   aaTTTGATTATCACTATCCACACACTTACTCCAATTTTATTCATTGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTggaacgg             <  2:1540267/95‑1 (MQ=255)
                                                                                    aTTTGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCt       >  2:334802/1‑100 (MQ=255)
                                                                                     tttGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCtt      >  1:309508/1‑100 (MQ=255)
                                                                                      ttGATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCttt     >  1:1487145/1‑100 (MQ=255)
                                                                                        tattatcACTATCCACACACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTg   <  2:145395/99‑1 (MQ=255)
                                                                                        gATTATCACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTg   >  1:485109/1‑100 (MQ=255)
                                                                                         attatcACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGc  >  2:17590/1‑100 (MQ=255)
                                                                                         attatcACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTAAAATAATCTTATGTATAATTGGAACAGAACGCTTTTg   >  2:184454/1‑99 (MQ=255)
                                                                                              |                                                                                              
GTCAGAAGGTTAAGATCAATATTTCGTAAGGTTGTAAATATTTGTGGTTTACCTTCTGACGACTTTCTGCTTAAACAGGGTTCAATTTGATTATTACTATCCACGCACTTACTCCAATTTTATTCATGGAAAAATAATATTTAAAAAATTACAATAATCTTATGTCTAATTGGAACGGAACGCTTTTGC  >  NC_000913/661275‑661463

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