Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,741,382 G→A 7.4% N191N (AAC→AAT PP_1538 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,741,3820GA 7.4% 151.8 / 3.1 54N191N (AAC→AATPP_1538hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base A (2/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.53e-01

CAGACCCAGGGGTTGGCGTTCCAGGAGTCGGCCCCGTTGATGGATACCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGAC  >  NC_002947/1741233‑1741523
                                                                                                                                                     |                                                                                                                                             
cagacccagGGGTTGGCGTTCCAGGAGTCGGCCCCGTTGATGGATACCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCgg                                                                                                                                               <  2:1245080/150‑1 (MQ=255)
 agacccagGGGTTGGCGTTCCAGGAGTCGGCCCCGTTGATGGATACCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGt                                                                                                                                              <  2:961778/150‑1 (MQ=255)
     ccagGGGTTGGCGTTCCAGGAGTCGGCCCCGTTGATGGATCAACAGGGATAACGACAAGCTTCGAATGGAGCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTAt                                                                                                                                          >  1:989288/1‑150 (MQ=255)
      cagGGGTTGGCGTTCCAGGAGTCGGCCCCGTTGATGGATACCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTTAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATg                                                                                                                                         >  2:283890/1‑150 (MQ=255)
           gTTGGCGTTCCAGGAGTCGGCCCCGTTGATGGATACCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCACGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTc                                                                                                                                    <  1:414069/150‑1 (MQ=255)
              ggCGTTCCAGGAGTCGGCCCCGTTGATGGATACCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCggg                                                                                                                                 <  1:596492/150‑1 (MQ=255)
               gCGTTCCAGGAGTCGGCCCCGTTGATGGATACCCAGAGGTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTAAACCATCGCCGTGGTGAATGCGGTTTATGCCCTCgggg                                                                                                                                >  2:743279/1‑150 (MQ=255)
                   tCCAGGAGTCGGCCCCGTTGATGGATACCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGt                                                                                                                            <  2:645908/150‑1 (MQ=255)
                        gAGTCGGCCACGTTGATGGATACCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAGTGGGTGGTAGTGGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCgg                                                                                                                       >  1:30495/1‑150 (MQ=255)
                            cGGCCCCGCTGATCGCTACCCAGCGCTCGCGCCACGCGTGGAAGGGTCCGGCGCAGTTTACCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCttt                                                                                                                   <  2:67157/150‑1 (MQ=255)
                                       aTGGATACCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAAt                                                                                                             >  2:296164/1‑145 (MQ=255)
                                       aTGGATACCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAAt                                                                                                             <  1:296164/145‑1 (MQ=255)
                                         ggATACCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGtctc                                                                                                      >  1:419289/1‑150 (MQ=255)
                                            tACCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTctcgcc                                                                                                   >  1:976651/1‑150 (MQ=255)
                                             aCCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTctcgccc                                                                                                  <  1:1081037/150‑1 (MQ=255)
                                              cgcggggCTCGCGCCAAGCGTCGAAGAGATCGTCCCAGTGTCCCGGCCCGGGATCGCTCTTGAGTGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAGTCCGTctcgccct                                                                                                 <  2:30495/144‑1 (MQ=255)
                                               ccAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGGATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCgg                                                                                                                                               >  2:115947/1‑103 (MQ=255)
                                               ccAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGGATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCgg                                                                                                                                               <  1:115947/103‑1 (MQ=255)
                                                   agCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCa                                                                                            <  1:394827/150‑1 (MQ=255)
                                                        gcgcCAAGCGTCGATCCTAGCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTg                                                                                       >  1:1008665/1‑150 (MQ=255)
                                                         cgcCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGc                                                                                      <  2:153084/150‑1 (MQ=255)
                                                             aaGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGgc                                                                                  <  2:419289/150‑1 (MQ=255)
                                                               gCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTc                                                                                                                  <  1:676925/116‑1 (MQ=255)
                                                               gCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTc                                                                                                                  >  2:676925/1‑116 (MQ=255)
                                                                    aaatggATCGGCCCTGTTTCCCGTTCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCTATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGac                                                                           <  1:218910/146‑1 (MQ=255)
                                                                     aaGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGaca                                                                          >  2:842762/1‑150 (MQ=255)
                                                                         gATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGAcacgca                                                                      >  2:48891/1‑150 (MQ=255)
                                                                                 cAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCCTGCCATCCTGCAGGCGCTCTACACGCACGTCGGTg                                                              <  1:743279/150‑1 (MQ=255)
                                                                                         cGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCt                                                                                                                                     >  2:1266583/1‑71 (MQ=255)
                                                                                         cGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCt                                                                                                                                     <  1:1266583/71‑1 (MQ=255)
                                                                                         cGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCcag                                                      >  1:177908/1‑150 (MQ=255)
                                                                                         cGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCcag                                                      >  2:630849/1‑150 (MQ=255)
                                                                                           gCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTc                                                                                                                  >  1:799425/1‑88 (MQ=255)
                                                                                           gCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTc                                                                                                                  <  2:799425/88‑1 (MQ=255)
                                                                                              cAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGATTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGa                                                 >  2:801501/1‑150 (MQ=255)
                                                                                              cAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGATTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGa                                                 >  2:801526/1‑150 (MQ=255)
                                                                                                       tctTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCggg                                                                                                                                 <  2:1221915/61‑1 (MQ=255)
                                                                                                       tctTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCggg                                                                                                                                 >  1:1221915/1‑61 (MQ=255)
                                                                                                       tctTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCCGGATTcggccgg                                        >  2:518596/1‑150 (MQ=255)
                                                                                                         ttGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCAAAGAAGATTAggctggcc                                      >  2:913461/1‑150 (MQ=255)
                                                                                                              tGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTc                                                                                                                  >  1:515636/1‑69 (MQ=255)
                                                                                                              tGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTc                                                                                                                  <  2:515636/69‑1 (MQ=255)
                                                                                                                    gtagtaGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCAt                           <  2:989288/150‑1 (MQ=255)
                                                                                                                    gtagtaGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCAt                           <  1:1221073/150‑1 (MQ=255)
                                                                                                                      agtagtGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCAtgt                         <  1:48891/150‑1 (MQ=255)
                                                                                                                          gtgtTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCgg                                <  2:760872/139‑1 (MQ=255)
                                                                                                                          gtgtTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCgg                                >  1:760872/1‑139 (MQ=255)
                                                                                                                             ttagCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATc                                                                                                            >  2:601727/4‑60 (MQ=255)
                                                                                                                             ttCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATc                                                                                                            <  1:601727/60‑1 (MQ=255)
                                                                                                                             ttCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCt                  >  2:1001778/1‑150 (MQ=255)
                                                                                                                               cGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTgg                <  1:510031/150‑1 (MQ=255)
                                                                                                                                 ccATCGCCGTGGTGAATGCGATTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGc              <  1:801526/150‑1 (MQ=255)
                                                                                                                                 ccATCGCCGTGGTGAATGCGATTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGc              <  1:801501/150‑1 (MQ=255)
                                                                                                                                    tCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTc           <  2:976651/150‑1 (MQ=255)
                                                                                                                                         gtggtgAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGc      >  2:263883/1‑150 (MQ=255)
                                                                                                                                             tgAATGCGGTTTATTCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGAc  >  1:651732/1‑150 (MQ=255)
                                                                                                                                             tgAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGAc  <  1:263883/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
CAGACCCAGGGGTTGGCGTTCCAGGAGTCGGCCCCGTTGATGGATACCCAGAGCTCGCGCCAAGCGTCGAAGGGATCGGCCCAGTTTCCCGGCCCAGGATCGCTCTTGAATGGGTGGTAGTAGTGTTCGCCATCGCCGTGGTGAATGCGGTTTATGCCCTCGGGGACGTAGCGGCTTTCGAAATCCGTCTCGCCCTCGCCATCCTGCAGGCGCTCGACACGCACGTCGGTGACCTCCAGCAGGATTCGGCCGGCCACGCGGGGCATGTGGATGCTGGGCTTCCATGCTGAC  >  NC_002947/1741233‑1741523

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: