Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,646,189 G→C A155G (GCA→GGA)  pepP ← M24B family Xaa‑Pro dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,646,1890GC96.5% 90.7 / NA 29A155G (GCA→GGA) pepPM24B family Xaa‑Pro dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (11/17);  minor base T (1/0);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATA  >  CP000730/1646056‑1646321
                                                                                                                                     |                                                                                                                                    
aTCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTac                                                                                                                               >  2:179802/1‑141 (MQ=255)
  cACTTGCATCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaa                                                                                                                             <  2:278358/141‑1 (MQ=255)
  cACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTGATATTGCCTTTAATTCTTTTTGAGTCGTGCCTCCTTTTacta                                                                                                                             <  2:104052/141‑3 (MQ=255)
           cacCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTaaaa                                                                                                                    <  1:174641/141‑1 (MQ=255)
              cATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaac                                                                                                                            >  2:76916/1‑130 (MQ=255)
              cATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTacaac                                                                                                                            <  1:76916/130‑1 (MQ=255)
                        tGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGt                                                                                                       >  2:182854/1‑141 (MQ=255)
                         gCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGtt                                                                                                      <  1:272067/141‑1 (MQ=255)
                              tctATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTtcatc                                                                                                 <  1:15730/141‑1 (MQ=255)
                                tATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATcaa                                                                                               <  2:140537/141‑1 (MQ=255)
                                   gACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTAATATGTTTCATcaacaa                                                                                            >  2:230028/1‑141 (MQ=255)
                                         aTGCTACACTCGTATCGAATGATGGTCCATCTGCTTCTAATTCTAGCATTTTGCTTTCTAATATTGCCGTTAATTCTTTTACAGTCATGCCTCCTTTTCCAACAGTTAAAATATCTGCATATGTTTCATCTACAATATTAg                                                                                      <  2:285348/141‑1 (MQ=255)
                                         aTGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAg                                                                                      <  1:182854/141‑1 (MQ=255)
                                           gCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgct                                                                                    <  1:145328/141‑1 (MQ=255)
                                                         gAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTaa                                                                      >  2:12313/1‑141 (MQ=255)
                                                         gAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTaa                                                                      <  2:324453/141‑1 (MQ=255)
                                                                      tCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCa                                                         <  2:117932/141‑1 (MQ=255)
                                                                           tCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATc                                                    <  2:16922/141‑1 (MQ=255)
                                                                             cTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTTTTTTTTTTCTTCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTTAATTAAAGCAATTTCGTCAGCATCtt                                                  >  2:229150/1‑141 (MQ=255)
                                                                                     aGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtct                                          >  1:17515/1‑141 (MQ=255)
                                                                                       cATTTTGCTTTCTAATATTTCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAGCAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTctct                                        <  2:152750/141‑1 (MQ=255)
                                                                                                 tCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGc                                                       <  2:60858/116‑1 (MQ=255)
                                                                                                 tCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGc                                                       >  1:60858/1‑116 (MQ=255)
                                                                                                                aaTTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAg                                                                                      >  2:9664/1‑70 (MQ=255)
                                                                                                                aaTTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAg                                                                                      <  1:9664/70‑1 (MQ=255)
                                                                                                                     tttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgc                                                                                  <  1:163376/69‑1 (MQ=255)
                                                                                                                     tttttCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgctgc                                                                                  >  2:163376/1‑69 (MQ=255)
                                                                                                                           aGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTCGCTGCTTTTTTAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGa    >  2:164936/1‑141 (MQ=255)
                                                                                                                             tCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGata  >  1:328843/1‑141 (MQ=255)
                                                                                                                                     |                                                                                                                                    
ATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGATA  >  CP000730/1646056‑1646321

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