Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,614,088 T→C 7.6% intergenic (+42/‑40) PP_mr40 → / → pet TPP regulatory region/protein Pet18

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,614,0880TC 7.6% 147.9 / 2.2 53intergenic (+42/‑40)PP_mr40/petTPP regulatory region/protein Pet18
Reads supporting (aligned to +/- strand):  ref base T (26/22);  new base C (2/2);  total (29/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

ACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGATTGTGCGACGGCAGCCTGCCT  >  NC_002947/3613941‑3614232
                                                                                                                                                   |                                                                                                                                                
accaAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAAtt                                                                                                                                                 >  6:324591/1‑149 (MQ=255)
accaAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTCACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAAtt                                                                                                                                                 <  5:324591/149‑1 (MQ=255)
  caAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAAtt                                                                                                                                                 >  5:287347/1‑147 (MQ=255)
  caAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAAtt                                                                                                                                                 <  6:287347/147‑1 (MQ=255)
  caAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGc                                                                                                                                               >  4:261507/1‑149 (MQ=255)
   aaGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGc                                                                                                                                               >  8:189895/1‑148 (MQ=255)
       ggAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAgg                                                                                                                                          <  6:8055/149‑1 (MQ=255)
                 aaTGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTc                                                                                                                                <  6:46776/149‑1 (MQ=255)
                          agaAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCg                                                                                                                                     >  1:13490/1‑135 (MQ=255)
                          agaAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCg                                                                                                                                     <  2:13490/135‑1 (MQ=255)
                          agaAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCtcatc                                                                                                                       <  8:17726/149‑1 (MQ=255)
                           gaAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCtcatca                                                                                                                      <  3:261507/149‑1 (MQ=255)
                                               cgcgGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGaa                                                                                                  <  7:189895/149‑1 (MQ=255)
                                                    tGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGtt                                                                                              >  1:19202/1‑148 (MQ=255)
                                                      aCCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGttt                                                                                             >  1:171820/1‑147 (MQ=255)
                                                        ccTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGc                                                                                         >  6:137347/1‑149 (MQ=255)
                                                        ccTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGAAAGCAACATGACCGAACGTTTCAGc                                                                                         >  6:322332/1‑149 (MQ=255)
                                                          ttCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCatc                                                                                         >  5:178396/1‑145 (MQ=255)
                                                          ttCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGc                                                                                         >  3:288178/1‑147 (MQ=255)
                                                            cGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCa                                                                                       <  2:19202/147‑1 (MQ=255)
                                                             gAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTc                                                                                    >  1:299708/1‑149 (MQ=255)
                                                               aCCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTccc                                                                                  <  5:137347/149‑1 (MQ=255)
                                                                 cTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCaa                                                                                                                                            >  6:180986/1‑89 (MQ=255)
                                                                 cTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCaa                                                                                                                                            <  5:180986/89‑1 (MQ=255)
                                                                    aTCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGTGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTgcgc                                                                             >  1:19872/1‑149 (MQ=255)
                                                                        ggATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCa                                                                                           <  3:270724/131‑1 (MQ=255)
                                                                        ggATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCa                                                                                           >  4:270724/1‑131 (MQ=255)
                                                                                     gAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGcctgc                                                            <  6:194701/149‑4 (MQ=255)
                                                                                      aaGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGAcc                                                           <  4:288178/149‑1 (MQ=255)
                                                                                          gATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGt                                                       >  2:87169/1‑149 (MQ=255)
                                                                                            tGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTcg                                                     >  5:266078/1‑149 (MQ=255)
                                                                                                  ttGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGcc                                               >  5:15485/1‑149 (MQ=255)
                                                                                                    gTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGcc                                               <  2:299708/147‑1 (MQ=255)
                                                                                                     tAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGtt                                            <  6:178396/149‑1 (MQ=255)
                                                                                                      aaCCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGAc                                                            <  5:310463/132‑1 (MQ=255)
                                                                                                      aaCCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGAc                                                            >  6:310463/1‑132 (MQ=255)
                                                                                                         cTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGcc                                        >  4:289506/1‑149 (MQ=255)
                                                                                                         cTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAACTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGTGGCCGTTg                                           <  1:63274/146‑1 (MQ=255)
                                                                                                         cTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAACTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGTGGCCGTTg                                           >  2:63274/1‑146 (MQ=255)
                                                                                                         cTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAACTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGTGGCCGTTg                                           <  1:67788/146‑1 (MQ=255)
                                                                                                         cTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAACTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGTGGCCGTTg                                           >  2:67788/1‑146 (MQ=255)
                                                                                                            ccTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCAcc                                     >  7:190696/1‑149 (MQ=255)
                                                                                                              tttttGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCt                                                                                 <  5:148486/103‑1 (MQ=255)
                                                                                                              tttttGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTAGCTGTCACCCTCATCAAGGACAGAAACATGACCGAACGTTTCAGCCAGTCCCt                                                                                 >  6:148486/1‑103 (MQ=255)
                                                                                                              tttttGCCGGCCTCACTGCGAGGCCAGTCCATGTCAAGTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGt                                   >  5:56384/1‑149 (MQ=255)
                                                                                                                  tGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTc                                <  6:266078/148‑1 (MQ=255)
                                                                                                                    ccGGCCTCACTGTGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGt                              <  2:19872/148‑1 (MQ=255)
                                                                                                                    ccGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGt                              <  5:295944/148‑1 (MQ=255)
                                                                                                                         cTCACTCCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCTCCCTCATCAACGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCg                        <  8:229066/149‑1 (MQ=255)
                                                                                                                              tGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAgagcg                                                                   >  7:122694/1‑101 (MQ=255)
                                                                                                                              tGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAgagcg                                                                   <  8:122694/101‑1 (MQ=255)
                                                                                                                                cGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGATtgt                   <  7:239346/147‑1 (MQ=255)
                                                                                                                                  aGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGATTGTGCGa               >  2:228851/1‑149 (MQ=255)
                                                                                                                                       aGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGATTGTGCGACgg            >  1:183337/1‑147 (MQ=255)
                                                                                                                                         tacATGTCAATTGCCAAGGTCTCGCTCTCAACATAATCAAGGAAAGCAACATGACCGAACGTTTCAGAAAGTCACTGCGCCGCAAGAGAGAGACGAACTGATCGGCGGACGTTGGACACCCTTTCGTGACACGATAGTGAGACGGCAg         >  2:234822/3‑148 (MQ=255)
                                                                                                                                              gTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGATTGTGCGACGGCAgcctgcc   >  2:233007/1‑149 (MQ=255)
                                                                                                                                                cAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGATTGTGCGACGGCAgcctgcct  >  6:170379/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ACCAAGGGGAGCCCGGCAATGGGCTGAGAAACCGCTGAACGCTACAGCGCGGTGACCCTTCGAACCTGATCCGGATCATGCCGGCGAAGGGATGGGGCTTGTAACCTGCCTTTTTGCCGGCCTCACTGCGAGGCCAGTCCATGTCAATTGCCAAGGTCTCGCTGTCACCCTCATCAAGGACAGCAACATGACCGAACGTTTCAGCCAGTCCCTGCGCCGCCAGAGCGAGCCGACCTGGTCGGCGGCCGTTGGCCACCGTTTCGTGACCCGATTGTGCGACGGCAGCCTGCCT  >  NC_002947/3613941‑3614232

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: