Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,632,296 G→A 14.2% intergenic (‑87/+595) 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,2960GA14.2% 126.4 / 3.2 53intergenic (‑87/+595)fnbB/fnbAfibronectin‑binding protein B/fibronectin‑binding protein A
Reads supporting (aligned to +/- strand):  ref base G (21/23);  new base A (2/6);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCCGCTCCCAATTTGTGTTTTCTTATGCCGTATCTAAGATTGCTTTTCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGAACACCTTGATATAACAGGCTTGTATTTTTTATACTTA  >  CP000730/2632158‑2632436
                                                                                                                                          |                                                                                                                                            
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                                                                               >  2:131853/1‑141 (MQ=255)
                                                             ctcCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTttat                                                                               >  1:574591/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)
                                                                                                                                          gAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGAACACCTTGATATAACAGGCTTGTATTTTTTAtactta  >  2:149507/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                            
AGGCCGCTCCCAATTTGTGTTTTCTTATGCCGTATCTAAGATTGCTTTTCACTATAATATTCTCCCTTAAATGCAAAATTCATTTATTTTTAAAACTCAATAAATGCAATTCTATATTGTTCGGTTTTTAAAAGCAATGAAAAAAAGCGAGTTAATAAAAAGTTAAGATTGTTGTTAACTTTATGTATAATGAGTTTTTTATTATTTGAAACTCACATATATATTGCATACAAAGCTCTTGAACACCTTGATATAACAGGCTTGTATTTTTTATACTTA  >  CP000730/2632158‑2632436

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: ‑