Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,003,070 G→C 10.2% V990L (GTG→CTG)  PP_1786 → glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,003,0700GC10.2% 135.7 / 4.9 47V990L (GTG→CTG) PP_1786glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (21/21);  new base C (3/2);  total (24/23)
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

TAGCTAATGGCGTCTCCAGTGATATTCCGTGGACGCAGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTCGCGTAAGCCAGAAAGGAATCGTAGACCTTTAT  >  NC_002947/2002945‑2003213
                                                                                                                             |                                                                                                                                               
tAGCTAATGGCGTCTCCAGTGATATTCCGTGGACGCAGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAAcc                                                                                                                         <  1:172759/150‑1 (MQ=255)
      aTGGCGTCTCCAGTGATATTCCGTGGACGCAGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGc                                                                                                                                              <  2:21426/123‑1 (MQ=255)
      aTGGCGTCTCCAGTGATATTCCGTGGACGCAGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGc                                                                                                                                              >  1:21426/1‑123 (MQ=255)
       tGGCGTCTCCAGTGATATTCCGTGGACGCAGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAAc                                                                                                                          <  2:180345/142‑1 (MQ=255)
       tGGCGTCTCCAGTGATATTCCGTGGACGCAGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTGCAGAGTCAGCCTGAGAAc                                                                                                                          >  1:180345/1‑142 (MQ=255)
            tctcCAGTGATATTCCGTGGACGCAGAAAGATACCTCTTACACCAAACGGGTTATATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCaga                                                                                                                                        >  2:842991/1‑123 (MQ=255)
            tctcCAGTGATATTCCGTGGACGCAGAAAGATACCTCTTACACCAAACGGGTTATATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCaga                                                                                                                                        <  1:842991/123‑1 (MQ=255)
                 aGTGATATTCCGTGGACGCAGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGGATCATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGAttgtt                                                                                                        >  2:312625/1‑150 (MQ=255)
                        ttCCGTGGACGCAGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCAtt                                                                                                 >  1:832914/1‑150 (MQ=255)
                         tCCGTGGACGCAGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTc                                                                                                <  1:748734/150‑1 (MQ=255)
                            gTGGACGCAGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCaa                                                                                             <  2:575823/150‑1 (MQ=255)
                                    aGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAg                                                                                                                                                 >  2:341391/1‑90 (MQ=255)
                                    aGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAg                                                                                                                                                 <  1:341391/90‑1 (MQ=255)
                                    aGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCa                                                                                                                                    <  1:179717/103‑1 (MQ=255)
                                    aGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCa                                                                                                                                    >  2:179717/1‑103 (MQ=255)
                                     gAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAggg                                                                                    <  1:602641/150‑1 (MQ=255)
                                              ctctTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATCATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAAca                                                                           <  1:312625/150‑1 (MQ=255)
                                                     acCAAACGGGTTGTATGCGGGCATAACGGCGGAATGCCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCCCTGCTTCAGAGTCAGCCTGAGAAACTTGAGTCTTTGATTGTTGACCATTCAAAGGAAGAGGGATATGAACATCATtct                                                                    >  1:747540/1‑150 (MQ=255)
                                                       cAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCa                                                                  >  1:295309/1‑150 (MQ=255)
                                                        aaaCGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGt                                                                                                                                      >  2:833530/1‑81 (MQ=255)
                                                        aaaCGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGt                                                                                                                                      <  1:833530/81‑1 (MQ=255)
                                                              gTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTggggg                                                           >  2:536693/1‑150 (MQ=255)
                                                                  tATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAg                                                       >  2:344773/1‑150 (MQ=255)
                                                                    tGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGAc                                                                                                     >  2:474162/1‑102 (MQ=255)
                                                                    tGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGAc                                                                                                     <  1:474162/102‑1 (MQ=255)
                                                                           atatCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAACGTGGCAGTgg                                              >  2:48590/1‑149 (MQ=255)
                                                                             atCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCAATAAACGGAATAGGGGAATGAACATCATTCGCAGTGGGGGACAGTGCCGGTGAcc                                            >  2:663645/1‑150 (MQ=255)
                                                                                         tCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTcgc                                <  2:832914/150‑1 (MQ=255)
                                                                                         tCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTcgc                                <  2:295309/150‑1 (MQ=255)
                                                                                            gAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTCGCGTa                             <  1:536693/150‑1 (MQ=255)
                                                                                                 caAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTCGCGTAAGCCa                        >  1:422547/1‑150 (MQ=255)
                                                                                                 caAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTCGCGTAAGCCa                        <  1:344773/150‑1 (MQ=255)
                                                                                                        tattatCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGAtt                                                                                                           <  2:146257/60‑1 (MQ=255)
                                                                                                        tattatCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGAtt                                                                                                           >  1:146257/1‑60 (MQ=255)
                                                                                                        tattatCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTCGCGTAAGCCAGAAAGGa                 >  2:360480/1‑150 (MQ=255)
                                                                                                               tGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTCGCGTAAGCCAGAAAGGAATCGTAg          >  1:737595/1‑150 (MQ=255)
                                                                                                                  ttAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTCGCGTAAGCCAGAAAGGa                 >  2:299224/1‑140 (MQ=255)
                                                                                                                  ttAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTCGCGTAAGCCAGAAAGGa                 <  1:299224/140‑1 (MQ=255)
                                                                                                                    tctgAAGCACTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAAcatca                                                                        <  1:621744/80‑1 (MQ=255)
                                                                                                                    tctgAAGCACTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAAcatca                                                                        >  2:621744/4‑83 (MQ=255)
                                                                                                                    tctgAAGCACTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATc                                      >  1:312671/4‑117 (MQ=255)
                                                                                                                    tctgAAGCACTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATc                                      <  2:312671/114‑1 (MQ=255)
                                                                                                                      ggAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGt                                                                >  2:740871/1‑89 (MQ=255)
                                                                                                                      ggAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGt                                                                <  1:740871/89‑1 (MQ=255)
                                                                                                                       gAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTgg                                                              >  2:358298/1‑90 (MQ=255)
                                                                                                                       gAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTgg                                                              <  1:358298/90‑1 (MQ=255)
                                                                                                                       gAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTCGCGTAAGCCAGAAAGGAATCGTGGACCTTTAt  <  2:422547/150‑1 (MQ=255)
                                                                                                                            aGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCg                                     >  2:431917/1‑110 (MQ=255)
                                                                                                                            aGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCg                                     <  1:431917/110‑1 (MQ=255)
                                                                                                                             gTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTCGCGTAAGCCAGAAAGGAATCGTa           <  2:261562/135‑1 (MQ=255)
                                                                                                                             gTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTCGCGTAAGCCAGAAAGGAATCGTa           >  1:261562/1‑135 (MQ=255)
                                                                                                                             |                                                                                                                                               
TAGCTAATGGCGTCTCCAGTGATATTCCGTGGACGCAGAAAGATACCTCTTACACCAAACGGGTTGTATGCGGGCATATCGGCGGAATGTCGGAGCACAAGGGGTATTATCTGCTTAAGGAAGCAGTGCTTCAGAGTCAGCCTGAGAACCTTGAGTTTTTGATTGTTGACCATTCCAAGGAAGAGGGGTATGAACATCATTCTCAGTGGGGGAAAGTGCCAGTGACCTTTATCGGTCGCGTAAGCCAGAAAGGAATCGTAGACCTTTAT  >  NC_002947/2002945‑2003213

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: