Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,601,431 C→T 7.4% G154D (GGC→GAC)  PP_4073 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,601,4310CT 7.4% 153.2 / 2.3 54G154D (GGC→GAC) PP_4073hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base T (2/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

TCGTCGATAAGCGGGTCCAGCCCCTGAGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAACCTCGAGATGTTCTGCACGCAGCTGATGATCAGG  >  NC_002947/4601284‑4601571
                                                                                                                                                   |                                                                                                                                            
gcggcgATAAGCGGGTCCAGCCCCGGAGCGGCAAAAACGACGCCATCAACGGCAGGGGTGGCGCGCTCTACGGGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGcc                                                                                                                                             <  1:152910/145‑1 (MQ=255)
   tcgATAAGCGGGTCCAGCCCCTGAGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGtt                                                                                                                                          >  6:113016/1‑149 (MQ=255)
       tAAGCGGGTCCAGCCCCTGAGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTcc                                                                                                                                      <  8:83825/149‑1 (MQ=255)
       tAAGCGGGTCCAGCCCCTGAGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTcc                                                                                                                                      <  2:138125/149‑1 (MQ=255)
           cGGGTCCAGCCCCTGAGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCa                                                                                                                                  >  4:276636/1‑149 (MQ=255)
                ccAGCCCCTGAGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGc                                                                                                                             >  8:159075/1‑149 (MQ=255)
                 cAGCCCCTGAGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGc                                                                                                                             >  7:231228/1‑148 (MQ=255)
                     cccTGAGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTACTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACgg                                                                                                                        >  7:191329/1‑149 (MQ=255)
                       cTGAGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGgtg                                                                                                                      >  2:126293/1‑149 (MQ=255)
                          aGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGgtgt                                                                                                                     <  1:126293/147‑1 (MQ=255)
                          aGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGtgtg                                                                                                                    <  8:231228/148‑1 (MQ=255)
                          aGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTgg                                                                                                                   <  1:280412/149‑1 (MQ=255)
                          aGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTgg                                                                                                                   >  2:280412/1‑149 (MQ=255)
                           gCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTgg                                                                                                                   <  3:90188/148‑1 (MQ=255)
                                     cGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAgcg                                                                                                        >  5:129858/1‑149 (MQ=255)
                                        gTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGccc                                                                                                     >  1:42904/1‑149 (MQ=255)
                                        gTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGccc                                                                                                     <  1:191532/149‑1 (MQ=255)
                                         tCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGcccc                                                                                                    <  2:42904/149‑1 (MQ=255)
                                           cATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGc                                                                                                  >  6:30767/1‑149 (MQ=255)
                                            aTCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCt                                                                                                 >  3:188800/1‑149 (MQ=255)
                                             tCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTc                                                                                                <  5:142585/149‑1 (MQ=255)
                                               aaCGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTccc                                                                                              >  4:178500/1‑149 (MQ=255)
                                                       ggTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAAc                                                                                      <  4:188800/149‑1 (MQ=255)
                                                         ttGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCg                                                                                    >  5:15578/1‑149 (MQ=255)
                                                                tctctACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTg                                                                              >  3:238354/1‑148 (MQ=255)
                                                                    tACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTc                                                                         >  2:271348/1‑149 (MQ=255)
                                                                    tACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTc                                                                         >  4:99935/1‑149 (MQ=255)
                                                                          aaTGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTgg                                                                   >  8:182019/1‑149 (MQ=255)
                                                                          aaTGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTg                                                                    <  4:238354/148‑1 (MQ=255)
                                                                            tGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCa                                                                 >  6:158045/1‑149 (MQ=255)
                                                                              cGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCAt                                                                <  3:178500/148‑1 (MQ=255)
                                                                                aTGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGa                                                             >  5:231999/1‑149 (MQ=255)
                                                                                    gCCCCGATCAATTGCTTGATCTACTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACggg                                                         <  8:191329/149‑1 (MQ=255)
                                                                                     ccccGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGTCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGt                                                        >  5:161919/1‑149 (MQ=255)
                                                                                             aaTTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCgg                                                 <  6:15578/148‑1 (MQ=255)
                                                                                             aaTTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGt                                                <  5:113016/149‑1 (MQ=255)
                                                                                                 gCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGAc                                             <  5:158045/148‑1 (MQ=255)
                                                                                                 gCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGAc                                             <  7:205458/148‑1 (MQ=255)
                                                                                                  cTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGAcc                                            <  5:30767/148‑1 (MQ=255)
                                                                                                  cTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCa                                           <  3:99935/149‑1 (MQ=255)
                                                                                                         tGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCaa                                    >  3:208969/1‑149 (MQ=255)
                                                                                                           cTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGTCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAAc                                   >  3:284049/1‑148 (MQ=255)
                                                                                                           cTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGTCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAAc                                   <  4:284049/148‑1 (MQ=255)
                                                                                                           cTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAAcc                                  >  2:163380/1‑149 (MQ=255)
                                                                                                                tACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAACCTCgag                             <  4:208969/149‑1 (MQ=255)
                                                                                                                   ccGGCAGGCCTCTGCCGGCGTCGAGCAGTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAACCTCGAGATg                          >  5:266704/1‑149 (MQ=255)
                                                                                                                      gCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAACCTCGAGATGTTc                       <  3:276636/149‑1 (MQ=255)
                                                                                                                                 ccGGCGTCGAGCACTTGGTCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAACCTCGAGATGTTCTGCACGCAGCt            <  6:161919/149‑1 (MQ=255)
                                                                                                                                  cGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCACCCTCGAGATGGTCTGCACGCAGCt            >  2:211176/1‑148 (MQ=255)
                                                                                                                                    gCGTCGATAACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGTGTTGGGCGGTGACCAGTTGCAACCTCGAGATGTTCTGCACGCAGCtgat         >  5:97008/1‑149 (MQ=255)
                                                                                                                                    gCGTCGAGCACTTGGCCGTTTTCCCTCACTTGCCACTGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAACCTCGAGATGTTCTGCACGCAGCtgat         <  2:250870/149‑1 (MQ=255)
                                                                                                                                    gCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAACCTCGAGATGTTCTGCACGCAGCtgat         <  1:271348/149‑1 (MQ=255)
                                                                                                                                           gCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAACCTCGAGATGTTCTGCACGCAGCTGATGATCAgg  >  6:272680/1‑149 (MQ=255)
                                                                                                                                              cTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAACCTCga                              >  7:215325/1‑118 (MQ=255)
                                                                                                                                              cTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAACCTCga                              <  8:215325/118‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
TCGTCGATAAGCGGGTCCAGCCCCTGAGCGTCAAAAACGAGTCCATCAACGGCAGGGTTGTCGCTCTCTACGTGAATGCGATGGGCCCCGATCAATTGCTTGATCTGCTCGATACCCGGCAGGCCTCTGCCGGCGTCGAGCACTTGGCCGTTTTCCATCACTTGCCACGGTGTGGGTACAGTAGCGCCCCGCTCCCCTCTCAACCGCACCTGCACTCGGCTGGCATCGACGGGTTGGGCGGTGACCAGTTGCAACCTCGAGATGTTCTGCACGCAGCTGATGATCAGG  >  NC_002947/4601284‑4601571

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