Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,381,918 A→C 6.5% intergenic (+16/‑29) ycjV → / → ompG pseudogene; putative ATP‑binding component of a transport system/outer membrane porin G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,381,9180AC 6.5% 219.8 / 7.3 62intergenic (+16/‑29)ycjV/ompGpseudogene; putative ATP‑binding component of a transport system/outer membrane porin G
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base C (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCG  >  NC_000913/1381797‑1382043
                                                                                                                         |                                                                                                                             
tGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCa                                                                                                                               <  3:647880/122‑1 (MQ=255)
tGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCa                                                                                                                               >  4:647880/1‑122 (MQ=255)
           tAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGTTGATGAAACTGAAATAGCATTTAGCTAAATATAGGTGGAAGGCATTCCCCCAGTTTAATACAAGGAACAATAATGAAAAAg                                                                                          >  4:112332/1‑148 (MQ=255)
                 aTTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAgg                                                                                                                     <  1:1061117/115‑1 (MQ=255)
                 aTTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAgg                                                                                                                     >  2:1061117/1‑115 (MQ=255)
                    aTCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGCAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGa                                                                                                         <  1:326337/124‑1 (MQ=255)
                    aTCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGa                                                                                                         >  2:326337/1‑124 (MQ=255)
                     tCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGa                                                                                                         >  4:1010775/1‑123 (MQ=255)
                     tCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGa                                                                                                         <  3:1010775/123‑1 (MQ=255)
                      cATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTccc                                                                                                                          <  4:818172/105‑1 (MQ=255)
                      cATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTccc                                                                                                                          >  3:818172/1‑105 (MQ=255)
                          cAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACaataa                                                                                                  <  1:505951/125‑1 (MQ=255)
                          cAGGAGAAAATATCACTATTCATTTTGATATGAAGAAATGTGATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACaataa                                                                                                  >  2:505951/1‑125 (MQ=255)
                             gagaAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGa                                                                                               <  3:637567/125‑1 (MQ=255)
                             gagaAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGa                                                                                               >  4:637567/1‑125 (MQ=255)
                                    tatCACTTTTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGa                                                                 <  3:112332/148‑1 (MQ=255)
                                         cTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTc                                                                                                                            <  2:184423/84‑1 (MQ=255)
                                         cTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTc                                                                                                                            >  1:184423/1‑84 (MQ=255)
                                           aTTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAATGCCTTccccctg                                                                                                                      >  4:693985/1‑86 (MQ=255)
                                           aTTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAATGCCTTccccctg                                                                                                                      <  3:693985/88‑3 (MQ=255)
                                             tCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTccccc                                                                                                                        <  1:902390/84‑1 (MQ=255)
                                             tCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTccccc                                                                                                                        >  2:902390/1‑84 (MQ=255)
                                                 tttGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACaata                                                                                                   <  3:1072377/101‑1 (MQ=255)
                                                 tttGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACaata                                                                                                   >  4:1072377/1‑101 (MQ=255)
                                                       atGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACaataa                                                                                                  <  1:869601/96‑1 (MQ=255)
                                                       atGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACaataa                                                                                                  >  2:869601/1‑96 (MQ=255)
                                                        tGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTAcc                                                                                  <  1:142316/111‑1 (MQ=255)
                                                        tGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTAcc                                                                                  >  2:142316/1‑111 (MQ=255)
                                                          aCGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGataat                                                                                                                >  2:574466/1‑79 (MQ=255)
                                                          aCGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGataat                                                                                                                <  1:574466/79‑1 (MQ=255)
                                                           cGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAATGCCTTccccc                                                                                                                        >  1:415945/1‑70 (MQ=255)
                                                           cGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAATGCCTTccccc                                                                                                                        <  2:415945/70‑1 (MQ=255)
                                                           cGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACt                                                                      <  2:965/120‑1 (MQ=255)
                                                           cGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACt                                                                      >  1:965/1‑120 (MQ=255)
                                                                   cATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACaataa                                                                                                  >  3:805422/1‑84 (MQ=255)
                                                                   cATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACaataa                                                                                                  <  4:805422/84‑1 (MQ=255)
                                                                     tttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGa                                                      >  4:197506/1‑126 (MQ=255)
                                                                     tttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGa                                                      <  3:197506/126‑1 (MQ=255)
                                                                      ttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACa                                                                                                      >  3:1107379/1‑77 (MQ=255)
                                                                      ttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACa                                                                                                      >  4:1014991/1‑77 (MQ=255)
                                                                      ttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACa                                                                                                      <  4:1107379/77‑1 (MQ=255)
                                                                      ttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACa                                                                                                      <  3:1014991/77‑1 (MQ=255)
                                                                      ttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGtt                                                                                        <  2:346869/91‑1 (MQ=255)
                                                                      ttCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGtt                                                                                        >  1:346869/1‑91 (MQ=255)
                                                                          ttGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGaaaa                                                                                            <  1:830760/83‑1 (MQ=255)
                                                                          ttGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGaaaa                                                                                            >  2:830760/1‑83 (MQ=255)
                                                                                      cGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGataat                                                                                                                >  1:411520/1‑51 (MQ=255)
                                                                                      cGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGataat                                                                                                                <  2:411520/51‑1 (MQ=255)
                                                                                        gAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGAtg                                                               >  2:257035/1‑98 (MQ=255)
                                                                                        gAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGAtg                                                               <  1:257035/98‑1 (MQ=255)
                                                                                                  ttCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTgg                                                                    >  2:915323/1‑83 (MQ=255)
                                                                                                  ttCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTgg                                                                    <  1:915323/83‑1 (MQ=255)
                                                                                                   tCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaa                              >  3:848177/1‑120 (MQ=255)
                                                                                                   tCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaa                              <  4:848177/120‑1 (MQ=255)
                                                                                                       tAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTAc                                                                                   <  2:182675/63‑1 (MQ=255)
                                                                                                       tAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTAc                                                                                   >  1:182675/1‑63 (MQ=255)
                                                                                                               gggggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaa                              <  2:9679/108‑1 (MQ=255)
                                                                                                               gggggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGaggaa                              >  1:9679/1‑108 (MQ=255)
                                                                                                                  ggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTgc                                             <  3:906117/90‑1 (MQ=255)
                                                                                                                  ggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTgc                                             >  4:906117/1‑90 (MQ=255)
                                                                                                                  ggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGgcg  >  3:664607/1‑133 (MQ=255)
                                                                                                                  ggAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCAAAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGgcg  <  4:664607/133‑1 (MQ=255)
                                                                                                                    aaGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAgg                                      >  2:916652/1‑95 (MQ=255)
                                                                                                                    aaGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGTGCGCAgg                                      <  1:916652/95‑1 (MQ=255)
                                                                                                                         |                                                                                                                             
TGCTGGTGCGTTAAATGATTATCATGCAGGAGAAAATATCACTATTCATTTTGATATGACGAAATGTCATTTCTTTGATGCAGAAACGGAAATAGCAATTCGCTAAATACAGGGGGAAGGCATTCCCCCAGGATAATACAAGGAACAATAATGAAAAAGTTATTACCCTGTACCGCACTGGTGATGTGTGCGGGAATGGCCTGCGCACAGGCCGAGGAAAGGAACGACTGGCACTTTAATATCGGCG  >  NC_000913/1381797‑1382043

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: