Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,632,292 C→A 15.6% intergenic (‑83/+599) fnbB ← / ← fnbA fibronectin‑binding protein B/fibronectin‑binding protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,632,2920CA15.6% 121.0 / 5.2 54intergenic (‑83/+599)fnbB/fnbAfibronectin‑binding protein B/fibronectin‑binding protein A
Reads supporting (aligned to +/- strand):  ref base C (23/22);  new base A (2/7);  total (25/29)
Fisher's exact test for biased strand distribution p-value = 1.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACTGAGGCCGCTCCCAATTTGTGTTTTCTTATGCCGTATCTAAGATTGCTTTTCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGAACACCTTGATATAACAGGCTTGTATTT  >  CP000730/2632152‑2632426
                                                                                                                                            |                                                                                                                                      
aTACTGAGGCCGCTCCCAATTTGTGTTTTCTTATGCCGTATCTAAGATTGCTTTTCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTGCGGTTTTTAAACGc                                                                                                                                        >  1:227936/1‑141 (MQ=255)
   cTGAGGCCGCTCCCAATTTGTGTTTTCTTATGCCGTATCTAAGATTGCTTTTCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAAt                                                                                                                                     >  2:313842/1‑141 (MQ=255)
      gtgCCGCTCCCAATTTGCTTTTTCTTATGCCGTATCTAAGATTGCTTTTCTCTATAATATTCTACCTTAAATGCAAAATTTATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTCTTAAAAGAAATGaa                                                                                                                                  <  2:519207/139‑1 (MQ=255)
      atccctcttccAATTTGTGTTTTCTTATTCCGTATCTAAGATTGCTTTTCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGaa                                                                                                                                  <  1:627479/132‑1 (MQ=255)
       ggCCGCTGCGAGTTTGTGTTGTCTTATGCCGTATCTAAGATTGCTTTTCACGATAATATTCTCCCTTAAATGCAACATTCATGTATTTTTAAAACTCAATATATGCAATTCTATGTCGTTCGGTTTTTAAAAGCAATGaaa                                                                                                                                 <  2:720273/141‑1 (MQ=255)
           gCTCCCAATTTGTGTTTTCTTATGCCGTATCTAAGATTGCTTTTCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGaaaaaaa                                                                                                                             >  1:404293/1‑141 (MQ=255)
             cctcAATTTGTGTTTTCTTATGCCGTATCTAAGATTGCTTTTCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGc                                                                                                                           <  2:37960/138‑1 (MQ=255)
                 aaTTTGTGTTTTCTTATGCCGAATCTAAGAATGCTTTTCACTATAATATTCTCCCTTAAATGCAAACTTCATTTATTTTTAAAACTCAATAAATTCAATTCTATATTGTTCGGTTTTTAAAAAAAAAAAAAAAAAGCGAGt                                                                                                                       <  2:569857/141‑1 (MQ=18)
                             cTTATGCCGTATCTAAGATTGCTTTTCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGtt                                                                                                           <  1:500701/141‑1 (MQ=255)
                             attaagcaatATCTAAGATTTATTTTCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAAAAAAAAAAAAAAGCGAGTTAATAAAAAGtt                                                                                                           <  2:34329/132‑1 (MQ=17)
                             aataaaccaaaTCAAAGATTGCTTTTCATTATAATATTCTTCATTCAATGCAAAATTCTTTTATTTTTAAGAACATTTATGTGGATATCAACATCGTTCAGTTTTTAAAAAAAAAAAAAAAAAGCGAGTTAATAAAAAGtt                                                                                                           <  2:345095/131‑1 (MQ=255)
                                        tattaacttGCTTTTCACACTAATCTCCACCCTTCAATGTAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATtgttgt                                                                                                <  1:245556/134‑1 (MQ=18)
                                           aaGATTGCTTTTCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTaa                                                                                             >  2:532032/1‑141 (MQ=255)
                                            aGATTGCTTTTCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAACACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAAc                                                                                            >  1:252707/1‑141 (MQ=255)
                                              aTTGCTTTTCTATATAATATGCTTCCTTAAGTGCAAAATTTATTTATTTTTTAAACAATATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAACAAGTGCGTTAAGAAAAAGTTAAGATTGTTGTTAACtt                                                                                          <  2:501915/141‑1 (MQ=17)
                                                   ttttCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATAAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGt                                                                                     >  2:90732/1‑141 (MQ=255)
                                                       cACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTGTGTATaa                                                                                 >  1:614765/1‑141 (MQ=255)
                                                       aaatataATATTCTCCCTTAAAGGCAAAATTCATTTATTTTTAAAACTCAATAATTGCAATTCTATATTGTTCGGTTTTTAAAAAAAAAAAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATaa                                                                                 <  2:165498/138‑1 (MQ=21)
                                                        aCTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAAt                                                                                <  1:69779/141‑1 (MQ=255)
                                                          tataATATTCTCCCTTACATGAAAAATTCATTTATTTTTAAAACTCTATAAGTGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAATAGTTAAGATTGTTGTTAACTTTATGTATAATGa                                                                              <  2:219871/141‑1 (MQ=21)
                                                            taataTTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGACATCACAAAATGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGt                                                                            >  1:21598/1‑141 (MQ=25)
                                                            aactacttccccTAAGATGCAAACTTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGt                                                                            <  1:481423/132‑1 (MQ=21)
                                                                atTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGttttt                                                                        >  2:420475/1‑141 (MQ=255)
                                                                   ctctctTAGATGCAAATTTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATCGTTCTGTTTTGAAAAAAAATAAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTttat                                                                     <  2:269724/141‑1 (MQ=255)
                                                                   ctcCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTttat                                                                     >  1:574591/1‑141 (MQ=255)
                                                                   ctcCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTttat                                                                     >  2:131853/1‑141 (MQ=255)
                                                                   ctcCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGGTTTttat                                                                     >  1:665607/1‑141 (MQ=255)
                                                                   ctcCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGAAATGAAAAAAAGCGAGTAAATAAAAAGTTATGATTGTTGTTAACTTTATGTAGAAGGAGTTGTttat                                                                     >  2:400133/1‑141 (MQ=21)
                                                                       cTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTattatt                                                                 <  1:368951/141‑1 (MQ=255)
                                                                          aaaTGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGa                                                              <  1:347495/141‑1 (MQ=255)
                                                                                 aaaTTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTcac                                                       >  1:657190/1‑141 (MQ=255)
                                                                                      cATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACatata                                                  <  1:689371/141‑1 (MQ=255)
                                                                                            tttttAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGc                                            <  1:119461/141‑1 (MQ=255)
                                                                                                tAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCACTAAAACAAACCGATGTCATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGATTTTTTTATTATTTGAAACTCACATATATATTGCATAc                                        >  1:177723/1‑141 (MQ=25)
                                                                                                 aaaaCTCAATAAATGTAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAAATCACATATATATTGCATACa                                       >  1:562775/1‑141 (MQ=255)
                                                                                                 aaaaCTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACa                                       <  1:237771/141‑1 (MQ=255)
                                                                                                 aaaaCTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACAGATATATTGCATACa                                       >  1:643457/1‑141 (MQ=255)
                                                                                                         ataaccagaaTTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCtt                               <  1:744029/133‑1 (MQ=255)
                                                                                                          tAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTg                              >  1:454584/1‑141 (MQ=255)
                                                                                                           taaTGCAATTTTATATTGTTCGATTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGa                             <  2:669453/140‑1 (MQ=255)
                                                                                                           aagtGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGa                             <  2:744401/138‑1 (MQ=255)
                                                                                                            aatgagaTTCTATCTTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAATATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGaa                            <  2:121809/135‑1 (MQ=255)
                                                                                                            aaTGCAATTCTATATTGGTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGGTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTAAAATGAGTTTTTTATTATTTGAAACTCAGATATATATTGCATACAAAGCTCTTGaa                            >  1:49634/1‑141 (MQ=255)
                                                                                                              tGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTCTGTATAATGAGTTTTTTATTATTTGCAACTAACATATATTTTGCATACAAAGCTCTTGAaca                          >  1:94176/1‑141 (MQ=255)
                                                                                                              tGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAGAATCAACTATATATGGCTTACAAAGCTCTTGAaca                          >  2:726989/1‑141 (MQ=25)
                                                                                                                 aaTTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGAACACCt                       >  2:583804/1‑141 (MQ=255)
                                                                                                                   ttCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTATTTTATTATTTGAAACTCACATATATATTGCATACAAAGCACTTTAACACCtta                     >  2:496754/1‑140 (MQ=255)
                                                                                                                    tCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGAACACCTTGa                    <  1:523494/141‑1 (MQ=255)
                                                                                                                             ttCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTCAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGAACACTTTGATATAACAgg           >  2:196250/1‑141 (MQ=255)
                                                                                                                             ttCGGTTTTTAAAAGCAAAGAAGAAAAGCGAGTTAATAAAAAGTTAAGGTTGTTGTTAACTTTATGTATAATGAGCTTTTTATTATTTGAAAGGTACATATATCTTGCATACAAAGCTCTTGAACCCCTTGCTATAACAgg           >  2:276794/1‑141 (MQ=255)
                                                                                                                              tCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGAACACCTTGATATAACAGGc          <  1:420475/141‑1 (MQ=255)
                                                                                                                               cGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGAACACCTTGATATAACAGGCt         <  1:131853/141‑1 (MQ=255)
                                                                                                                                ggTTTTAAAAAAAAAAAAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGAACACCTTGATATAACAGGCtt        <  2:150437/141‑1 (MQ=255)
                                                                                                                                   ttttAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGAACACCTTGATATAACAGGCTTGTa     <  1:532032/141‑1 (MQ=255)
                                                                                                                                      tAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTTTTGTTGACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGAACACCTTGATATAACAGGCTTGTAttt  >  1:332890/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                      
ATACTGAGGCCGCTCCCAATTTGTGTTTTCTTATGCCGTATCTAAGATTGCTTTTCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGAACACCTTGATATAACAGGCTTGTATTT  >  CP000730/2632152‑2632426

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