Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,474,639 C→G A9129P (GCA→CCA)  ebh ← extracellular matrix binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,474,6390CG97.5% 136.3 / NA 40A9129P (GCA→CCA) ebhextracellular matrix binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (16/23);  minor base A (1/0);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACCTTTAAATCCATCTTGTTGAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGCTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATCTACAGTTTGCTCAACGTTAGCATGTTCAT  >  CP000730/1474504‑1474773
                                                                                                                                       |                                                                                                                                      
gACCTTTAAATCCATCTTGTTGAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGAt                                                                                                                                   <  1:371641/141‑1 (MQ=255)
gACCTTTAAATCCATCTTGTTGAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGAt                                                                                                                                   <  1:154315/141‑1 (MQ=255)
        aaTCCATCTTGTTGAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACttgt                                                                                                                           <  2:361239/141‑1 (MQ=255)
          tCCATCTTGTTGAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTACTGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACttgttg                                                                                                                         >  2:233359/1‑141 (MQ=255)
                     gAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCAt                                                                                                              <  1:115674/141‑1 (MQ=255)
                      agagttGTTTAAGGGATTCAATTGGGCGATGTCATGTTTACCTTTTGTTTTGGCTAGTGCTAGTCGTTAATCTCGATTTATTGCTTTTTTTGCATTTAAAATATTTTGAAGTAGTGGTTCCTCTTTTTGCTTATTCGCAtt                                                                                                             <  2:192908/137‑1 (MQ=255)
                              ttAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCa                                                                                                     >  1:71944/1‑141 (MQ=255)
                              ttAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCa                                                                                                     >  2:170662/1‑141 (MQ=255)
                                 aTGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATa                                                                                                  >  1:130178/1‑141 (MQ=255)
                                    cATTCAATTGGTTGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACAt                                                                                               <  2:46951/141‑1 (MQ=255)
                                    cATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACAt                                                                                               <  1:226618/141‑1 (MQ=255)
                                     aTTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACAtt                                                                                              <  2:205987/141‑1 (MQ=255)
                                                       aTGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGc                                                                            >  1:405527/1‑141 (MQ=255)
                                                        tGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCt                                                                           <  1:283788/141‑1 (MQ=255)
                                                        tGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCt                                                                           <  1:9831/141‑1 (MQ=255)
                                                                 tttGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATaa                                                                  >  2:181116/1‑141 (MQ=255)
                                                                     ttttAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCAt                                                              >  2:320993/1‑141 (MQ=255)
                                                                      tttAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCAtt                                                             >  1:6296/1‑141 (MQ=255)
                                                                      tttAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCAtt                                                             <  1:224360/141‑1 (MQ=255)
                                                                       ttAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCAttt                                                            <  1:397549/141‑1 (MQ=255)
                                                                         aGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTg                                                          <  1:166697/141‑1 (MQ=255)
                                                                                gCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTTGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCACCATCAGCAATAGCTTGTTTATAAGCATTTTATTTATCt                                                   >  2:243869/1‑141 (MQ=255)
                                                                                gCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCt                                                   >  2:167645/1‑141 (MQ=255)
                                                                                     acGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTCTTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATc                                              <  1:17286/141‑1 (MQ=255)
                                                                                     acGTTCATCACCATTTAATGCTTGTATTGCATTTAAAATATTTTGAAGTGATTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGATTGTTTATAAGCATTTTGTATATCTGAATc                                              >  2:162970/1‑141 (MQ=255)
                                                                                        ttCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGc                                           <  1:31824/141‑1 (MQ=255)
                                                                                          catcaCCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCtt                                         >  2:126775/1‑141 (MQ=255)
                                                                                              tccATTTAATGATTGTTTTGCAGTTAAAATATTTTGAAGTGGTTGAGCCACTTGTTGCTTATTCTCATTTTATTTCAATACATTTTCAGCATCAGCATTAGCCTGCTTATCGGCATTTTGTCTATCTGGATCTGCTtgtgt                                     <  2:294062/140‑1 (MQ=255)
                                                                                              aCCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATCAATTTTCAGCATCAGTAGTCGCTTGTTTATAAGAATTTTGTTTATCTGAATCTGCTtgtgt                                     >  1:349540/1‑141 (MQ=255)
                                                                                               ccATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTtgtgtg                                    <  1:102345/141‑1 (MQ=255)
                                                                                                cATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACGTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTGgtgtgt                                   >  1:398485/1‑141 (MQ=255)
                                                                                                    tAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATc                               >  2:314874/1‑141 (MQ=255)
                                                                                                        gCTTGTTTTGCATTTAAAATATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATCTACa                           >  1:34315/1‑141 (MQ=255)
                                                                                                                         aaTATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATCTACAGTTTGCTCAACGTTAGc          <  1:388357/141‑1 (MQ=255)
                                                                                                                         aaTATTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATCTACAGTTTGCTCAACGTTAGc          <  2:29563/141‑1 (MQ=255)
                                                                                                                            atTTTGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATCTACAGTTTGCTCAACGTTAGCATg       <  2:128507/141‑1 (MQ=255)
                                                                                                                             ttttGAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTTAATCTGCTTGTGTGTAATCTACAGTTTGCTCAACGTTAGCATGt      <  2:203087/141‑1 (MQ=255)
                                                                                                                                 gAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATCTACAGTTTGCTCAACGTTAGCATGTTCAt  <  2:300149/141‑1 (MQ=255)
                                                                                                                                 gAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATCTACAGTTTGCTCAACGTTAGCATGTTCAt  >  2:186273/1‑141 (MQ=255)
                                                                                                                                 gAAGTGGTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATCTACAGTTTGCTCAACGTTAGCATGTTCAt  <  2:189976/141‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GACCTTTAAATCCATCTTGTTGAGCATTGTTTAATGCATTCAATTGGTCGATGTCATGTTTACCATTTGTTTTAGCAAGTGCTACACGTTCATCACCATTTAATGCTTGTTTTGCATTTAAAATATTTTGAAGTGCTTGATCCACTTGTTGCTTATTCGCATTTTGTTTCAATACATTTTCAGCATCAGCAATAGCTTGTTTATAAGCATTTTGTTTATCTGAATCTGCTTGTGTGTAATCTACAGTTTGCTCAACGTTAGCATGTTCAT  >  CP000730/1474504‑1474773

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