Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,054,193:1 +T 8.9% coding (767/1320 nt) PP_3573 → monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,054,1931.T 8.9% 209.1 / 14.8 56coding (767/1320 nt)PP_3573monooxygenase
Reads supporting (aligned to +/- strand):  ref base . (21/30);  new base T (3/2);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGG  >  NC_002947/4054048‑4054339
                                                                                                                                                   |                                                                                                                                                  
cGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCAt                                                                                                                                                 <  1:167568/149‑1 (MQ=255)
        cgACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTgg                                                                                                                                         <  2:88795/149‑1 (MQ=255)
         gACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGt                                                                                                                                        <  6:34841/149‑1 (MQ=255)
              ctgctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGc                                                                                                                                    <  1:255613/148‑1 (MQ=255)
               tgctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGAGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCt                                                                                                                                  <  6:34068/149‑1 (MQ=255)
                 ctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGa                                                                                                                                <  8:245648/149‑1 (MQ=255)
                 ctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGa                                                                                                                                <  6:110933/149‑1 (MQ=255)
                    gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCGACCAGGGTCATGCCCTGGTGGGcgaggg                                                                                                                               <  7:204177/147‑6 (MQ=255)
                        cTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTc                                                                                                                         <  3:128092/149‑1 (MQ=255)
                         tGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCTACCAGGGTCATGCCCTGGTggg                                                                                                                                     <  1:277775/136‑1 (MQ=255)
                         tGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCTACCAGGGTCATGCCCTGGTggg                                                                                                                                     >  2:277775/1‑136 (MQ=255)
                            gCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGAcc                                                                                                                              <  5:292271/140‑1 (MQ=255)
                            gCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGAcc                                                                                                                              >  6:292271/1‑140 (MQ=255)
                                gCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGaa                                                                                                                  >  6:193769/1‑148 (MQ=255)
                                            ggACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTcgcgc                                                                                                     >  5:180472/1‑149 (MQ=255)
                                              aCGCGGAAATTGCCAAGCCCGGTCGCGAGCATGGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCAGCCTGACCCAGTCGTTCGAAAACGGc                                                                                                            >  1:43068/1‑141 (MQ=255)
                                              aCGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCgcgcgc                                                                                                   <  4:303342/149‑1 (MQ=255)
                                              aCGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCgcgcgc                                                                                                   <  3:349826/149‑1 (MQ=255)
                                              aCGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCAGCCTGACCCAGTCGTTCGAAAACGGc                                                                                                            <  2:43068/140‑1 (MQ=255)
                                                gcgGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCtt                                                                                                 >  6:316468/1‑149 (MQ=255)
                                                gcgGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCtt                                                                                                 >  1:139590/1‑149 (MQ=255)
                                                 cgGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTc                                                                                                <  6:41521/149‑1 (MQ=255)
                                                     aaTTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGTCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAgg                                                                                            <  5:193769/149‑1 (MQ=255)
                                                     aaTTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGAcc                                                                                                                              >  6:191384/1‑115 (MQ=255)
                                                     aaTTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGAcc                                                                                                                              <  5:191384/115‑1 (MQ=255)
                                                          ccAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAAc                                                                                       <  5:256665/149‑1 (MQ=255)
                                                           cAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAAcc                                                                                      >  7:43861/1‑149 (MQ=255)
                                                           cAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAAcc                                                                                      >  8:18204/1‑149 (MQ=255)
                                                               cccGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGt                                                                                  <  1:229599/149‑1 (MQ=255)
                                                                ccGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTc                                                                                 <  2:139590/149‑1 (MQ=255)
                                                                    tCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCgaga                                                                             <  1:206022/149‑1 (MQ=255)
                                                                        gAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACAt                                                                         >  3:89326/1‑149 (MQ=255)
                                                                               gtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATcg                                                                   >  1:153241/1‑149 (MQ=255)
                                                                                tgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATcgc                                                                  >  4:26284/1‑149 (MQ=255)
                                                                                  aCCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcgg                                                                >  5:180630/1‑149 (MQ=255)
                                                                                    cATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcggc                                                               >  5:359233/1‑148 (MQ=255)
                                                                                     aTCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTGCCGTGCCCTCGCCCACCAGGGCATGACCCTCGACGGCCTGACCCAGTCGTACGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcga                                                             >  8:253783/1‑149 (MQ=255)
                                                                                                           gtggCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCctgct                                        <  1:75806/148‑1 (MQ=255)
                                                                                                            tggCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTgg                                      <  6:180630/149‑1 (MQ=255)
                                                                                                              gCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGAt                                    <  5:222476/149‑1 (MQ=255)
                                                                                                                     cggggACTTCCGTGCCCGCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACTGCGTCGCGCGCTTCCATGACAACCTGGTCGAGAACATCAATCgcgg                                                                <  4:37741/111‑1 (MQ=255)
                                                                                                                     cGTGGACTTACGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGCCCCAGTCGTGCGAAAACGGAGTCGCGCGCTTCCAGGACAACCGGGTCGAGAACATCAATCgcgg                                                                >  3:37741/1‑114 (MQ=255)
                                                                                                                      gTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACg                                                                                                              <  4:235884/67‑1 (MQ=255)
                                                                                                                      gTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACg                                                                                                              >  3:235884/1‑67 (MQ=255)
                                                                                                                       tGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGAt                           <  3:26284/149‑1 (MQ=255)
                                                                                                                          aCTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGcc                        <  7:18204/149‑1 (MQ=255)
                                                                                                                           cTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCt                       <  6:359233/149‑1 (MQ=255)
                                                                                                                                gTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTgg                                               <  3:28968/120‑1 (MQ=255)
                                                                                                                                gTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGCTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATc                  <  3:230528/149‑1 (MQ=255)
                                                                                                                                gTGCCCTCGCCCACCAGGGCATGACACTGGTCGGCCTGACCCAGTCGATCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTgg                                               >  4:28968/1‑120 (MQ=255)
                                                                                                                                    ccTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAAc              >  3:241145/1‑149 (MQ=255)
                                                                                                                                    ccTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAAc              >  7:284044/1‑149 (MQ=255)
                                                                                                                                          cgaccaGGGTCATGCCCTGGTGGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACg        >  2:128337/3‑149 (MQ=255)
                                                                                                                                          cgaccaGGGTCATGCCCTGGTGGGCCTGACCCAGTCGTTCGAAAACGGCGACGCGCGCATCCAGGACAACCTGGTCGAGAACAACAATCGAGGCGACGACCACTACCTGGACCTGCAGCATGCGGCCGATGCCTACATCGAACGCaa          >  4:112288/3‑147 (MQ=255)
                                                                                                                                          ccaccaGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcgacga                                                          >  6:156754/1‑99 (MQ=255)
                                                                                                                                          ccaccaGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcgacga                                                          <  5:156754/99‑1 (MQ=255)
                                                                                                                                                gggCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTgg  <  4:241145/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGG  >  NC_002947/4054048‑4054339

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: