Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,680,805 G→A 9.3% Y150Y (TAC→TAT)  PP_2350 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,680,8050GA 9.3% 121.2 / 3.2 45Y150Y (TAC→TAT) PP_2350LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (22/18);  new base A (2/2);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

GCGGCATCATGCTGGCGGCGGAAGGCCTCGCCCCAGTCGATGTAGATGTACGGCTCGGGCCGCTCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAGGGCACTTACCCAGTCGAGCATCATCGGG  >  NC_002947/2680680‑2680949
                                                                                                                             |                                                                                                                                                
gcggcATCATGCTGGCGGCTGAAGGCCTCGCCCCAGTCGATGTAGATGTACGGCTCGGGCCGCTCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAgcgc                                                                                                                           <  1:247282/149‑1 (MQ=255)
    catcatGCTGGCGGCGGAAGGCCTCGCCCCAGTCGATGTAGATGTACGGCTCGGGCCGCTCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCgcggcg                                                                                                                       <  7:230079/149‑1 (MQ=255)
          gCTGGCGGCGGAAGGCCTCGCCCCAGTCGATGTAGATGTACGGCTCGGGCCGCTCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAg                                                                                                                  >  8:58070/1‑148 (MQ=255)
          gCTGGCGGCGGAAGGCCTCGCCCCAGTCGATGTAGATGTACGGCTCGGGCCGCTCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGAGCGGCGTCCAg                                                                                                                  >  7:122851/1‑148 (MQ=255)
            tGGCGGCGGAAGGCCTCGCCCCAGTCGATGTAGATGTACGGCTCGGGCCGCTCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTTGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCaa                                                                                                               >  8:70619/1‑149 (MQ=255)
            tGGCGGCGGAAGGCCTCGCCCCAGTCGATGTAGATGTACGGCTCGGGCCGCTCGACCCGGCGCACACGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGc                                                                                                                 >  5:17782/1‑147 (MQ=255)
                                      gatgtagatgtaCGGCTCGGGCCGCTCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTa                                                                                                                                                <  5:238818/90‑1 (MQ=255)
                                      gatgtagatgtaCGGCTCGGGCCGCTCGACCCGCCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTa                                                                                                                                                >  6:238818/1‑90 (MQ=255)
                                            gatgtaCGGCTCGGGCCGCTCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATtctc                                                                               >  7:210937/1‑149 (MQ=255)
                                               gtaCGGCTCGGGCCGCTCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGt                                                                            >  6:245244/1‑149 (MQ=255)
                                                gaCGGCTCGGGCCGCTCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGt                                                                            <  2:307779/147‑1 (MQ=255)
                                                  cGGCTCGGGCCGCTCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGc                                                                         >  3:26716/1‑149 (MQ=255)
                                                        gggCCGCTCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGTCCTTGCGCAGCAGCGATTCTCCGTCGCTGACCt                                                                   >  6:167175/1‑149 (MQ=255)
                                                               tCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTcc                                                                              >  6:54555/1‑131 (MQ=255)
                                                                cGACCCGGCGCACGCGAATCAGGTTCTCTTCCATCAACTGCTACACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCGCCGTCGCAGACCTCGCTGcc                                                            >  7:199367/1‑147 (MQ=255)
                                                                   cccGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGatg                                                        <  8:210937/149‑1 (MQ=255)
                                                                     cGGCGCACGCGAATCAGCTGCTCTTCCAGCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCCGCGCGGCGTCCAGCAAGCCCATAGAGCCCTTGTGCATCAGAGATTAGCCGTCGCTGACCTCGAAGCGGatgg                                                       >  7:276025/1‑148 (MQ=255)
                                                                           aCGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCt                                                >  8:212631/1‑149 (MQ=255)
                                                                            cgcgAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTg                                               <  7:212631/149‑1 (MQ=255)
                                                                             gcgAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTTGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTgg                                              <  7:70619/149‑1 (MQ=255)
                                                                             acgAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTgg                                              <  6:17782/148‑1 (MQ=255)
                                                                              cgAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTggg                                             <  5:245244/149‑1 (MQ=255)
                                                                                  tCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGt                                         <  6:316336/149‑1 (MQ=255)
                                                                                     gctgctCTGCCATCAAATGCGACAGCTGACGGCCTGGACCGTAGGTCGGCTGGTATCCCAGAGCGGCGTCCAGCAAGCCTATCTCGCCCTTGTGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCa                                      <  8:276025/149‑1 (MQ=255)
                                                                                     gCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCCGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGAc                                                                     >  1:124179/1‑118 (MQ=255)
                                                                                          tctTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGccc                                                                                                           >  7:179072/1‑75 (MQ=255)
                                                                                          tctTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGccc                                                                                                           <  8:179072/75‑1 (MQ=255)
                                                                                            ttCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTggg                                             <  8:240557/135‑1 (MQ=255)
                                                                                            ttCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTggg                                             >  7:240557/1‑135 (MQ=255)
                                                                                            ttCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTg                                 >  3:80670/1‑147 (MQ=255)
                                                                                            ttCCATCAACTGCTCCACCTGCAGGCCTGGACCATAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTg                                 >  4:300932/1‑147 (MQ=255)
                                                                                              ccATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAgg                             >  4:86068/1‑149 (MQ=255)
                                                                                              ccATCAACTGCTCCACCTGCAGGCCTGGACCATAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCa                               <  3:300932/147‑1 (MQ=255)
                                                                                                    aCTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCATCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAGGGCACt                        >  5:244149/1‑148 (MQ=255)
                                                                                                           cACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAgg                             <  8:100618/136‑1 (MQ=255)
                                                                                                           cACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAgg                             >  7:100618/1‑136 (MQ=255)
                                                                                                            ccccgCAGGGCTTGACCCGTGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCATCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAGGGCACTTACCCAGTc               <  8:199367/145‑1 (MQ=255)
                                                                                                              cTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCTTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAGGGCACTTACCCAGTCGa             <  6:244149/149‑1 (MQ=255)
                                                                                                                gCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAGGGCACTTAccc                   <  7:299105/141‑1 (MQ=255)
                                                                                                                gCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAGGGCACTTAccc                   >  8:299105/1‑141 (MQ=255)
                                                                                                                 cAGGCCTGGACCATAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAgg                             <  5:244423/130‑1 (MQ=255)
                                                                                                                 cAGGCCTGGACCATAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAgg                             >  6:244423/1‑130 (MQ=255)
                                                                                                                      gtGGACCGTAGGGCGGCTGGTAGACCACCGCGGCGTCCACCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCGGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAGGGCACTTAGCCAGTCGAGCatcatc     <  8:122851/148‑1 (MQ=255)
                                                                                                                         gACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAGGGCACTTACCCAGTCGAGCATCATCggg  >  7:27090/1‑149 (MQ=255)
                                                                                                                         gACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAGGGCACTTACCCAGTCGAGCATCATCg    <  7:58070/147‑1 (MQ=255)
                                                                                                                             |                                                                                                                                                
GCGGCATCATGCTGGCGGCGGAAGGCCTCGCCCCAGTCGATGTAGATGTACGGCTCGGGCCGCTCGACCCGGCGCACGCGAATCAGCTTCTCTTCCATCAACTGCTCCACCTGCAGGCCTGGACCGTAGGTCGGCTGGTAGACCAGCGCGGCGTCCAGCAAGCCCATCTCGACCTTGCGCAGCAGCGATTCTCCGTCGCTGACCTCGCTGCGGATGGCATGGCTGGGCAGTTCACGGTGCAGGGCACTTACCCAGTCGAGCATCATCGGG  >  NC_002947/2680680‑2680949

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