Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,107,580 T→C 7.5% Y41H (TAC→CAC)  PP_3615 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,107,5800TC 7.5% 148.2 / 3.1 53Y41H (TAC→CAC) PP_3615hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/21);  new base C (2/2);  total (30/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 = 7.35e-01

GCGCCATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCAAGGCCGTCCCCACTACCGGGGAAGGCGGCCTGGCCTG  >  NC_002947/4107439‑4107728
                                                                                                                                             |                                                                                                                                                    
gcgcCATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTggg                                                                                                                                                <  8:232330/148‑1 (MQ=255)
    cATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGGAGACCAAACGCGTACTGGGCCCTg                                                                                                                                           <  7:50669/149‑1 (MQ=255)
    cATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTAc                                                                                                                                                    <  7:335115/140‑1 (MQ=255)
    cATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTAc                                                                                                                                                    >  8:335115/1‑140 (MQ=255)
                    ggTGCGCTTAGCAAAGGAGCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTggg                                                                                                                                                <  7:160108/128‑1 (MQ=255)
                    ggTGCGCTTAGCAAAGGAGCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTggg                                                                                                                                                >  8:160108/1‑128 (MQ=255)
                        cgcTTAGCAAAGGAGCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTgg                                                                                                                                          <  6:10166/130‑1 (MQ=255)
                        cgcTTAGCAAAGGAGCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTgg                                                                                                                                          >  5:10166/1‑130 (MQ=255)
                             aGCAAAGGAGCCCCGTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTg                                                                                                                   <  8:351900/148‑1 (MQ=255)
                                aaaGGAGCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAAAGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGctg                                                                                                               >  1:77706/1‑148 (MQ=255)
                                 aaGGAGCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGctct                                                                                                              <  7:86138/149‑1 (MQ=255)
                                    gAGCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGc                                                                                                           >  1:207118/1‑149 (MQ=255)
                                      gCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCg                                                                                                         >  1:100537/1‑149 (MQ=255)
                                         ccTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATATGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGgct                                                                                                      >  7:82693/1‑149 (MQ=255)
                                          cTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGgctg                                                                                                     >  6:64034/1‑149 (MQ=255)
                                          cTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGgctg                                                                                                     >  3:72734/1‑149 (MQ=255)
                                                 ccccGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGc                                                                                              <  4:121002/149‑1 (MQ=255)
                                                      aTAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATATGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGc                                                                                         <  8:82693/149‑1 (MQ=255)
                                                             aCCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCaa                                                                                   <  6:156639/148‑1 (MQ=255)
                                                             aCCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAg                                                                                  >  5:293271/1‑149 (MQ=255)
                                                             aCCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATAAATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCGGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAg                                                                                  >  2:297550/1‑149 (MQ=255)
                                                                       tGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTacg                                                                        >  3:209080/1‑149 (MQ=255)
                                                                           aTGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTacgccac                                                                    >  5:202839/1‑149 (MQ=255)
                                                                            tGAAATATCCTGCAAGGCATTTGCTGGCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCACCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTacgccacg                                                                   >  1:198704/1‑149 (MQ=255)
                                                                             gAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTacgccacgc                                                                  >  3:211146/1‑149 (MQ=255)
                                                                              aaaTATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTacgccacgcc                                                                 >  3:278024/1‑149 (MQ=255)
                                                                               aaTATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCt                                                                <  4:278024/149‑1 (MQ=255)
                                                                                 tatCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGCACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGc                                                                                              <  6:163871/117‑1 (MQ=255)
                                                                                 tatCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGCACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGc                                                                                              <  6:158365/117‑1 (MQ=255)
                                                                                 tatCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGCACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGc                                                                                              >  5:158365/1‑117 (MQ=255)
                                                                                 tatCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGCACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGc                                                                                              >  5:163871/1‑117 (MQ=255)
                                                                                   tCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCa                                                             <  3:307414/148‑1 (MQ=255)
                                                                                     cTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCACCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGt                                                           <  7:175595/148‑1 (MQ=255)
                                                                                      tGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTgg                                                         >  4:202442/1‑149 (MQ=255)
                                                                                      tGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCAAGCCTGCAGTgg                                                         >  3:207997/1‑149 (MQ=255)
                                                                                       gCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTggg                                                        >  1:33240/1‑149 (MQ=255)
                                                                                         aaGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCt                                                      >  1:135343/1‑149 (MQ=255)
                                                                                           ggCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTc                                                     >  6:135154/1‑148 (MQ=255)
                                                                                              aTTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAAc                                                 <  4:72734/149‑1 (MQ=255)
                                                                                                   ctgcCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCtt                                            >  1:67759/1‑149 (MQ=255)
                                                                                                   ctgcCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCt                                             <  3:202442/148‑1 (MQ=255)
                                                                                                      ccAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCttt                                           <  4:207997/147‑1 (MQ=255)
                                                                                                           cgcAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCAAgg                                    <  8:262288/149‑1 (MQ=255)
                                                                                                                  ggggAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCAAGGCCGTccc                             >  5:348284/1‑149 (MQ=255)
                                                                                                                     gAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCAAGGCCGTCCCCAc                          >  5:60147/1‑149 (MQ=255)
                                                                                                                      aaTACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCAAGGCCGTCCCCAc                          >  7:170922/1‑148 (MQ=255)
                                                                                                                        tACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCaa                                      >  3:219254/1‑134 (MQ=255)
                                                                                                                        tACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCaa                                      <  4:219254/134‑1 (MQ=255)
                                                                                                                           aTGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCAAGGCCGTCCCCACTACCgg                    <  4:211146/149‑1 (MQ=255)
                                                                                                                                cAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCAAGGCCGTCCCCACTACCGGGGAAg               >  6:286256/1‑149 (MQ=255)
                                                                                                                                 aGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCAAGGCCGTCCCCACTACCGGGGAAgg              <  3:138147/149‑1 (MQ=255)
                                                                                                                                      aaaCGCGTACTGGGCACTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCAAGGCCGTCCCCACTACCGGGGAAGGCggcc         <  2:100537/149‑1 (MQ=255)
                                                                                                                                       aaCGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCAAGGCCGTCCCCACTACCGGGGAAGGCggcct        >  6:181295/1‑149 (MQ=255)
                                                                                                                                             tACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCAAGGCCGTCCCCACTACCGGGGAAGGCGgcctggcctg  <  2:135343/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
GCGCCATCAGGTAAGCGCAGGGTGCGCTTAGCAAAGGAGCCCCTTGACGCCCCGATAGTTCACCGTTTCGCTGTAATGAAATATCCTGCAAGGCATTTGCTGCCAAGCGCAGAAGGGGAATACATGAGCAGACCAAACGCGTACTGGGCCCTGGCCACCCTGCTCAGCCTCAGCCTGGCTCTCGCCGGCTGCTCCAGCAACGCCAGCAAGGCACGCTACGCCACGCCTGCAGTGGGCTCCAACTGCTTTGCCAAGGCCGTCCCCACTACCGGGGAAGGCGGCCTGGCCTG  >  NC_002947/4107439‑4107728

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