Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,156,730 G→A 9.7% I243I (ATC→ATT)  PP_2772 ← monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,156,7300GA 9.7% 107.0 / 3.7 41I243I (ATC→ATT) PP_2772monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (21/16);  new base A (2/2);  total (23/18)
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.19e-01

GGCGCACGCCCTCCAGACCGACCCAGAACTCGCCGTAGTGGCGCAGGATGTAAGGGT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCTG  >  NC_002947/3156585‑3156876
                                                                                                                                                  |                                                                                                                                                  
ggCGCACGCCCTCCAGACCGACCCAGAACTCGCCGTAGTGGCGCAGGATGTAAGGGT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATa                                                                                                                                                 <  7:30495/149‑1 (MQ=255)
   gcACGCCCTCCAGACCGACCCAGAACTCGCCGTAGTGGCGCAGGATGTAAGGGT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTAATAcc                                                                                                                                               >  7:263240/1‑148 (MQ=255)
   gcACGCCCTCCAGACCGACCCAGAACTCGCCGTAGTGGCGCAGGATGTAAGGGT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTAATAcc                                                                                                                                               >  5:169647/1‑148 (MQ=255)
               gaccgaccCAGAACTCGCCGTAGTGGCGCAGGATGTAAGGGT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTg                                                                                                                                  <  3:284262/149‑1 (MQ=255)
                      ccAGAACTCGCCGTAGTGGCGCAGGATGTAAGGGT‑CTTCGCTGCCATGTCGGGCTGTGGAGCGGACCCACGGGCGGCCTTCCTGCAGGGTGTACTGGCGAGCCTCTGCCAGAGAGCCTTCGGTGATACCCAGGAACACGTTGGCGAAAt                                                                                                                           >  7:309539/1‑149 (MQ=255)
                        aGAACTCGCCGTAGTGGCGCAGGATGTAAGGGTGCTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATg                                                                                                                          >  1:50095/1‑149 (MQ=255)
                          aaCTCGCCGTAGTGGCGCAGGATGTAAGGGT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCaa                                                                                                                       <  1:200523/149‑1 (MQ=255)
                          aaCTCGCCGTAGTGGAGCAGGATGTAAGGGT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCTAAATGCaa                                                                                                                       <  8:192830/149‑1 (MQ=255)
                          aaCTCGCCGTAGGGGCGCAGGATGTAAGGGT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCNAAATGCaa                                                                                                                       <  8:198137/149‑1 (MQ=255)
                            cTCGCCGTAGTGGCGCAGGATGTAAGGGT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACt                                                                                                                     >  6:14630/1‑149 (MQ=255)
                                ccGTAGTGGCGCAGGATGTAAGGGT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTgcgc                                                                                                                 >  8:258935/1‑149 (MQ=255)
                                                tGTAAGGGT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAgg                                                                                                 <  6:310488/149‑1 (MQ=255)
                                                    aGGGT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTg                                                                                              <  2:362927/148‑1 (MQ=255)
                                                     gggT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGc                                                                                            >  7:292051/1‑149 (MQ=255)
                                                       gT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCGTACTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCg                                                                                           >  1:78053/1‑148 (MQ=255)
                                                          cTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAg                                                                                        >  3:186445/1‑149 (MQ=255)
                                                           ttCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCATCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCTGGCGCAGGCTGGCGAAgg                                                                                       >  5:98404/1‑149 (MQ=255)
                                                             cGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAggg                                                                                      <  2:50095/148‑1 (MQ=255)
                                                             cGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAgg                                                                                       <  8:292051/147‑1 (MQ=255)
                                                                       tCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCg                                                                           <  5:14630/149‑1 (MQ=255)
                                                                       tCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCg                                                                           >  7:90284/1‑149 (MQ=255)
                                                                              gtgGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCggg                                                                    >  7:273300/1‑149 (MQ=255)
                                                                                   gCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTc                                                                 >  1:6216/1‑147 (MQ=255)
                                                                                          cACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCAt                                                        >  6:149921/1‑149 (MQ=255)
                                                                                            cGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCg                                                                      >  6:24105/1‑133 (MQ=255)
                                                                                            cGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCg                                                                      <  5:24105/133‑1 (MQ=255)
                                                                                                    cttccttcAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGc                                                                            >  7:120037/1‑119 (MQ=255)
                                                                                                    cttccttcAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGc                                                                            <  8:120037/119‑1 (MQ=255)
                                                                                                         ttcAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGACGGGGTCGAGAAGCATCTAGTTCTCTTCGac                                         >  8:208261/1‑149 (MQ=255)
                                                                                                                tgtAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGGACACGTTGGCGAAATGCAACATCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTACGCGGGCCGGGGTCGAGAAGCATATCCTTCTCTTCGacgcgcacg                                  >  8:222015/1‑149 (MQ=255)
                                                                                                                tgtAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTtctc                                               <  5:162238/136‑1 (MQ=255)
                                                                                                                tgtAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTtctc                                               >  6:162238/1‑136 (MQ=255)
                                                                                                                  tAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGacgcgcacgcg                                >  4:135953/1‑149 (MQ=255)
                                                                                                                       gCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCAtctc                                                     <  5:273717/123‑1 (MQ=255)
                                                                                                                       gCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCAtctc                                                     >  6:273717/1‑123 (MQ=255)
                                                                                                                         gAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGaaa                         <  7:222015/149‑1 (MQ=255)
                                                                                                                             ctctTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTcg                     >  7:104497/1‑149 (MQ=255)
                                                                                                                                     gagaGCCTTCGGTAATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCg             <  6:169647/149‑1 (MQ=255)
                                                                                                                                     gagCGCCTTCGGTAATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCg             <  8:263240/149‑1 (MQ=255)
                                                                                                                                            ttCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCgg      >  8:203726/1‑149 (MQ=255)
                                                                                                                                                gTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCTg  <  2:78053/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGCGCACGCCCTCCAGACCGACCCAGAACTCGCCGTAGTGGCGCAGGATGTAAGGGT‑CTTCGCTGACATGTCGGGCTGTGGAGCGGAACCACGGGCGGCCTTCCTTCAGGGTGTAATGGCGAGCCTCTTCCAGAGCGCCTTCGGTGATACCCAGGAACACGTTGGCGAAATGCAACTGCGCGATCAGCGGGCGCAGGCTGGCGAAGGGGGTGCTCAGCGGGCCGGGGTCGAGAAGCATCTCGTTCTCTTCGACGCGCACGCGCTCGAAAGTCGCGCTGCCGCTGTCGGTCTG  >  NC_002947/3156585‑3156876

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