Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,484,946 G→A 7.6% L26L (CTG→CTA)  PP_3977 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,484,9460GA 7.6% 149.1 / 2.9 54L26L (CTG→CTA) PP_3977hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/25);  new base A (2/2);  total (27/27)
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

GTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATCCAGAGCCTTCCCCAGCTGAATGATGTCACTTCCGCCTTCGCGTTGACCTCCGCCG  >  NC_002947/4484797‑4485094
                                                                                                                                                     |                                                                                                                                                    
gTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGGCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATActg                                                                                                                                                      >  1:224584/1‑150 (MQ=255)
      aTTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGt                                                                                                                                                <  2:2013211/150‑1 (MQ=255)
          aaaCCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCACAATCTTCAGGATTTCCGCAAGCGTGCAGTCGACTCCATACTGCTGAGTGCCa                                                                                                                                            >  2:655194/1‑150 (MQ=255)
           aaCCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGCCGACTGCATACTGCTGAGTGCCAg                                                                                                                                           <  1:1584914/150‑1 (MQ=255)
                      tAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAgatcat                                                                                                                                <  2:889818/150‑6 (MQ=255)
                               acctacAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAgg                                                                                                                       <  2:1270421/150‑1 (MQ=255)
                                   acAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTaa                                                                                                                   >  1:831497/1‑150 (MQ=255)
                                    cAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTaaa                                                                                                                  <  1:1342955/150‑1 (MQ=255)
                                         ttAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGCGTAAatggat                                                                                                             <  2:1160683/150‑1 (MQ=255)
                                                  gATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGa                                                                                                                                   >  1:811339/1‑119 (MQ=255)
                                                  gATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGa                                                                                                                                   <  2:811339/119‑1 (MQ=255)
                                                     cGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCAcgc                                                                                                 <  1:1753106/150‑1 (MQ=255)
                                                     cGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTATTCCGCGCAAGCTTCAGGATTTCCACAAGCGTGGAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAGAGGGTAAATGGATTGCAGATCAcgc                                                                                                 >  2:503400/1‑150 (MQ=255)
                                                       gTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCAcgc                                                                                                 <  1:669307/148‑1 (MQ=255)
                                                       gTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCAcgc                                                                                                 >  2:669307/1‑148 (MQ=255)
                                                            tCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTACTGAGTGc                                                                                                                                              >  1:1902590/1‑98 (MQ=255)
                                                            tCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTACTGAGTGc                                                                                                                                              <  2:1902590/98‑1 (MQ=255)
                                                            tCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTACTGAGTGc                                                                                                                                              <  2:1902622/98‑1 (MQ=255)
                                                            tCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTACTGAGTGc                                                                                                                                              >  1:1902622/1‑98 (MQ=255)
                                                              gACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGT‑AATGGATGGCAGAGCACGCGAGAGAAtgc                                                                                       >  1:168165/1‑150 (MQ=255)
                                                               aCGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAAtgc                                                                                       <  2:1411721/150‑1 (MQ=255)
                                                                cGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATActgct                                                                                                                                                    <  2:468027/88‑1 (MQ=255)
                                                                cGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATActgct                                                                                                                                                    >  1:468027/1‑88 (MQ=255)
                                                                    cTTTATGGACTGTAAGCTAGCAAAACCGGGTATTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGAAAATGGATGGCAGAGCACGCGAGAGAATGCTGCaa                                                                                  >  1:1022353/1‑150 (MQ=255)
                                                                       tATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTGTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGcc                                                                               >  1:511074/1‑150 (MQ=255)
                                                                        aTGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCt                                                                              >  1:194880/1‑150 (MQ=255)
                                                                         tGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCCGAGCACGCGAGAGAATGCTGCAAGCCTc                                                                             <  1:2066440/150‑1 (MQ=255)
                                                                          ggACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGTAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCa                                                                            >  1:1122596/1‑150 (MQ=255)
                                                                              tGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAg                                                                                                                                           >  2:1500253/1‑83 (MQ=255)
                                                                              tGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAg                                                                                                                                           <  1:1500253/83‑1 (MQ=255)
                                                                                   gCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGc                                                                   >  1:2020483/1‑150 (MQ=255)
                                                                                        cAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGTAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACa                                                              <  1:503400/150‑1 (MQ=255)
                                                                                             ccGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGc                                                                   <  2:1281553/140‑1 (MQ=255)
                                                                                             ccGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGc                                                                   >  1:1281553/1‑140 (MQ=255)
                                                                                             ccGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATc                                                         <  2:831497/150‑1 (MQ=255)
                                                                                               gggTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATCCa                                                       <  2:2102056/150‑1 (MQ=255)
                                                                                               gggTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTACACAACGTGTAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGCATCCa                                                       >  2:572822/1‑150 (MQ=255)
                                                                                                ggTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGT‑AATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATCCaga                                                     <  2:168165/150‑1 (MQ=255)
                                                                                                    ttCCGCGCAAGCTTCAGGATTTCCGCAATCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATCCAGAGcc                                                  >  1:315909/1‑150 (MQ=255)
                                                                                                       cgcgCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATCCAGAGCCTTc                                               <  2:224584/150‑1 (MQ=255)
                                                                                                        gcgcAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATCCAGAGCCTTc                                               <  2:2241479/149‑1 (MQ=255)
                                                                                                        gcgcAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATCCAGAGCCTTc                                               >  1:2241479/1‑149 (MQ=255)
                                                                                                              gCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATc                                                         >  2:1327524/1‑133 (MQ=255)
                                                                                                              gCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATc                                                         <  1:1327524/133‑1 (MQ=255)
                                                                                                                ttCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATCCAGAGCCTTCCCCAGCTGa                                      >  1:1842262/1‑150 (MQ=255)
                                                                                                                 tCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCAcgc                                                                                                 >  2:765293/1‑90 (MQ=255)
                                                                                                                 tCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCAcgc                                                                                                 <  1:765293/90‑1 (MQ=255)
                                                                                                                     gATTTCCGCAAGCGTGCAGTCGACTGCATACTCCCGTGTGCCTGGCCAAGGAGCGTTTTAAAGGGTAAATGGATGGCAGAGTACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATCCAGAGCCTTCCCCAGCTGAatgat                                 >  2:1154352/1‑150 (MQ=255)
                                                                                                                                       gTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATCCAGAGCCTTCCCCAGCTGAATGATGTCACTTCCGCCTTcgcg               >  1:353040/1‑150 (MQ=255)
                                                                                                                                       gTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATCCAGAGCCTTCCCCAGCTGAATGATGTCACTTCCGCCTTcgcg               >  1:2213008/1‑150 (MQ=255)
                                                                                                                                        tCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACaa                                                                     >  1:1857806/1‑95 (MQ=255)
                                                                                                                                        tCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACaa                                                                     <  2:1857806/95‑1 (MQ=255)
                                                                                                                                         cGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATCCAGAGCCTTCCCCAGCTGAATGATGTCACTTCCGCCTTCGCGtt             <  1:2249890/150‑1 (MQ=255)
                                                                                                                                               cATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGaa                                                                <  1:1453952/93‑1 (MQ=255)
                                                                                                                                               cATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGaa                                                                >  2:1453952/1‑93 (MQ=255)
                                                                                                                                                  aCTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGc                                                                   <  2:724351/87‑1 (MQ=255)
                                                                                                                                                  aCTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGc                                                                   >  1:724351/1‑87 (MQ=255)
                                                                                                                                                    tgctcAGTGCCAGGACAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATACAGAGCCTTCCCCAGCTGAATGATGTAACTTCCGCCTTCGCGTTTACCTcagccg  >  2:145649/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTAGAGATTCAAACCCTTATCCTAGGCTCCTACCTACAGACTTAATTTCTGATCGGTAGCTCGACGTCCTTTATGGACTGTAAGCTAGCAAAACCGGGTCTTCCGCGCAAGCTTCAGGATTTCCGCAAGCGTGCAGTCGACTGCATACTGCTGAGTGCCAGGCCAAGGACAGTTTTAAAGGGTAAATGGATGGCAGAGCACGCGAGAGAATGCTGCAAGCCTCAGAAACAAGCGAACAGGATCCAGAGCCTTCCCCAGCTGAATGATGTCACTTCCGCCTTCGCGTTGACCTCCGCCG  >  NC_002947/4484797‑4485094

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑