Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,457,153 G→A 11.2% Y111Y (TAC→TAT)  PP_1275 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,457,1530GA11.2% 132.6 / 8.6 54Y111Y (TAC→TAT) PP_1275hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base A (3/3);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

AGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGT  >  NC_002947/1457011‑1457297
                                                                                                                                              |                                                                                                                                                
gtgaAACTCATCAATTGGCAAACGTACCTTGCAGGTGTTCGGCTCGGAATTCCAGGGTCTGCGGCCGATCGCCGGAACGTAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATcgcg                                                                                                                                           <  1:239811/148‑1 (MQ=255)
  gaAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCAGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCggg                                                                                                                                         <  2:855903/150‑1 (MQ=255)
      cacaTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGATACCGCgggccag                                                                                                                                     <  2:860350/150‑5 (MQ=255)
       acacCAATTGGCAAACGTACCTTGCAGGTGCGCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAAGGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGtt                                                                                                                                    <  1:998378/150‑1 (MQ=255)
       acaTCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCTGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCACCGGCTGCAATTCACCTGGCCCTCGGTGTCGCGGGGCGtt                                                                                                                                    >  2:40156/1‑150 (MQ=255)
          tCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTa                                                                                                                                 <  2:1044655/150‑1 (MQ=255)
            aaTTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATATCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCgggg                                                                                                                                        >  2:818522/1‑141 (MQ=255)
            aaTTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATATCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCgggg                                                                                                                                        <  1:818522/141‑1 (MQ=255)
                ggCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGTCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGtt                                                                                                                           >  1:726678/1‑150 (MQ=255)
                   aaaCGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTtctc                                                                                                                        >  1:379336/1‑150 (MQ=255)
                   aaaCGTACATTGCAGGTGATCGGCTCGGAATTCAAGGGGCTGCGGCCGAGAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTtctc                                                                                                                        >  2:1030875/1‑150 (MQ=255)
                    aaCGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCg                                                                                                                       <  2:440433/150‑1 (MQ=255)
                    aaCGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCg                                                                                                                       <  1:841320/150‑1 (MQ=255)
                           cTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGt                                                                                                                >  1:721298/1‑150 (MQ=255)
                           cTTGCAGGTGATCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGt                                                                                                                >  2:376038/1‑150 (MQ=255)
                            ttGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGtg                                                                                                               >  1:612934/1‑149 (MQ=255)
                              gCAGGTGCTCGGCTCGGAATTCCAGGCTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCCCCTGGCCCGCCGTACGGCGGGGCGTTGTAGTGGTGCTCTGCCAGGGTcc                                                                                                             >  2:560346/1‑149 (MQ=255)
                                    gctcggctcgGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTAt                                                                                                       <  1:121924/150‑1 (MQ=255)
                                     ctcggctcggAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATc                                                                                                      <  2:898426/150‑1 (MQ=255)
                                              gAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAg                                                                                                                  >  1:391195/1‑129 (MQ=255)
                                              gAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAg                                                                                                                  <  2:391195/129‑1 (MQ=255)
                                              gAATTCCAGGGTCTGCGGCCGATACCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGc                                                                                             <  1:72094/150‑1 (MQ=255)
                                                 ttCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTa                                                                                          <  2:859293/150‑1 (MQ=255)
                                                 ttCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGCCCCGCGATACCGCgggccag                                                                                                                                     <  1:806889/107‑5 (MQ=255)
                                                 ttCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGCCCCGCGATACCGCgggccag                                                                                                                                     >  2:806889/1‑103 (MQ=255)
                                                            gCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCAcgc                                                                               <  2:379336/150‑1 (MQ=255)
                                                                   aTAGCCGGAACGCTCCCGGGGCAATGGCAGGCAGTCTTGCCACTCGGCCGCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGGTGTAGGTGTTCTCGGCCAGGGTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAgcg                                                                        >  1:587615/1‑150 (MQ=255)
                                                                     aGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCa                                                                      <  1:707945/150‑1 (MQ=255)
                                                                                      ggcAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCAcg                                                     >  1:1111807/1‑150 (MQ=255)
                                                                                      ggcAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCAcg                                                     >  1:222139/1‑150 (MQ=255)
                                                                                        cAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGCACCGCGTTACCAAGCCGAGTTTTGGGTTTTCTCGGCCAGTTGCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCAcgtc                                                   <  1:523568/150‑1 (MQ=255)
                                                                                                       ttGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCAt                                    >  1:242799/1‑150 (MQ=255)
                                                                                                         gCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGc                                  <  2:721298/150‑1 (MQ=255)
                                                                                                            aCTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTa                               >  2:604186/1‑150 (MQ=255)
                                                                                                             cTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGAATGCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTAc                              <  1:1030875/150‑1 (MQ=255)
                                                                                                              tCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTAcc                             <  1:160830/150‑1 (MQ=255)
                                                                                                                   ccccGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCaa                        <  2:726678/150‑1 (MQ=255)
                                                                                                                     ccGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAgg                      >  1:265544/1‑150 (MQ=255)
                                                                                                                      cGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCt                            >  1:955369/1‑143 (MQ=255)
                                                                                                                      cGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCt                            <  2:955369/143‑1 (MQ=255)
                                                                                                                        gCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCAc                   <  1:560346/150‑1 (MQ=255)
                                                                                                                          tGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTCCCTGCAAGGCACgg                 >  2:177117/1‑150 (MQ=255)
                                                                                                                          tGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACgg                 >  2:1085611/1‑150 (MQ=255)
                                                                                                                           gCAATTCACCTGGCCCGCGGTATCTCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCa                                                                           >  1:145706/1‑91 (MQ=255)
                                                                                                                           gCAATGCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCa                                                                           <  2:145706/91‑1 (MQ=255)
                                                                                                                             aaTTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTcg              >  1:303515/1‑150 (MQ=255)
                                                                                                                               ttCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTcggc            >  2:732663/1‑150 (MQ=255)
                                                                                                                                  aCCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCggg         >  2:912891/1‑150 (MQ=255)
                                                                                                                                   ccTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGa        <  1:376038/150‑1 (MQ=255)
                                                                                                                                   aacgccccGCGATACCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCAt                                    <  1:158959/117‑1 (MQ=255)
                                                                                                                                   aacgccccGCGATACCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCAt                                    >  2:158959/6‑122 (MQ=255)
                                                                                                                                    acgccccGCGATACCGCGGGGCGTTGTAACACTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGAc       >  1:29682/5‑150 (MQ=255)
                                                                                                                                     tGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGATAAGGTTCCAGCTACCCCCCCGGGCGTACGCCGCCACGCGCCAGCGCAGGTCGGTCAGCGGCACGTCGTTGCCGTTCCTCATGCGTACCTGCAAGGCACGGTCGGCGGGACa      <  2:242799/150‑1 (MQ=255)
                                                                                                                                         ccGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAgtgt  >  2:981834/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                
AGGAAACACATCAATTGGCAAACGTACCTTGCAGGTGCTCGGCTCGGAATTCCAGGGTCTGCGGCCGATAGCCGGAACGCAGCGGCGGCAATGGCAGGCAGTCTTGCCACTCGGCCCCCGGCTGCAATTCACCTGGCCCGCGGTATCGCGGGGCGTTGTAGGTGTTCTCGGCCAGGTTCACGGTATCGCCCGGTGCGTACGCCGCCACGCGCCAGCGCAGTTCGGTCAGCGGCACGTCGTTGCCGTTCTTCATGCGTACCTGCAAGGCACGGTCGGCGGGACAGTGT  >  NC_002947/1457011‑1457297

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