Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,971,159 G→A 12.2% intergenic (+26/‑112) ygdQ → / → ygdR UPF0053 inner membrane protein YgdQ/DUF903 domain‑containing lipoprotein YgdR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,971,1590GA12.2% 97.8 / 6.8 49intergenic (+26/‑112)ygdQ/ygdRUPF0053 inner membrane protein YgdQ/DUF903 domain‑containing lipoprotein YgdR
Reads supporting (aligned to +/- strand):  ref base G (20/23);  new base A (4/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

CTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCT  >  NC_000913/2971023‑2971297
                                                                                                                                        |                                                                                                                                          
cgggAAAGGTGCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCCCTCCCGCAGGAGAGCGTac                                                                                                                                      <  1:119312/141‑1 (MQ=255)
    atggTTTCGAAATCCGCGGACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACCTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactcac                                                                                                                                  <  2:37912/140‑1 (MQ=255)
    aaaGTTTCGACATCCACGTACCGAAAGCTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTactcac                                                                                                                                  <  1:509950/143‑1 (MQ=255)
          tCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCttt                                                                                                                            <  1:62609/143‑1 (MQ=255)
                   aCGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTcac                                                                                                                   >  1:377373/1‑143 (MQ=255)
                    cGTACCGAAAGGTTACATTTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACTCTcaca                                                                                                                  <  1:435229/143‑1 (MQ=255)
                           aagggTTAGATCTACTGCGCGATGTTCTTCTCTACTGCGGTTGACAGCCTCAACTTGATTCGCAACAAAAAGAAGCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAg                                                                                                           <  2:503244/140‑1 (MQ=255)
                           aaaGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAg                                                                                                           >  1:161227/1‑143 (MQ=255)
                           aaaGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAg                                                                                                           >  2:379031/1‑143 (MQ=255)
                            aaGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGa                                                                                                          <  1:588357/143‑1 (MQ=255)
                               gTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAttt                                                                                                       >  1:17102/1‑143 (MQ=255)
                               gTTACATCTACTTCCTGATTTTCTTCGCTATTGCGGTGGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGCTACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCGCTCCTTTTACCCTCACAGTTTAAGAttt                                                                                                       <  2:372989/143‑1 (MQ=255)
                                ttACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGAtttt                                                                                                      >  1:187249/1‑143 (MQ=255)
                                    aTCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTg                                                                                                  <  1:444970/143‑1 (MQ=255)
                                     tCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATAGTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTGAAGATCTTCCTGc                                                                                                 >  2:649220/1‑143 (MQ=255)
                                           tCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCaa                                                                                           <  1:160393/143‑1 (MQ=255)
                                               gATGTTATTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCAAAACATTTCTATTTATTGTTTAAGATTTTCCTGCTTTCAAGata                                                                                       <  1:316991/143‑1 (MQ=255)
                                                    tcttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATaa                                                                                  <  2:365133/143‑1 (MQ=255)
                                                      ttctCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTCCTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACg                                                                                >  1:241384/1‑143 (MQ=255)
                                                      ttctCTATTCCTGTTGCCGGCCTTAACTTGATTCGCAACTACATGCAGCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGTCGTACTCCCTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACg                                                                                <  2:504775/143‑1 (MQ=255)
                                                        ctctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGTTTTCAAGATATATAACGTc                                                                              >  1:455888/1‑143 (MQ=255)
                                                        ctctATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATCTTCCCGCTTTGAAGATATATACCGTc                                                                              >  1:278162/1‑143 (MQ=255)
                                                           tATTGCAGTGGAAAGCCTCAACGTGATTCGCATCGAAAATAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGt                                                                           <  2:138372/143‑1 (MQ=255)
                                                               gCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTCCCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATCTATAACGTCGGTTtat                                                                       >  2:103961/1‑143 (MQ=255)
                                                                  gTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGCTTTTCCTGCTTTCAAGATATATAACGTCGGTATATaaa                                                                    >  2:148589/1‑143 (MQ=255)
                                                                      gaaGCCTCAACTTTGTTCGGAACAAAAAGAATCCGCTCGGCTACTCCGTACGCTCTCCTGCGGGAGGGCTTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGACATATAACGTCGGTTTATAAACAGa                                                                <  1:271483/142‑1 (MQ=255)
                                                                       aaGCCTCACCTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGAc                                                               <  1:396337/143‑1 (MQ=255)
                                                                        aGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATAATCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACt                                                              <  2:168659/143‑1 (MQ=255)
                                                                         gCCTCAACTTGATTCGCAACAAAAACACTCCGCGCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACta                                                             <  1:341522/143‑1 (MQ=255)
                                                                                 ttGATTCTCACCATAACGAATCCGCTCTGATCCGCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCTCAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCAta                                                     <  2:184849/143‑1 (MQ=255)
                                                                                        gtcacaaAGAGAATCCGCTCTGATTCTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATACCGTCGGTTTATAAACAGACTATTATCATAGGTAGGc                                              <  2:605080/140‑1 (MQ=255)
                                                                                         caacaaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTCCTCACTCCTTTTACCCTCACCGTTTAAGATTTTCCTGCTTTCACGATATATAACGTCGGTTTATAAACAGACTATTATCATACGTAGGcc                                             >  2:425310/1‑143 (MQ=255)
                                                                                            caaAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAAGTAGGCCaaa                                          >  2:187857/1‑143 (MQ=255)
                                                                                             aaaaaGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCaagt                                         >  1:51841/1‑141 (MQ=255)
                                                                                                     tCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATACCGTCCGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATaa                                 >  1:222047/1‑143 (MQ=255)
                                                                                                           ctGATACTCCGTACGCTCTCCTGCCGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGAt                           <  2:439741/143‑1 (MQ=255)
                                                                                                           ctGACACTCCGTACGCTCTCCCGCGGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATACACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGAt                           <  1:137521/143‑1 (MQ=255)
                                                                                                            tGATACTCCGTACGCTCTCCTGCGGGAGGGCGTTCACACGCCTGTTGCCTTTACCGTTTAAGATGTTCCTGCTTTCAAGGTAGATAACGTCGGTTTATAAACAGACTATTATCAGAGCTCGGCCAGATGAGGATAACCGGacg                          >  2:449206/1‑141 (MQ=255)
                                                                                                             gttACTCCGTGCGCTGTGCTGCGGGAGGTCGTAGTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGAGATATAACGTCTGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGa                         <  2:5378/141‑1 (MQ=255)
                                                                                                             gATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGa                         <  1:379031/143‑1 (MQ=255)
                                                                                                              aTACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATCACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGAGGaa                        >  1:9888/1‑143 (MQ=255)
                                                                                                                 cTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGaaaaa                     >  2:11985/1‑143 (MQ=255)
                                                                                                                    agTACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATg                  >  2:295016/2‑143 (MQ=255)
                                                                                                                     gTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTGTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGAATAACAGGAGGAAGAAATgg                 >  2:19444/1‑143 (MQ=255)
                                                                                                                      tACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAGGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATggg                >  1:502304/1‑143 (MQ=255)
                                                                                                                      tACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATGTAACGTCGGTTTATAAACAGACTGTTATCATAGGTAGCACAAGTGAGGCTAACCGGAAGAACAAATggg                >  1:243929/1‑143 (MQ=255)
                                                                                                                      tACGCCCTCCCGCAGGAGAGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGACAAATGAGGATAACCGGATGAAAAAATggg                >  1:32766/1‑143 (MQ=255)
                                                                                                                         gCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGGTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAg             >  1:83924/1‑143 (MQ=255)
                                                                                                                                    ggAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCt  <  2:187249/143‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                          
CTGGAAAGTTTCGACATCCACGTACCGAAAGGTTACATCTACTTCGCGATGTTCTTCTCTATTGCGGTTGAAAGCCTCAACTTGATTCGCAACAAAAAGAATCCGCTCTGATACTCCGTACGCTCTCCTGCGGGAGGGCGTACTCACTCCTTTTACCCTCACAGTTTAAGATTTTCCTGCTTTCAAGATATATAACGTCGGTTTATAAACAGACTATTATCATAGGTAGGCCAAATGAGGATAACCGGATGAAAAAATGGGCAGTAATAATTTCT  >  NC_000913/2971023‑2971297

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