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.4 / NA 29A155G (GCA→GGA) pepPM24B family Xaa‑Pro dipeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (14/14);  minor base T (1/0);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATTTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGAT  >  CP000730/1646049‑1646320
                                                                                                                                            |                                                                                                                                   
taatTTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTac                                                                                                                              <  1:171641/148‑1 (MQ=255)
taatTTTATCACTTGCAACACCATGCGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTac                                                                                                                              >  2:43244/1‑148 (MQ=255)
        tCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTa                                                                                                                      >  1:132600/1‑148 (MQ=255)
              gCAACACCATGTTGTAATGCACCTCTCTGTCCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatat                                                                                                                <  2:223982/148‑1 (MQ=255)
                aaCACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAatatat                                                                                                              >  1:31320/1‑148 (MQ=255)
                                           aCCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgct                                                                                   <  1:76285/148‑1 (MQ=255)
                                           aCCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAgct                                                                                   >  2:9563/1‑148 (MQ=255)
                                                    tACACTCTTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTACTATTGCCTTTAATTCTTTTCCAGACATGCCTCTCTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGa                                                                          <  2:81189/148‑1 (MQ=255)
                                                    tACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGa                                                                          <  2:400576/148‑1 (MQ=255)
                                                         tCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCATTCATTCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTaa                                                                     >  1:82935/1‑148 (MQ=255)
                                                         tCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTaa                                                                     >  2:327376/1‑148 (MQ=255)
                                                         tCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTaa                                                                     >  2:30896/1‑148 (MQ=255)
                                                         tCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTATCTAATATTGCCTTTAATTCTTTTACATTCATGCCCCCTTTTACAACACTTAAAATTTTTTAATATTTTTCCTCAACAATTTTGCTTGGTTTTTGAATTaa                                                                     >  2:60917/1‑148 (MQ=255)
                                                           gTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATATTGCCTTCTAATATTGCCATTAATTCTTTTTCAGTCGTGCCTCCTTTTACAACAGTTCAAATATATTCACATGTTTCACCAACAATATTCGCTGCTTTTTGAATTaaac                                                                   >  2:150113/1‑147 (MQ=255)
                                                                                    cTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGtc                                          <  2:31320/148‑1 (MQ=255)
                                                                                       aTTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTctct                                       >  2:16056/1‑148 (MQ=255)
                                                                                       aTTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTctct                                       >  2:59215/1‑148 (MQ=255)
                                                                                       aTTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTctct                                       >  2:144002/1‑148 (MQ=255)
                                                                                         tCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTTATTTTTTTTTCCTTTTACAATAGTTAAAATATATTTATATGTTTATTCAACAATATTATCTTCTTTTTCTATTAATTCAATTTCTTCACATTTTTTAACGTTTCTaa                                     >  2:359807/1‑148 (MQ=255)
                                                                                         tCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTaa                                     >  1:267106/1‑148 (MQ=255)
                                                                                             gCATTTTGCTTTCTATTATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAAtttt                                 <  2:135052/148‑1 (MQ=255)
                                                                                             gCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAAtttt                                 <  2:241783/148‑1 (MQ=255)
                                                                                                   tGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTCTCTAc                           <  2:127384/148‑1 (MQ=255)
                                                                                                               ttGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAAt               <  2:384429/148‑1 (MQ=255)
                                                                                                               ttGCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAAt               <  1:247548/148‑1 (MQ=255)
                                                                                                                 gCCTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGc             <  2:235912/148‑1 (MQ=255)
                                                                                                                   cTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCtt           <  2:255735/148‑1 (MQ=255)
                                                                                                                   cTTTAATTCTTTTTCAGTCATGCCTCCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCtt           <  2:281572/148‑1 (MQ=255)
                                                                                                                            tttttCAGTCATTCCTTCTTTTACAACAGTTAAAATATATTCATATGTTTCATCATCAATATTAGCTGCTTTTTGAATTAATTCAATTTCTTCAGCATATTTGACGTCTCTAATTTTATCTATAGTATTAGAAATGATTATTAATGat  >  2:164593/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                   
TAATTTTATCACTTGCAACACCATGTGGTAATGCACCTCTATGACCAGATGCTACAATCGTATCGAATGATGGTCCATCTGCTCCTAATTCTAGCATTTTGCTTTCTAATATTGCCTTTAATTCTTTTTCAGTCATGCCTGCTTTTACAACAGTTAAAATATATTCATATGTTTCATCAACAATATTAGCTGCTTTTTGAATTAAAGCAATTTCGTCAGCATCTTTGACGTCTCTAATTTTATCTACAGTATTAGAAATGCTTATTAATGAT  >  CP000730/1646049‑1646320

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