Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 614,181 T→C 19.6% D831D (GAT→GAC sdrC → Ser‑Asp rich fibrinogen/bone sialoprotein‑binding protein SdrC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730614,1810TC19.6% 65.1 / 2.9 39D831D (GAT→GACsdrCSer‑Asp rich fibrinogen/bone sialoprotein‑binding protein SdrC
Reads supporting (aligned to +/- strand):  ref base T (19/10);  new base C (2/6);  total (22/17)
Fisher's exact test for biased strand distribution p-value = 5.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAACGACTC  >  CP000730/614041‑614321
                                                                                                                                             |                                                                                                                                            
gactcagatccaGAGAGCGACTCAGACTCAGATACCGATTCGGACTCAGACACTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTgattcagat                                                                                                                                              <  2:267121/141‑1 (MQ=21)
   tcagactcagaTAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTgattcagattca                                                                                                                                           <  1:33033/141‑1 (MQ=255)
     agacccagaTAGCGATGCAAATTCAGACAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGAttcagattcaga                                                                                                                                         <  1:20362/141‑1 (MQ=255)
      gactcagaTAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATCCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGAGTCAGATTCAGAc                                                                                                                                        >  1:207469/1‑141 (MQ=255)
               aGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACtca                                                                                                                               >  1:132762/1‑141 (MQ=255)
                gCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGGGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGGGATTCAGATTCCGACAGCGCCtcag                                                                                                                              >  1:67099/1‑141 (MQ=255)
                 cGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTAAGATTCAGACTCAGCTAGCGACTCGGATTCTGACGGCGATTCAGACTCAGACAGCGACTCAGAGTCAGACAGTGATTCAGAGTCCGACAGCGACtgaga                                                                                                                             >  1:108671/1‑141 (MQ=21)
                      cagattcagatCGCGACTCGGACTCAGACTGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACtcagattcag                                                                                                                        <  2:171246/141‑1 (MQ=255)
                         attcagatAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcagata                                                                                                                     >  1:55247/1‑141 (MQ=255)
                                      tccgcaCTCAGACAGTGATTCAGATTCACACACAGATAGCCACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcagatagtgactcagactc                                                                                                        <  2:259855/136‑1 (MQ=17)
                                         ggACTCAGACCGTGACTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcagatagcgactcagactcaga                                                                                                     <  1:292712/141‑1 (MQ=255)
                                          gacccAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACagcgactcagactcagatagcgactcagactcagat                                                                                                    >  2:25168/5‑141 (MQ=255)
                                           aCTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcagatagcgacccagactcagata                                                                                                   >  1:181085/1‑141 (MQ=255)
                                           aCTCAGACAGAGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGTTtcagatagcgactcagactcagaaa                                                                                                   >  1:244751/1‑139 (MQ=255)
                                                gACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGCCAGCGATTCAGACTCAGACAGCGACTCAGAATCAGACAGTGATTCAGATGCAGACAGCGACTCAGATtcagatagcgactcagactccgatagccgc                                                                                              >  1:25974/1‑138 (MQ=17)
                                                   aGTGATTCAGATTCAGCCTCAGATAGCGACTCAGATTCTGATAGCGAGTCAGACTCAGACAGCGACCCAGACCCAGACAGGGATTCAGATTCAGACAGCGACTCAGATtcagatagcgactcagactcagatagcgactca                                                                                           >  2:316402/1‑141 (MQ=17)
                                                                   aCTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcagatagcgactcagactcagatagcgacccagactcagatagcgacc                                                                           >  2:1163/1‑140 (MQ=255)
                                                                           agggACTCAGATTCGGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCATACTCAGATAGCGACTCAGACTCg                                                                   <  1:291566/138‑1 (MQ=255)
                                                                            gCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATtcagatagcgactcagactcagatggagaatcagactcaggtagcgattcagtacggg                                                                  >  2:192434/1‑135 (MQ=17)
                                                                                ggcAGATTCTGATAGCGATTCAGACTCTGATCGCGAGTCAGGCTCAGACAGTGATTCAGCTTCAGACAGCGACTCAGCTTCAGATAGCGATTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAg                                                              <  2:2652/139‑1 (MQ=14)
                                                                                      ttCTGACAGCGATTCAGACTCAGACAGCGACTCCGACTCAGACAGTGCTTAACGCTACTACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGTCTCAGATAGCGACTCAGAACCGGATAGCgattc                                                        >  1:63327/1‑141 (MQ=17)
                                                                                        cTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACCGTGCTTAGAGCCCATACAGCGCATATGCCGCAGATAGCGACTCAGACTCAGATAACGACCCAGACTCAGATAGCGACTCAGACTCGGATAGCgattcag                                                      >  1:329924/1‑141 (MQ=17)
                                                                                        cTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCgattcag                                                      >  2:214529/1‑141 (MQ=255)
                                                                                          gACAGCGATTCAGACTCAGATAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCgattcagat                                                    >  1:216616/1‑141 (MQ=255)
                                                                                               cGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCAAttgagattcaga                                               >  1:56925/1‑141 (MQ=255)
                                                                                               cGATTCAGACTCAGAAAGTGTTTCGGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGAttcagagtcaga                                               >  2:322901/1‑141 (MQ=14)
                                                                                                  ttttttttttgacAGCGACTCAGCCTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCCGACTCAGATAGCGACTCAGACTCGGAAAGCGATTCACATGCAGACAg                                            >  2:204209/11‑141 (MQ=17)
                                                                                                        ctcagacAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACtc                                      <  1:150122/141‑1 (MQ=255)
                                                                                                        atcagacTCCCACCCAGACCCAGACAGCGACTCAGACTCAGACAGCGACACAGATTCAAACACCGACTCAGACTCAGATAGCGACTCAGACTCAGACAGCGACTCAGACTCTGATAGCGATTCAGATTCAGACAGCGACtc                                      <  2:328674/140‑1 (MQ=11)
                                                                                                        accagacACCGAACCACACTCAGACAGCGACCCAGACTCAGACACCGACTCAGATTCAGATACCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACtc                                      <  2:328796/139‑1 (MQ=11)
                                                                                                          cagacAGCGACTCAGACTCAGACCGAGATTCAGACTCAGACAGCGACTCAGATTCAGATAGCACCTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACtcag                                    <  2:230444/141‑1 (MQ=14)
                                                                                                          cagacAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACtcag                                    >  2:331759/1‑141 (MQ=255)
                                                                                                          cagacAGCGACTCAGACTCAGACAGCGACTCCGACTCAGACAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGCCTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACtcag                                    <  2:230539/141‑1 (MQ=14)
                                                                                                              aaGCGACTCCGACTGAGACAGTGATTCCGATTCAGACAGCGACTCCGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGGTAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACtcagattc                                <  2:22817/140‑1 (MQ=17)
                                                                                                               cgccacccagactcagacAGGGACTCAGACTCAGACACCGACTCAGGTTCAGATAGCGACTCAGACTCAGACAGCGACTCAGACTCAGATAGCTACTCAGCCTCGGATAGCGATTCAGATTCAGACAGCGACtcagattca                               <  2:216616/134‑1 (MQ=14)
                                                                                                                      cagactcagacAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGCCCCAGATTCAGATAGca                        >  2:29855/1‑140 (MQ=255)
                                                                                                                       agactcagacAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGa                       <  1:70195/141‑1 (MQ=255)
                                                                                                                        cgcccagacAGGTATTCAGCTCCAGACAGCGACTCAGACTCAGAGAGCGACTCCGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGAt                      <  2:148453/137‑1 (MQ=14)
                                                                                                                        ccctccgacAGCGACTCCCACTCCGACAGCGACTCAGAATCAGACAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGAt                      <  2:39265/139‑1 (MQ=14)
                                                                                                                          ctcagatCGTGACTCTGAATCAGACGGCGACCCAGACTCAGACAGCGACTCAGACTCAGATAGCGACTCAGACTCCGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTc                    <  2:74611/141‑1 (MQ=11)
                                                                                                                                  gcgattcagattgaGACAGCGACCCCGCTTCAGATCGCGACTCCGACGCAGATAGCGACGCAGACTCAGATAGCGATTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAg            <  2:65/139‑1 (MQ=17)
                                                                                                                                  gTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATGGCGATTCGGACTCAg            >  2:49845/1‑141 (MQ=255)
                                                                                                                                   tGATTCAGATTCAGACAGCGACTCAGTTTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACACGGATAGCGATTCCGATTCAGACAGCGACTCAGATTCCGATAGCGATTCGGACTCAGa           <  1:202660/141‑1 (MQ=255)
                                                                                                                                           attcagaCAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAACGACt   >  1:62450/1‑141 (MQ=255)
                                                                                                                                             ttcagaCCGCGACTCATAATCTTATATCGACTCAGACTCAGACTGCTCCCCCGACCCAGTCATCCACTAAGATTCGGATAGCGATTCAGGTTCAGACCGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAGCGACtc  <  2:142415/141‑1 (MQ=2)
                                                                                                                                             |                                                                                                                                            
GACTCAGACTCAGATAGCGATTCAGATTCAGATAGCGATTCGGACTCAGACAGTGATTCAGATTCAGACTCAGATAGCGACTCAGATTCTGACAGCGATTCAGACTCAGACAGCGACTCAGACTCAGACAGTGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCAGATAGCGACTCAGACTCGGATAGCGATTCAGATTCAGACAGCGACTCAGATTCAGATAGCGATTCGGACTCAGACAACGACTC  >  CP000730/614041‑614321

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