Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,554,354 G→A 8.7% A181V (GCG→GTG)  PP_4886 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,554,3540GA 8.7% 126.2 / 2.4 46A181V (GCG→GTG) PP_4886hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/20);  new base A (2/2);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

AACAGTGTTGACGACATGCGCATTCAGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTCTTGATGCGCCCGTAGTTCAGCTGCACGACTTCCACTGGCACCCCAGCGCTGGCGCTTTG  >  NC_002947/5554212‑5554502
                                                                                                                                              |                                                                                                                                                    
aacaGTGTTGACGACATGCGCATTCAGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATAct                                                                                                                                                >  1:201062/1‑149 (MQ=255)
aacaGTGTTGACGACATGCGCATTCAGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATAct                                                                                                                                                >  7:248948/1‑149 (MQ=255)
  caGTGTTGACGACATGCGCATTCAGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAAAGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATActtc                                                                                                                                              >  2:29546/1‑149 (MQ=255)
  aaGTGTTGACGACATGCTCATTCAGAACATCAAACTTTGCCCTCTGATTCGCTTTCATGGTCGTATCGGCCTTCAGCGAATCATACGACTGCGCATGCTGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATActtc                                                                                                                                              <  3:233535/148‑1 (MQ=18)
    gtgtTGACGACATGCGCATTCAGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATActtctt                                                                                                                                            <  6:267452/149‑1 (MQ=255)
      gtTGACGACATGCGCATTCAGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCttg                                                                                                                                           <  6:203345/148‑1 (MQ=255)
         gacgacATGCGCATTCAGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCttgttg                                                                                                                                        <  6:99184/148‑1 (MQ=255)
             acATGCGCATTCAGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAAAGAAGCTGCGCGCCTCGGACATGGAAGTACCTCACGCATACTTCTTGTTGTTGAt                                                                                                                                   <  1:29546/149‑1 (MQ=255)
                   gcATTCAGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGc                                                                                                                                                     <  6:308589/125‑1 (MQ=255)
                   gcATTCAGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGc                                                                                                                                                     >  5:308589/1‑125 (MQ=255)
                     aTTCAGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATActtc                                                                                                                                              <  1:60236/130‑1 (MQ=255)
                     aTTCAGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATActtc                                                                                                                                              >  2:60236/1‑130 (MQ=255)
                         aGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATcc                                                                                                                       >  3:239435/1‑149 (MQ=255)
                              aTCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCg                                                                                                                  >  1:300280/1‑149 (MQ=255)
                                 aaaCTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTc                                                                                                                             <  8:281918/135‑1 (MQ=255)
                                 aaaCTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTc                                                                                                                             >  7:281918/1‑135 (MQ=255)
                                  aaCTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCttg                                                                                                                                           >  2:280836/1‑120 (MQ=255)
                                  aaCTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCttg                                                                                                                                           <  1:280836/120‑1 (MQ=255)
                                   aCTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACg                                                                                                             <  2:201062/149‑1 (MQ=255)
                                          ccGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTgccgc                                                                                                      <  4:239435/149‑1 (MQ=255)
                                                       ttCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTc                                                                                         <  2:300280/149‑1 (MQ=255)
                                                             ggCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGc                                                                                    <  7:54728/148‑1 (MQ=255)
                                                              gCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGt                                                                                  >  7:22929/1‑149 (MQ=255)
                                                                 tATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTc                                                                               <  8:22929/149‑1 (MQ=255)
                                                                  aTCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCt                                                                              >  1:315109/1‑149 (MQ=255)
                                                                        cTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAg                                                                                                   >  2:40754/1‑122 (MQ=255)
                                                                        cTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAg                                                                                                   <  1:40754/122‑1 (MQ=255)
                                                                              cGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCATTTCTGCCGCgtg                                                                      >  3:29129/1‑145 (MQ=255)
                                                                              cGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACACATACTTCTTGTTGTTGATGTGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCcgc                                                                         <  5:266449/142‑1 (MQ=255)
                                                                              cGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACACATACTTCTTGTTGTTGATGTGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCcgc                                                                         >  6:271331/1‑142 (MQ=255)
                                                                              cGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACACATACTTCTTGTTGTTGATGTGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCcgc                                                                         >  6:266449/1‑142 (MQ=255)
                                                                              cGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACACATACTTCTTGTTGTTGATGTGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCcgc                                                                         <  5:271331/142‑1 (MQ=255)
                                                                                  tCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGa                                                                  <  4:22566/145‑1 (MQ=255)
                                                                                  tCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGa                                                                  >  3:22566/1‑145 (MQ=255)
                                                                                    aTACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTc                                                             <  6:295155/148‑1 (MQ=255)
                                                                                    aTACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTc                                                             >  5:295155/1‑148 (MQ=255)
                                                                                    aTACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTCt                                                            >  4:301991/1‑149 (MQ=255)
                                                                                    aTACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTCt                                                            >  7:267879/1‑149 (MQ=255)
                                                                                       cGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCATTTCTGCCGCGTGTAGATTGTCTTg                                                          <  4:29129/148‑1 (MQ=255)
                                                                                                aTGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTCTTGATGCGCCCg                                                 >  3:45808/1‑148 (MQ=255)
                                                                                                 tGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTCt                                                            <  6:289084/136‑1 (MQ=255)
                                                                                                 tGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTCt                                                            >  5:289084/1‑136 (MQ=255)
                                                                                                          aTGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTCTTGATGCGCCCGTAGTTCAGCTg                                      >  3:139905/1‑149 (MQ=255)
                                                                                                             aaGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTCTTGATGCGCCCGTAGTTCAGCTGCa                                    <  8:267879/148‑1 (MQ=255)
                                                                                                                             aCATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTCTTGATGCGCCCGTAGTTCAGCTGCACGACTTCCACTGGCAcc                   >  2:318573/1‑149 (MQ=255)
                                                                                                                                          tCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCATGTAGATTGTCTTGATGCGCCCGTAGTTCAGCTGCACGACTTCCACTGGCACCCCAGCGCTGgcgc      >  1:188719/1‑149 (MQ=255)
                                                                                                                                              gCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTCTTGATGCGCCCGTAGTTCAGCTGCACGACTTCCACTGGCACCCCAGCGCTGGCGCTTTg  >  5:304360/1‑149 (MQ=255)
                                                                                                                                              gCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTCTTGATGCGCCCGTAGTTCAGCTGCACGACTTCCACTGGCACCCCAGCGCTGGCGCTTTg  <  3:45354/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
AACAGTGTTGACGACATGCGCATTCAGCACATCAAACTTGGCCCGCTGATTCGCGTTCATGGGCGTATCGGCCTTCAGCGACTCATACGACTGCGCATGCGGATTAATGAAGCTGCGCGCCTCGGACATGGTAGTACCTCACGCATACTTCTTGTTGTTGATGCGGTCCCATCCGCCCGTGACGTTGCCGCCAGCGCTGCCGTCGAGGCGTTTCTGCCGCGTGTAGATTGTCTTGATGCGCCCGTAGTTCAGCTGCACGACTTCCACTGGCACCCCAGCGCTGGCGCTTTG  >  NC_002947/5554212‑5554502

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