Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 5,107,530 G→A 8.2% V25V (GTC→GTT)  scpA ← Segregation and condensation protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1435,107,5300GA 8.2% 160.3 / 3.8 49V25V (GTC→GTT) scpASegregation and condensation protein A
Reads supporting (aligned to +/- strand):  ref base G (30/15);  new base A (2/2);  total (32/17)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GATATCGATATTCTGCTTGCGGATCAGGTACAGCAACAGGTCCAGCGGGCCTTCGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGTCGCCACAACCTTCGCGGATCAGTTCCGGTTCG  >  ACB143/5107398‑5107661
                                                                                                                                    |                                                                                                                                   
gATATCGATATTCTGCTTGCGGATCAGGTACAGCAACAGGTCCAGCGGGCCTTCGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGAcc                                                                                                                                  <  4:38892/136‑1 (MQ=255)
     cGATATTCTGCTTGCGGATCAGGTACAGCAACAGGTCCAGCGGGCCTTCGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAgcgc                                                                                                                             >  6:149622/1‑136 (MQ=255)
             tGCTTGCGGATCAGGTACAGCAACAGGTCCAGCGGGCCTTCGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAAATg                                                                                                                       >  7:238164/1‑134 (MQ=255)
                  gCGGATCAGGTACAGCAACAGGTCCAGCGGGCCTTCGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGTCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGc                                                                                                                >  4:166647/1‑136 (MQ=255)
                  gCGGATCAGGTACAGCAACAGGTCCAGCGGGCCTTCGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGc                                                                                                                >  5:214031/1‑136 (MQ=255)
                   cGGATCAGGTACAGCAACAGGTCCAGCGGGCCTTCGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGcc                                                                                                               <  5:267246/136‑1 (MQ=255)
                   cGGATCAGGTACAGCAACAGGGCCAGCGGGCCTTCGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGc                                                                                                                <  5:65596/135‑1 (MQ=255)
                                   aCAGGTCCAGCGGGCCTTCGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCg                                                                                               >  8:106068/1‑136 (MQ=255)
                                       gTCCAGCGGGCCTTCGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTcg                                                                                           >  7:113286/1‑136 (MQ=255)
                                               ggCCTTCGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCg                                                                                                              >  4:3176/1‑109 (MQ=255)
                                               ggCCTTCGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCg                                                                                                              <  3:3176/109‑1 (MQ=255)
                                                     cGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACATGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGACAACTGCAGTTTCCGAGGCTCCTCGCCCCGTACGGCCTGCGCCACAgg                                                                             >  2:288341/1‑136 (MQ=255)
                                                           cTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTc                                                                       <  3:306105/136‑1 (MQ=255)
                                                                 ggATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTg                                                                 >  1:93585/1‑136 (MQ=255)
                                                                       cTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCAc                                                           <  1:274389/136‑1 (MQ=255)
                                                                            aGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGtt                                                      <  3:69088/136‑1 (MQ=255)
                                                                            aGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGtt                                                      <  3:57693/136‑1 (MQ=255)
                                                                            aGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGtt                                                      <  1:155809/136‑1 (MQ=255)
                                                                                cGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAAACCAGCGCCAACTGCAGTTGc                                                                                                                >  7:81050/1‑74 (MQ=255)
                                                                                cGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAAACCAGCGCCAACTGCAGTTGc                                                                                                                <  8:81050/74‑1 (MQ=255)
                                                                                 gTCCGGCGGTATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCa                                                 <  3:167156/136‑1 (MQ=255)
                                                                                  tCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAg                                                >  2:162094/1‑136 (MQ=255)
                                                                                   ccGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAg                                                >  6:248888/1‑135 (MQ=255)
                                                                                     gtcggAATGTACAGGTCCAGAGTCATATCGTTAATGGCATCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCAAGCTGGCTCACGCGTTGACCAGGaa                                             >  8:207192/3‑136 (MQ=255)
                                                                                     ggcggAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGGGGCTCCTCGACCCGTACGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGaa                                             >  4:165085/1‑136 (MQ=255)
                                                                                     ggcggAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGaa                                             >  7:297360/1‑136 (MQ=255)
                                                                                     ggcggAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGaa                                             >  6:151342/1‑136 (MQ=255)
                                                                                     ggcggAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGaa                                             >  5:72100/1‑136 (MQ=255)
                                                                                     ggcggAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCGCAGGCGCTTCCCGCTGGCTCACGCGTGGACCAggca                                             >  2:46434/1‑134 (MQ=255)
                                                                                                    tCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGTCGc                               >  1:294507/1‑135 (MQ=255)
                                                                                                    tCCAGCGGCAGCTCGGTCAGGGCTTCACCATAAACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCAc                                                           <  1:152626/107‑1 (MQ=255)
                                                                                                    tCCAGCGGCAGCTCGGTCAGGGCTTCACCATAAACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCAc                                                           >  2:152626/1‑107 (MQ=255)
                                                                                                      cAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGTCGCca                             >  7:174016/1‑135 (MQ=255)
                                                                                                      cAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCACGCGCTACCCCCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGGCGCca                             >  8:46042/1‑135 (MQ=255)
                                                                                                           gCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGTCGCCACAACCt                       >  6:135305/1‑136 (MQ=255)
                                                                                                           gCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGTCGCCACAACCt                       >  4:296981/1‑136 (MQ=255)
                                                                                                           gCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGTCGCCACAACCt                       >  3:165648/1‑136 (MQ=255)
                                                                                                           gCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGTCGCCACAACCt                       >  3:118260/1‑136 (MQ=255)
                                                                                                                tCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGtt                                                      <  3:216184/100‑1 (MQ=255)
                                                                                                                tCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGtt                                                      >  4:216184/1‑100 (MQ=255)
                                                                                                                tCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGTCGCCACAACCTTcgcg                  >  2:96213/1‑136 (MQ=255)
                                                                                                                tCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGCAAAGGCGTGGGGTCGCCACAACCTTcgcg                  >  5:281793/1‑136 (MQ=255)
                                                                                                                  ggTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGTCGCCACAACCTTcgcg                  <  3:49754/134‑1 (MQ=255)
                                                                                                                  ggTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGGCGCCACAACCTTcgcg                  <  8:263936/134‑1 (MQ=255)
                                                                                                                            tCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGTCGCCACAACCTTCGCGGATCAGTTCCgg      >  6:298621/1‑136 (MQ=255)
                                                                                                                             cACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGTCGCCACAACCTTCGCGGATCAGTTCCg       <  3:200405/134‑1 (MQ=255)
                                                                                                                                cATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTgg                                     >  6:183230/1‑101 (MQ=255)
                                                                                                                                cATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTgg                                     <  5:183230/101‑1 (MQ=255)
                                                                                                                                cATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGTCGCCACAACCTTCGCGGATCAGTTCCGGTtcg  >  3:273188/1‑136 (MQ=255)
                                                                                                                                    |                                                                                                                                   
GATATCGATATTCTGCTTGCGGATCAGGTACAGCAACAGGTCCAGCGGGCCTTCGAAAGCTTCGAGGATGACTTCAAGGGCGTCCGGCGGAATGTACAGGTCCAGCGGCAGCTCGGTCAGGGCTTCACCATAGACCAGCGCCAACTGCAGTTGCCGAGGCTCCTCGACCCGTTCGGCCGGCGCCTCAGGCGCTTCCAGCTGGCTCACGCGTTGACCAGGAATGGCGTGGGGTCGCCACAACCTTCGCGGATCAGTTCCGGTTCG  >  ACB143/5107398‑5107661

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