Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,988,845 A→G 11.7% H32R (CAT→CGT)  PP_1780 → mannosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,988,8450AG11.7% 86.5 / 2.0 34H32R (CAT→CGT) PP_1780mannosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (15/15);  new base G (2/2);  total (17/17)
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

AAAGCCCCAATTTGTCGGCATGGCTTGGGTTCCAACAGCTGTGGAGTGCATCCGGTGAAAAAACCTTTCTATCGCGCTTTCGAGGATCGCTACCGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGACCTGGATGACGGCATGCTGGAGCAATGCCGGATCGAGG  >  NC_002947/1988697‑1988991
                                                                                                                                                    |                                                                                                                                                  
aaaGCCCCAATTTGTCGGCATGGCTTGGGTTCCAACAGCTGTGGAGTGCATCCGGTGAAAAAACCTTTCTATCGCGCTTTCGAGGATCGCTACCGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCa                                                                                                                                                    >  3:21759/1‑149 (MQ=255)
      caaaattGTCGGCATGGCTTGTGTTACAACAGCTGTGGAGTGCATCCGGTGAAAAAACCTATCTATCGCGCTTTCCAGGATCGCTACCGTGGCTCCACCAAACAGATCACCGACCGGCTCCAGGTCTATATGCCGTTCCTTCATCAGc                                                                                                                                               >  2:75271/6‑148 (MQ=255)
                                 aaCAGGTGTGGAGTGGATCCGGTGAAAAAACCTTTCTTTCGCGGTTTCGAGGATGGCTTCCGTGGCTCCACGTTACTGATCTCCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCTCTGTTGGCACTGGATGACGAGCATCGTGc                                                                                                                   <  1:75271/149‑1 (MQ=255)
                                 aaCAGCTGTGGAGTGCATCCGGTGAAAAAACCTTTCTATCGCGCTTTCGAGGATCGCTACCGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGc                                                                                                                   >  4:138457/1‑149 (MQ=255)
                                    aGCTGTGGAGTGCATCCGGTGAAAAAACCTTTCTATCGCGCTTTCGAGGATCGCTACCGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGAt                                                                                                                >  6:44318/1‑149 (MQ=255)
                                       tgtgGAGTGCATCCGGTGAAAAAACCTTTCTATCGCGCTTTCGAGGATCGCTACCGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGa                                                                                                             >  3:180223/1‑149 (MQ=255)
                                             gTGCATCCGGTGAAAAAACCTTTCTATCGCGCTTTCGAGGATCGCTACCGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTggg                                                                                                       <  3:138457/149‑1 (MQ=255)
                                             gTGCATCCGGTGAAAAAACCTTTCTATCGCGCTTTCGAGGATCGCTACCGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATAGACTTggg                                                                                                       <  3:192719/149‑1 (MQ=255)
                                                              aCCTTTCTATCGCGCTTTCGAGGATCGCTACCGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGt                                                                                      >  3:297905/1‑149 (MQ=255)
                                                                    cTATCGCGCTTTCGAGGATCGCTACCGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTg                                                                                 <  4:21759/148‑1 (MQ=255)
                                                                           gcTTTCGAGGATCGCTACCGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACt                                                                         <  3:270161/149‑1 (MQ=255)
                                                                                        gCTACCGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGc                                                            >  5:276298/1‑149 (MQ=255)
                                                                                          tACCGTGGCTCCACTAAACTGATCACCGAGCGGCTGCAGGTTTATATGCCGTTCCATCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAggg                                                             <  3:223366/146‑1 (MQ=255)
                                                                                            ccGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATg                                                                                                                                >  2:266632/1‑77 (MQ=255)
                                                                                            ccGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATg                                                                                                                                <  1:266632/77‑1 (MQ=255)
                                                                                               tGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGgt                                                                              >  5:101653/1‑124 (MQ=255)
                                                                                               tGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGgt                                                                              <  6:101653/124‑1 (MQ=255)
                                                                                                 gCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCAATGCa                                                   <  4:297905/149‑1 (MQ=255)
                                                                                                        gAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCa                                                        <  6:219091/137‑1 (MQ=255)
                                                                                                        gAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCa                                                        >  5:219091/1‑137 (MQ=255)
                                                                                                             tGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCGTCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGAc                                       >  5:163969/1‑149 (MQ=255)
                                                                                                             tGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCGTCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGAc                                       >  1:268111/1‑149 (MQ=255)
                                                                                                             tGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGtt                                                                                  >  4:194237/1‑106 (MQ=255)
                                                                                                             tGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGtt                                                                                  <  3:194237/106‑1 (MQ=255)
                                                                                                               aTCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGACCt                                     <  4:180223/149‑1 (MQ=255)
                                                                                                                    cGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGACCTGGAt                                 >  5:26579/1‑148 (MQ=255)
                                                                                                                       gCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGACCTgg                                   >  6:295246/1‑143 (MQ=255)
                                                                                                                       gCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCAGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGACCTgg                                   <  5:295246/143‑1 (MQ=255)
                                                                                                                              cAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGACCTGGATGACGGCATGc                       <  6:26579/148‑1 (MQ=255)
                                                                                                                                      tatGCCGTTCCTTCATCAGCTGTTGGCACTGGATGACCAGCATCTTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCATCAATGAACTTCGCGTAGACATGGa                                  >  1:302474/1‑129 (MQ=255)
                                                                                                                                          ccGTTCCTTCGTCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTTCTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGACCGGGATGACGGCATGCTGGAGCAATGcc           <  2:268111/148‑1 (MQ=255)
                                                                                                                                          ccGTTCCTTCGTCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGACCTGGATGACGGCATGCTGGAGCAATGcc           <  6:163969/148‑1 (MQ=255)
                                                                                                                                           cGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGACCTGGATGACGGCATGCTGGAGCAATGCCg          >  4:299741/1‑148 (MQ=255)
                                                                                                                                               ccTTCATCCGCTGTTGGCACTGGATGACGTGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGTGTGTACTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGACCTGGATGACGGCATGCTGGAGCAATGCCGGATc      <  7:2607/148‑1 (MQ=255)
                                                                                                                                                  tCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGACCTGGATGACGGCATGCTGGAGCAATGCCGGATCGAgg  >  7:284832/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AAAGCCCCAATTTGTCGGCATGGCTTGGGTTCCAACAGCTGTGGAGTGCATCCGGTGAAAAAACCTTTCTATCGCGCTTTCGAGGATCGCTACCGTGGCTCCACGAAACTGATCACCGAGCGGCTGCAGGTGTATATGCCGTTCCTTCATCCGCTGTTGGCACTGGATGACGAGCATCGTGCGATTGACTTGGGCTGTGGCCGTGGCGAGTGGTTGGGTGTACTGGCCGGGGAGGGCTTCAATGCAGTTGGCGTAGACCTGGATGACGGCATGCTGGAGCAATGCCGGATCGAGG  >  NC_002947/1988697‑1988991

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