Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,189,989 A→T 13.6% F221I (TTC→ATC) ‡ PP_2799 ← class III aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,189,9890AT13.6% 107.7 / 9.4 44F221I (TTC→ATC) ‡PP_2799class III aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (18/20);  new base T (4/2);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

ATTGTGCCAAGCCTGCCAAAGCCGGTGATGATTTCGTCGATCACCACGAGGATGTCGTAGCGTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTTCGGGCAGGTCAGGTGATGAACCGGGATGGCCGGC  >  NC_002947/3189844‑3190134
                                                                                                                                                   |                                                                                                                                                 
aTTGTGCCAAGCCTGCCAAAGCCGGTGATGATTTCGTCGATCACCACGAGGATGTCGTAGCGTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAAcgc                                                                                                                                                <  7:43165/149‑1 (MQ=255)
               ccAAAGCCGGTGATGATTTCGTCGATCACCACGAGGATGTCGTAGCGTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGCCGCTGGTTTCCgg                                                                                                                                 >  8:112261/1‑149 (MQ=255)
                      cGGTGATGATTTCGTCGATCACCACGAGGATGTCGTAGCGTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCg                                                                                                                          >  2:193459/1‑149 (MQ=255)
                           atgatTTCGTCGATCACCACGAGGATGTCGTAGCGTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAg                                                                                                                     >  5:6053/1‑149 (MQ=255)
                            tgatTTCGTCGATCACCACGAGGATGTCGTAGCGTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACGCCCCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCaa                                                                                                                      <  1:92644/149‑1 (MQ=255)
                             gatTTCGTCGATCACCACGAGGATGTCGTAGCGTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGTTCATCGGCGAGCCGGTAATCgcggcagg                                                                                                                                          >  6:65518/1‑123 (MQ=37)
                             gatTTCGTCGATCACCACGAGGATGTCGTAGCGTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGTTCATCGGCGAGCCGGTAATCgcggcagg                                                                                                                                          <  5:65518/126‑4 (MQ=37)
                                      gATCACCACGAGGATGTCGTAGCGTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGtt                                                                                                          >  5:290956/1‑149 (MQ=255)
                                           ccacGAGGATGTCGTAGCGTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAgt                                                                                                     <  4:305709/149‑1 (MQ=255)
                                                  gATGTCGTAGCGTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGttcgttcg                                                                                              >  4:103833/1‑149 (MQ=255)
                                                    tGTCGTAGCGTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGcc                                                                                            <  1:193459/149‑1 (MQ=255)
                                                             gTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGAcc                                                                                   >  5:163378/1‑149 (MQ=255)
                                                              tttGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCg                                                                                  <  6:6053/149‑1 (MQ=255)
                                                                tGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAATGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCAAGACGGACCGt                                                                                 >  1:147758/1‑148 (MQ=255)
                                                                tGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATAGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGt                                                                                 >  4:114995/1‑148 (MQ=255)
                                                                   acacCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAg                                                                             <  8:271114/149‑1 (MQ=255)
                                                                     acCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGTGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGTTTTCCGGCCCTTCGGCAAGTATATAGCGTTCCAGTTCGTTCGCCAGACGTACCGTAAAGGc                                                                           <  2:17175/149‑1 (MQ=255)
                                                                     acCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAgttc                                                                                                   <  4:209717/125‑1 (MQ=255)
                                                                     acCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAgttc                                                                                                   >  3:209717/1‑125 (MQ=255)
                                                                                cgccgcCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGAt                                                                 >  4:16269/1‑148 (MQ=255)
                                                                                      ccAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTg                                                          <  2:222319/149‑1 (MQ=255)
                                                                                      ccAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTg                                                          <  2:50258/149‑1 (MQ=255)
                                                                                               tgtgGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAg                                                                             <  8:43619/121‑1 (MQ=255)
                                                                                               tgtgGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAg                                                                             >  7:43619/1‑121 (MQ=255)
                                                                                               tgtgGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCa                                                  >  2:258059/1‑148 (MQ=255)
                                                                                                         tggATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCTATGGTTTCCGGCCCTTCGGCAATGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCTGTAGaa                                       <  7:182241/148‑1 (MQ=255)
                                                                                                          ggATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAg                                      >  6:276803/1‑149 (MQ=255)
                                                                                                           gATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGt                                     >  2:176829/1‑149 (MQ=255)
                                                                                                               aCG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTTCgg                                 >  8:112638/1‑149 (MQ=255)
                                                                                                               aCG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTTCg                                  <  4:242641/148‑1 (MQ=255)
                                                                                                                cG‑‑CCACCTGCCGCGATTACCTGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCg                                                                   <  5:563/114‑1 (MQ=255)
                                                                                                                cG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCg                                                                   >  6:563/1‑114 (MQ=255)
                                                                                                                cG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTACggg                                <  3:16269/149‑1 (MQ=255)
                                                                                                                     ccatcgCCGCGTTCATCGGCGAGCCGGTAATCGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTTCGGGcag                             >  3:163606/6‑149 (MQ=38)
                                                                                                                      catcgCCGCGTTCATCGGCGAGCCGGTAATCGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTTCGGGc                               >  4:247031/5‑146 (MQ=38)
                                                                                                                      catcgCCGCGTTCATCGGCGAGCCGGTAATCGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTTCGGGc                               <  3:247031/142‑1 (MQ=38)
                                                                                                                       atcgCCGCGTTCATCGGCGAGCCGGTAATCGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTTCGGGcaggt                           >  3:144539/4‑149 (MQ=38)
                                                                                                                           ccGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTTCGGGcaggtcag                        <  7:112638/148‑1 (MQ=255)
                                                                                                                                  tACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCtt                                                                      <  3:154388/95‑1 (MQ=255)
                                                                                                                                  tACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCtt                                                                      >  4:154388/1‑95 (MQ=255)
                                                                                                                                       gCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTTCgg                                 <  2:47611/127‑1 (MQ=255)
                                                                                                                                       gCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTTCgg                                 >  1:47611/1‑127 (MQ=255)
                                                                                                                                           gCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTTCGGGCAGGTCAGGTGATGAACCGGGATg        <  5:276803/148‑1 (MQ=255)
                                                                                                                                             cGATGAATGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGAGGGACCGTAAAGGGCTCTAGCGATTCTCCTGGCCGAGCAAAGCGGTAGAAGTTCGGGCAGGTCAGGTGATGAACAGGGATGGcc     <  2:147758/149‑1 (MQ=255)
                                                                                                                                                  aaCGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTTCGGGCAGGTCAGGTGATGAACCGGGATGGCCGGc  <  6:290956/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATTGTGCCAAGCCTGCCAAAGCCGGTGATGATTTCGTCGATCACCACGAGGATGTCGTAGCGTTTGCACACCGCCTGGATCGCCGCCCAGTAGCCTGTGGGAGGAGGGATCACG‑‑CCACCTGCCGCGATTACCGGCTCGCCGATGAACGCGGCGATGGTTTCCGGCCCTTCGGCAAGGATATAGCGTTCCAGTTCGTTCGCCAGACGGACCGTAAAGGCCTCTTGCGATTCTCCTGGCCGAGCAAATCGGTAGAAGTTCGGGCAGGTCAGGTGATGAACCGGGATGGCCGGC  >  NC_002947/3189844‑3190134

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