Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 129,029 A→G 8.7% L199L (TTA→CTA)  spa ← immunoglobulin G binding protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730129,0290AG 8.7% 117.9 / 4.0 46L199L (TTA→CTA) spaimmunoglobulin G binding protein A
Reads supporting (aligned to +/- strand):  ref base A (28/14);  new base G (2/2);  total (30/16)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

AGGTAAATGTAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAATTTGTTATCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTCGTTTAATTTTTTAGCTTCACCT  >  CP000730/128895‑129168
                                                                                                                                      |                                                                                                                                           
aGGTAAATGTAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAATTTGTTATCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTg                                                                                                                                       >  1:176252/1‑141 (MQ=255)
      aTGTAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAATTTGTTATCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGaa                                                                                                                                 >  1:278414/1‑141 (MQ=255)
               ttCATAGAAAGCATTTTGTTGTTCTTTGTTGAATTTGTTATCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCg                                                                                                                        <  2:287168/141‑1 (MQ=255)
                tCATAGAAAGCATTTTGTTGTTCTTTGTTGAATTTGTTATCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCgt                                                                                                                       <  1:345924/141‑1 (MQ=255)
                 cATAGAAAGCATTTTGTTGTTCTTTGTTGAATTTGTTATCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCgtt                                                                                                                      >  1:250419/1‑141 (MQ=255)
                          cATTTTGTTGTTCTTTGGTGAATTTGTTGTCCGCTTTCGGTGCTTGAGATTCATTTAGCTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAGGCTTTGGATGAAACCATTGCGTTGTTCTTCGt                                                                                                             <  2:96631/141‑1 (MQ=255)
                             tttGTTGTTCTTTGTTGAATTTGTTGTCCGCTTTTGGTGCTTGAGATTCATTTAGCTTTTTAGCTTCTGATAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTa                                                                                                          <  2:163367/141‑1 (MQ=255)
                              ttgttgTTCTTTGTTGAATTTGTTGTCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTaa                                                                                                         <  2:113726/141‑1 (MQ=255)
                              ttgttgTTCTTTGTTGAATTTGTTATCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTaa                                                                                                         >  2:1057/1‑141 (MQ=255)
                              ttgttgTTCTTTGTTGAATTTGTTATCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTaa                                                                                                         >  2:1055/1‑141 (MQ=255)
                              ttgttgTTCTTTGTTGAATTTGTTATCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTaa                                                                                                         >  2:64786/1‑141 (MQ=255)
                                      cTTTGTTGAATTTGTTGTCCGCTTTCGGTGCTTGAGATTCATTTAGCTTTTTAGCTTCTGCCAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCa                                                                                                 <  2:106503/141‑1 (MQ=255)
                                        ttgttgAATTTGTTATCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCata                                                                                               <  1:142162/141‑1 (MQ=255)
                                             gAATTTGTTATCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCaa                                                                                          >  2:196699/1‑141 (MQ=255)
                                                tttGTTATCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGAt                                                                                       >  1:88824/1‑141 (MQ=255)
                                                  tGTTATCCGCTTTCGGTGCTTGAGATTCATTTAGCTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGAttt                                                                                     <  2:241435/141‑1 (MQ=255)
                                                             ttCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCa                                                                          >  1:267031/1‑141 (MQ=255)
                                                               cGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCAtt                                                                        <  2:250419/141‑1 (MQ=255)
                                                                   gCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTttgt                                                                    >  2:114096/1‑141 (MQ=255)
                                                                           ttCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCttt                                                            <  1:113420/141‑1 (MQ=255)
                                                                             cATTTAACTTTTTAGCTTCTGACACTAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTttgt                                                          >  1:80141/1‑141 (MQ=255)
                                                                                    cTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAAtt                                                   >  1:125898/1‑141 (MQ=255)
                                                                                        ttAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTa                                               >  1:276966/1‑141 (MQ=255)
                                                                                        ttAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCAGCGTTTAAGTTAGGGATATTCAAGCTTTCATAGAAAGCATTTTGTTGTTCGTGGTGGTATTTGttg                                               >  2:129576/1‑140 (MQ=17)
                                                                                           gCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTCGGCATCTTAACGCTCTCCTCGAAAGCATTTTGTTGTTATTTGTTGAAATTGTTATCa                                            >  1:130343/1‑141 (MQ=18)
                                                                                           gCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCa                                            >  1:326205/1‑141 (MQ=255)
                                                                                            cTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAg                                           >  1:153937/1‑141 (MQ=255)
                                                                                             ttCTGAAAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAAGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGc                                          <  1:143616/141‑1 (MQ=255)
                                                                                              tCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAGGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCt                                         >  1:251121/1‑141 (MQ=255)
                                                                                                     aTAGGTTAGCATATCGGTTTGGGTCATCTTTTAAGCTTTGGGTTAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATGCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTg                                  <  2:272057/141‑1 (MQ=17)
                                                                                                      tAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGc                                 >  1:117278/1‑141 (MQ=255)
                                                                                                          ttAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTga                             >  2:326754/1‑141 (MQ=255)
                                                                                                          ttAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTga                             >  1:232886/1‑141 (MQ=255)
                                                                                                              cgcTTTGGCTTGGGTCATCTTTTAGGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGTATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGAtt                         >  2:32659/3‑141 (MQ=255)
                                                                                                               aCTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTc                        >  2:51076/1‑141 (MQ=255)
                                                                                                                     gCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTCGTTTaa                  <  1:214847/141‑1 (MQ=255)
                                                                                                                      cTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTAGGTGCTTGAGTTTCGTTTAAt                 >  2:273945/1‑141 (MQ=255)
                                                                                                                       ttGGGTCATCTTTTAGGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGTATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGaa                                                      >  2:138439/1‑103 (MQ=9)
                                                                                                                       ttGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGaa                                                      <  1:138439/103‑1 (MQ=255)
                                                                                                                       ttGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTCGTTTAAtt                >  1:197402/1‑141 (MQ=255)
                                                                                                                         gggTCATCTTTTAGGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTGAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTCGTTTAAtttt              <  2:53723/141‑1 (MQ=255)
                                                                                                                          ggTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTCGTTTAAttttt             >  1:365093/1‑141 (MQ=255)
                                                                                                                            tCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTCGTTTAATTTTTTa           >  1:19605/1‑141 (MQ=255)
                                                                                                                             cATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTCGTTTAATTTTTTAg          >  2:365791/1‑141 (MQ=255)
                                                                                                                              aTCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATCTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTCGTTTCATTTTTTAgt         >  2:14192/1‑140 (MQ=255)
                                                                                                                                   ttAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTCGTTTAATTTTTTAGCTTCAc    <  1:168653/141‑1 (MQ=255)
                                                                                                                                     aaGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAACGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTCGTTTAATTTTTTAGCTTCCCCt  <  2:88824/141‑1 (MQ=255)
                                                                                                                                     aaGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTCGTTTGATGTGTTAGCTTCACCt  >  1:304254/1‑141 (MQ=255)
                                                                                                                                     aaGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTCGTTTAATTTTTTAGCTTCACCt  <  1:114601/141‑1 (MQ=255)
                                                                                                                                     aaGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCGTTGGGTGCTTGAGATTCGTTTAATTTTTTAGCTTctgct  >  2:301790/1‑137 (MQ=255)
                                                                                                                                      |                                                                                                                                           
AGGTAAATGTAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAATTTGTTATCCGCTTTCGGTGCTTGAGATTCATTTAACTTTTTAGCTTCTGACAATAGGTTAGCACTTTGGCTTGGGTCATCTTTTAAGCTTTGGATGAAACCATTGCGTTGTTCTTCGTTTAAGTTAGGCATATTCAAGATTTCATAGAAAGCATTTTGTTGTTCTTTGTTGAAATTGTTATCAGCTTTCGGTGCTTGAGATTCGTTTAATTTTTTAGCTTCACCT  >  CP000730/128895‑129168

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